OKAMOTO'S SPACE FOR THE FIRST PAINLEVE EQUATION 
IN BOUTROUX COORDINATES 



J.J. DUISTERMAAT AND N. JOSHI 



Abstract. We study the completeness and connectedness of asymptotic behaviours of solutions of the first 
Painleve equation d 2 y/ dx 2 = 6y 2 + x, in the limit x — > oo, x £ C. This problem arises in various physical 
contexts including the critical behaviour near gradient catastrophe for the focusing nonlinear Schrodinger 
equation. We prove that the complex limit set of solutions is non-empty, compact and invariant under the 
flow of the limiting autonomous Hamiltonian system, that the infinity set of the vector field is a repellor for 
the dynamics and obtain new proofs for solutions near the equilibrium points of the autonomous flow. The 
results rely on a realization of Okamoto's space, i.e., the space of initial values compactified and regularized 
by embedding in CP2 through an explicit construction of nine blow-ups. 



1. Introduction 



d y a 2 , n -n 

-6y z + x, (1.1) 



In this paper, we consider the completeness and connectedness of the asymptotic behaviours of 
the first Painleve equation 

i 2 : 

dx 2 

in the limit x — > oo, x G C. The first Painleve equation arises in many physical contexts, as a 
reduction of the Korteweg-de Vries equation, in the double scaling limit of random matrix models 
and in the critical behaviour near the point of "gradient catastrophe" of the solution to the Cauchy 
problem for the focusing nonlinear Schrodinger equation [9]. 



The asymptotic limit x — > oo of Equation (1.1) was first studied in 1913 by Boutroux j2j, who 



provided a transformation of variables that makes the asymptotic behaviours explicit. It is known 



that all solutions of (1.1) are meromorphic in C with double movable poles, i.e., with locations 
that change with initial conditions. Boutroux found that locally, in each patch near infinity, the 
general solutions are given to leading-order by elliptic functions. More detailed results about how 
the local asymptotic behaviours of solutions change slowly as x moves near infinity were provided 
by Joshi and Kruskal [151 [16], who constructed a complex multiple-scales method to carry out 
asymptotic analysis along a large circle in the complex plane for the first and second Painleve 
equations. Such local behaviours were used in the Riemann-Hilbert method, which was applied to 
deduce connections between behaviours valid along special directions approaching infinity (see the 
review by Kitaev |18j). 
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In addition to the two-parameter solutions asymptotic to elliptic-function behaviours, Boutroux 
identified five one-parameter family of solutions asymptotic to algebraic power expansions in certain 
sectors of angle 4tt/5 in C. He called them tronquee or truncated solutions. In each family of 
tronquee solutions, there is a unique solution whose algebraic expansion is valid in a sector of angle 
8tt/5. Boutroux called these tritronquee or triply truncated solutions. In the literature, this term 
has come to be associated with the unique tritronquee solution that is real on the real line; each of 



the other four such solutions can be obtained from this one by a discrete symmetry of Equation (1.1) 
corresponding to rotating variables in C. This real tritronquee solution appears as a distinguished 
solution in various physical problems (see, e.g., Dubrovin et al. [9]). In the form Yu = 6Y 2 — t 
(for y(x) = y (—£)), Joshi and Kitaev [14J constructed a sequence of solutions that converge to the 
tritronquee solution on the whole positive real axis and proved that the tritronquee solution has no 
poles whatsoever on the positive real axis. O. and R. Costin [7j applied Borel-summation methods 
to deduce many complex properties of the tronquee solutions. 



However, while it is known that the solution space of Equation (1.1) is connected, through 
Okamoto's [21] compactification and regularization of the space of initial- values, there has been no 
investigation (to our knowledge) of the completeness or connectedness of the known asymptotic 
behaviours of the solutions of this equation (or of any of the six Painleve equations). We tackle this 
problem by undertaking asymptotic analysis in Okamoto's space. Our approach relies on explicit 
resolution of singularities in an asymptotic version of this space. We note that although we focus 
on the first Painleve equation in this paper, our approach can be extended to the other Painleve 
equations. 



1.1. Outline of the paper. The paper is organized as follows. In we recast Equation (1.1) 



as a Hamiltonian system, provide a rescaling of it under Boutroux's transformation of variables, 
and summarize known properties of solutions. The resolution of singularities of this Boutroux- 
Painleve system is provided explicitly in Appendix A where we carry out the sequence of changes 
of variables necessary first to compactify and then blow up the nine base points of the system in 
CP2. The construction shows that the vector field is infinite on the union / := L)f =0 Lj of nine 
complex projective lines. 

We show in £[3] that the Boutroux-Painleve vector field is regular and transverse to the last 
complex projective line, Lg, in Sg. It is shown here that the Taylor expansion of the flow around a 
point on this line provides us with the Laurent expansion of the solutions y(x) near a pole. In £j4j we 
consider the vector field near the infinity set I and show that this is a repellor for the flow. We also 
construct the limit set for each solution and show that it is a non-empty, compact and connected 
subset of Sg, which remains invariant under the autonomous flow. As a corollary, we prove that 



every solution of Equation ( 1.1 ) must have an infinite number of poles in the complex plane. Finally, 
in we consider the Boutroux-Painleve system near the equilibria of the autonomous limit system 
and prove several results about tronquee solutions, ending with a determination of their sequence 
of poles near the boundaries of pole-free sectors, by using classical methods. 

1.2. Notation. The Painleve flow is denoted (2/1,2/2)) where y = y±, dy/dx = y 2 - in Boutroux's 
coordinates, this flow is transformed to (ui,u 2 ), where u\{z) = x~ 1 / 2 yi(x), u 2 (z) = x~ 3//4 2/2(2), 
with (4/5) a; 5 / 4 = z. We embed the resulting Painleve-Boutroux flow into CP2 by identifying the 
affine coordinates (1*1,1*2) with homogeneous coordinates as follows 

[1 : i*i : u 2 ] = [u{ 1 : 1 : u^ 1 u 2 ] =: [«02i = 1 : ^022] 

= [1*2 1 : u\ U2 1 : 1] =: [uo31 : 1*032 : ^ 

The line at infinity given respectively by U021 = 0, 1*031 = is denoted by Lq. For < i < 8, 
corresponding to the *-th stage of the regularizing (blow-up) sequence, we denote a base point by 
bi, the exceptional line attached to that base point by Lj+i and the coordinates in the j-th chart 



of the i-th blowup by (itijij tty^)- I n each coordinate chart, the Jacobian of the coordinate change 
from (u\,U2) to (1^1,1^2) is denoted by 

_ dujji dui j2 _ dujji dui j2 
du\ du2 du2 du\ 

The last space Sg constructed by this sequence of blow-ups is Okamoto's "space of initial values." 
The regularization of the limiting autonomous system obtained as z — > 00 is the same as that of 
the related Weierstrass cubic curve. We add the superscript ell to the above notation for the base 
points and lines and a superscript to the vector field to refer to the corresponding objects in the 
resulting anti-canonical pencil. 

2. BOUTROUX SCALING 



.2 , n o/ K ^-l.._ (2-2) 



Equation (1.1) can be viewed, upon the substitutions y = y%, dy/dx = y2, as a Hamiltonian 
system dy±/ d x = dH/dy2, dyil dx = —dH/dyi with an x-dependent Hamiltonian function 

H = H(x, yi, y 2 ) := y 2 2 /2 -2 Vl 3 -x Vl . (2.1) 

The function H is a weighted homogeneous polynomial in the sense that if we substitute x = A 4 £, 
yi = X 2 ui, and 2/2 = A 3 U2, then H = A 6 (ii 2 2 /2 — 2 ui 3 — £ u±). We have £ = 1 if and only if x = A 4 , 
when A = x 1 / 4 , y\ = x 1 / 2 -^, u\ = x _1 / 2 yi, 2/2 = x 3 / 4 U2, and 112 = x _3 / 4 y2- If x = x(z) and a 
dot means differentiation with respect to z, then iii = x (-(l/2)x- 1 ui + x 1 / 4 u 2 ). If we choose 
(4/5) x 5 / 4 = z then xx 1 / 4 = 1, and the Painleve system takes the form 

iii = U2 — 2 (5 z)^ 1 ui, 
u 2 = 6m 2 + 1 - 3(5z)- 1 n 2 . 

This is an order z _1 perturbation of the Hamiltonian system with Hamiltonian function equal to 
the ^-independent energy function 

E := u 2 2 /2 - 2 m 3 - m, 
where (2.3) 
E = (5Z)- 1 (2ui + 12ui 3 - 3u 2 2 ) = -(5 z)" 1 (6 E + 4 Ul ). 

It implies Boutroux's second order differential equation 

ui = Qui 2 + 1 - z^iii + 4(5z)~ 2 m (2.4) 

for u\. These transformations have been used by Boutroux |2] in order to investigate the asymptotic 
behavior of the solutions of the Painleve equation when x — > 00. 

Because f°° z~ x dz = 00, we cannot straightforwardly conclude that solutions of the Boutroux 



system (2.2) converge to solutions of the autonomous Hamiltonian system when z — > 00. Actually 



they don't: we will see that each solution of (2.2 ) converges to different solutions of the autonomous 



limit system, depending on the path along which z runs to infinity. 

Remark 2.1. If a and b are n onze ro complex constant complex numbers, then the substitutions 



y{x) = ar)(£), £ = bx turn (1.1) into the differential equation d 2 7?/d£ 2 = arf + where 
a = 6a/b 2 and j3 = 1/ab 3 , or equivalently b = (6/a/3) 1//5 and a = ab 2 /6. 

Boutroux [2, p. 311] took the first Painleve equation with the constants a = 6 and f3 = —6, 
or equivalently b = (— 1/6) 1 / 5 and a = b 2 , as his point of departure, and applied the substitutions 
X = (4/5) £ 5 / 4 and rj = ^/ 2 Y in order to arrive at the differential equation d 2 Y j dX 2 = 6 Y® — 
6 — (1/X) dY/ dX + (4/(5 X) 2 ) Y. Therefore the translation from Boutroux's notation to ours is 
X = (-1/6)VS, Y = (-l/6)-V2 u= (-l/6)-V2 Ul) and y< = (_i/ 6 )-3/4 ^ = (_i/ 6 )-3/4 ^ _ 

2 (5 z)^ 1 tii). An expression which plays a central role in Boutroux [21 §7-11] is 

(Y') 2 - 4Y 3 + 12 Y = (-l/6)" 3/2 (u 2 - 4u 3 -2u). 
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Joshi and Kruskal |15j . |16j took the first Painleve equation with the constants a = 3/2 and 
P = —3/2, or equivalently b = (— 8/3) 1 / 5 and a = b 2 /4, as their point of departure, and applied 
the substitutions Z = (4/5) £ 5 / 4 and r/ = ^^U, where they actually wrote z and u instead of 
Z and U, respectively. Therefore the translation from their z and u to ours is Z = (— 8/3) 1//4 z, 
U = 4 (-8/3)" 1 / 2 u= ± i y/6 ui, and dU / dZ = 4 (-8/3)" 3 / 4 u = 4 (-8/3)" 4 / 3 (u 2 - 2 (5 z)- 1 ui). 
A central role is played in [15], [16] by the function 



£ := ((dU/ dZf -U 3 + 3U)/2 = 2^ 2 (-3) 3/2 (u 2 /2 - 2 < 



,3 



U 



The functions (Y 7 ) 2 — 4 1" 3 + 12 Y and 5 are closely related to the energy function E in ( 2.3 ), as 
n 2 /2 -2u 3 -u = E-2(5 z)- 1 mu 2 + 2 (5 z)- 2 Ul 2 . (2.5) 



Remark 2.2. The Boutroux substitutions x = ((5/4) z) with inverse z = (4/5) x 5//4 , and y(x) = 

x 1 / 2 u(z) = x 1 / 2 u((4/5) x 5 / 4 ) = ((5/4) z) 2 ^ 5 u(z) with inverse ■u(z) = x^ 1 / 2 y{x) are singular at 
a; = and correspondingly z = 0. These substitutions introduce multi-valuedness of the solutions 
u(z) of the Boutroux-Painleve equation (|2.4[) when z runs around the origin in the complex plane, 



where the solutions y(x) of the Painleve equation (1.1) are single-valued. 



More precisely, every local solution y(x) of ( 1.1 ) extends to a single- valued meromorphic function 



on the whole complex x-plane, where the poles are of order two and have leading coefficient equal 
to 1. This is the Painleve property in its strongest form; see \10\ Remark 1.1] for some remarks on 
its proofs in the literature. 

The equation u(z) = ((5/4) z)~ 2 ^ y(((5/4) z) 4 ^ 5 ), in combination with the single-valuedness of 
y(x), implies that the analytic continuation of u(z) along the path z e 1 e , avoiding the poles of u(z), 
returns to its opposite if 9 £ M runs from to 5/4 times 2ir. This may be expressed by the formula 

u l {ze^ i ' 2 )= - Ul (z), u 2 (ze^'l 2 ) = iu 2 (z), (2.6) 



where the second equation follows from the first, in view of the first equation in (2.2). This 
observation has been used in Joshi and Kruskal |15| Sec. 5] as a consistency check for their 
asymptotic results for u(z) for large \z\. 



Remark 2.3. Each solution y(x) of the Painleve equation (1.1) has a convergent Laurent expansion 

oo 

y( x ) = ynx n 



n=no 



for < |x| << 1, when u(z) = ((5/4) z) 2 ^ y(x) = ((5/4) z) 2//5 y(((5/4) z) 4 / 5 ) implies the conver- 
gent power series 

oo 

«(*)=£ y n (m) Z ) ( - 2+4n)/5 

n=no 

for < | jzr| << 1. We have the following cases. 

i) y(0) is finite, when no = 0. The Painleve equation y" = 6y 2 +x is equivalent to the recursive 
equations 

n-2 

n(n-l)y n = 6 y n -2~m Vm for n > 2, n / 3, (2.7) 
V3 = 2 y(m + l/6 (2.8) 



for the coefficients y n . The mapping which assigns to the solution y(x) the complex numbers 
2/0 = 2/(0) and y\ = y'(0) is bijective from the set of all regular solutions y(x) near x = 
onto C 2 . Subcases: 

ia) 2/0 = 2/1 = 0) when \2.7\ for n = 2 yields that 2/2 = 0, whereas (2.8) implies that ys = 1/6. 
An induction on n yields that y n = unless n £ 3 + 5Z, as n- 2 — m = 3 + 5 k and 
m = 3 + 5 / imply that n = 3 + 5 (fc + / + 1). Because -2 + 4(3 + 5 j) = 10(2 j + 1), 
it follows that u(z) = X]^=o 2/3+5 j ((5/4) z) 2 ( 2j+1 ). In particular this solution u(z) is 
single- valued. 

ib) t/o = and y x ± 0, when u{z) = y\ ((5/4) z) 2/5 + 0(z 2 ) as z — > 0. 

ic) y / 0, when u(z) = y ((5/4) z)~ 2 / 5 + 0(z 2 / 5 ) as z -> 0. 

y(x) has a pole at x = 0, when no = —2, y_2 = 1, y-i = 2/0 = 2/1 = 2/2 = 0, and 2/3 = —1/6. 
The mapping which assigns to the solution y(x) the coefficient y% is bijective from the set 
of solutions with a pole at x = onto C, see for instance |10|. the text following (11.3)]. For 
n > 5 the Painleve equation y" = 6 y 2 + x implies the recursive equations 



(n(n-l)-12)i/ n = 6 ^ 



(2.9) 



m=3 



Subcases: 

iia) t/4 = 0. As in ia) it follows from (2.9) by induction on n that y n = unless n G 3 + 5Z, 
and it follows that n(z) = X^=-i 2/3+5j ((5/4) z) 2 ( 2:,+1 ). In particular this solution 
ii(z) is single-valued. 

iib) y 4 / 0, when «(«) = ((5/4) z)- 2 - (1/6) ((5/4) z) 2 + y 4 ((5/4) z) 14 / 5 + 0(z 18 / 5 ) as z -»• 0. 



It follows that the solution u(z) of the Boutroux-Painleve equation (2.4) is not single- valued, unless 
we are in the cases ia) or iia). That is, the solution y(x) of (1.1) is either equal to the unique 
regular solution near x = for which y(0) = y'(0) = 0, or y(x) is the unique solution of (1.1) with 
a pole at x = such that m = 0. 



If y(x) is a solution of ( 1.1 ), and a G C is a fifth root of unity, that is, a = 1, then i + a y(a 



is also a solution of (1.1). The solutions in ia) and iia), corresponding to the single- valued 
solutions u(z) of (2.4), are exactly the solutions which are invariant under this five- fold symmetry, 
that is, which satisfy y(x) = a -1 y(a 2 x) for every fifth root of unity a, as this means that in the 
Laurent expansion of y(x) only the powers x 3 appear such that 2 j — 1 G 5Z 2j — l G 5 (2 Z + 1) 
44> j G 5 Z + 3. The solutions y(ic) in ia) and iia) appear in Boutroux [2 p. 336, 337]. 

Remark 2.4. If y(x) is a solution of the Painleve equation (|1.1[ ), and a G C is a fifth root of unity, 



that is, 



1, then x i — ^ a 



1 y(a 2 x 



is also a solution of (ILlj). 



The solutions y{x) in ia) and iia) of Remark 2.3, corresponding to the single- valued solutions 



u(z) of the Boutroux-Painleve equation (2.4), are exactly the solutions which are invariant under 
this five-fold symmetry, that is which satisfy y{x) = a -1 y(a 2 x) for every fifth root of unity a, as 
this means that in the Laurent expansion of y(x) only the powers x 3 appear such that 2j — l G 5Z 
O 2 j - 1 G 5 (2Z + 1) o j G 5Z + 3. This explains why Boutroux p. 336, 337] called the 



solutions y(x) in ia) and iia) of Remark 2.3 the symmetric solutions. 



3. The poles in Okamoto's space 

In this section, we consider the Boutroux-Painleve vector field in Okamoto's space S$ constructed 
explicitly in Appendix A. This construction shows that the vector field has no base points in Sg, 
is infinite along the configuration I := (Ji=o ^ °^ n i ne complex projective lines, and regular in 
Sg \ I. For this reason we will call the set I the infinity set of the vector field. 
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The change of coordinates to the coordinate chart (ugn, 1*912) is shown in Appendix A to be 
u 911 (z) = ui(z) -9 u 2 {z) (-32 Ul {z) 7 u 2 (z) - 4 Ul (zf u 2 {zf + u 2 (z) 7 + 256 (5 z)' 1 u^z) 8 ), 
u 9 \ 2 (z) = u 1 (z)u 2 (z)~ , 
ui(z) = u 9l2 (z)- 2 (4 + 32 u 912 {z) 4 + u 911 (z)u 9l2 {zf - 256(5 z)' 1 u 912 {zf)- 1 (3.1) 
u 2 {z) = ug 12 (z)~ 3 (4 + 32u 9 i 2 (z) 4 + Ugu(z) u 9X2 {zf - 256 (5 z)' 1 u 9l2 {zf)~ l 

The set of points in Sg\I which project to Lq, the set where («i, u 2 ) is infinite, is equal to L 9 \I. 
Because L 9 nl = L 9 PiL^ consists of one point, L 9 \I is isomorphic to the affine complex plane. The 
regular vector field in S 9 \ / is nonzero at and transversal to L 9 \I. A solution crosses the complex 

line L 9 \ at the time z = £, if and only if u(z) = u\{z) becomes infinite as z — > C- The whole 
set L 9 \ I is visible in the coordinate chart (itgn, W912), where it is the line M912 = 0, parametrized 



by M911 G C. Because u{z) = u\(z) given by Equation (3.1) is a rational expression in z, ugn(z), 
and 14922(20, and the solution z \-t (u 9 n(z), ug\ 2 (z)) of the regular non-autonomous system is a 
complex analytic function in a neighborhood of z = £ with ugi 2 (() = and a := ugn(C) £ C, it 
follows that the solution z \-t u(z) of the Boutroux-Painleve equation is a meromorphic function in 
a neighborhood of z = C, with a pole of order two. For this reason the line L 9 \ I is called the pole 
line. 



Consider the vector field (ttgn, ugi 2 ) which is given in Section A. 10 of Appendix A. We recall it 
here for ease of reference by the reader: 

ug U = (4: + 32ug 12 4 + u 9 nUg 12 fi - 256(5z)- 1 u 912 5 )- 1 

x [u 912 (-2 11 - 2 6 • 5 u 9U u 912 2 + 2 13 • 7 n 9 i 2 4 - 3 2 n 911 2 u 9l2 A 

+2 12 U911 n 912 6 + 2 16 • 3 U912 8 + 2 3 • 3 2 u 911 2 u 912 8 

+2 12 • 5 ug U u 912 W + 2 6 • 11 u 9ll 2 u 9 i2 12 + 2 3 n 911 3 u 912 lA ) 

-2 (5 z)~ l (2 2 • 3«9n - 2 12 • 3 2 u 9 i2 2 - 2 5 • 3 2 • 7u 9n u 912 A 

+2 15 -3-5 u 912 e + 3 n 9 n 2 u 912 6 + 2 10 • 17 ugii u 912 8 

+2 17 • 19 u 912 w + 2 13 • 3 • 7 ngii n 9 i2 12 + 2 7 • 23 u 911 2 u 912 u ) 

+2 9 (5 z)~ 2 u 9 i2 3 (-2 6 -3-5 + 3 ug U u 912 2 

+2 13 u 912 4 + 2 14 • 5 n 9 i2 8 + 2 8 • 11 U911 n 912 10 ) 

-2 24 - 7 (5zrW 2 ], 

n 9 i2 = -(4 + 32u9i2 4 + n 9 iiU9i2 6 - 256 (5z)- 1 M9i2 5 ) _1 

x [2 - 2 4 n 912 4 - u 911 n 912 6 + 2 s u 912 8 + 2 3 u 9U u 912 10 

+2 10 U 9 i2 12 + 2 6 U 911 U 9i2 14 + Ug U 2 Ug 12 16 

-(5 z)- 1 u 912 (2 2 - 2 5 • 7 n 9 i2 4 + u 911 u 912 6 
+2 11 u 912 8 + 2 14 n 912 12 + 2 9 ug Xl n 912 14 ) 
+2 8 (5z)- 2 n 912 6 (l + 2 8 U912 8 )]. 
It follows from the equation for u 9 \ 2 that the coefficie nts of (z — £)* in the Taylor expansion at 



z = C of 1(912(2) do not depend on a for 1 < i < 6, when (3.1 ) shows that the coefficients of {z — C,) 3 
in the Laurent expansion at z = C of u(z) do not depend on a for —2 < j < 3. Substitution of the 
Taylor expansion at z = C of order % of ug\ 2 {z) in the formula for ?i 9 i2 yields the Taylor expansion 
at z = C of order i of itgi 2 (z) of order i, hence the Taylor expansion at z = ( of order i+ 1 of tig^z), 
as long as i < 5. Then substitution of Mgn(C) = a in the formula for ugi2 yields the coefficient of 
(z — C) 6 in the Taylor e xpa nsion at z = £ of -ug^-z), hence of (z — () 7 in the Taylor expansion at 



z = ( of ugi 2 (z), when (3.1) yields the coefficients of (z — in the Laurent expansion at z = ( of 



u(z) for -2 < j < 4. This yields 
u 9n (z) = a + 0(z-(), 

u ^ z ) = -\ {z ~ c) ~ {z ~ ° 2 + 2 2 • 5^ • c 2 {z ~ ° 3 ~ 2 3 c 3 {z ~ ° 4 

/ 3-7-19 J_\ 5 / 3-7-19 41 \ 

+ V2 3 -5 5 -C 4 2-5^ 1 U V2-5 6 -C 5 2 2 -3-5 2 -Cy 1 U 
3-19-29 41 3o \ , an7 ^„ anSx 



2-5 7 -C 6 2-3-5 3 -C 2 2 9 -7 

1 3 31 

«M = (*-0^- r? (.-0- 1 + 3rr? - rF ^(.-0 

/ 19-283 _ J_\ _ 2 _ / 3- 11 -727 11 , 



2 4 -5 5 -C 4 2-57 v s/ V 2 4 -5 6 -C 5 2-4-5 2 -C 
/ 197-443 29 a \ , .. 4 _„ ^ , n nS 

+ (2675^ + 2 3. 3 - 52 -C 2 " 2BT7J (Z " C) + " C) )■ (3 ' 2) 

The anticanonical pencil has a base point at the point ^bf x determined by the equations 1x921 = 
256 (5z) _1 , U922 = which is the lift to Sg of the point 6| n of the anticanonical pencil in Sg- The 
blowing up of 5s in the point bs, which is not the base point of the anticanonical pencil, causes 
Ewq2 to be infinite along Lg \ the line determined by the equation 11921 = 0. In turn this 
forces the energy function E{z) to have a pole at the point z = C where the solution u(z) of the 
Boutroux-Painleve equation has a pole. 



The equations for wgi, Ewgi, and Ewgi (given in Section A. 10) imply in combination with 
the Taylor expansion for ugn(z) and 1*912(2) in (3.2) that E = — 2 2 (5£) -1 ( z ~ C) _1 + 2 _7 a — 



22 (5 C)" 2 + 0(z - C) and E = -2 2 (5C)" 1 {z - C)" 2 (1 + 0({z - () 2 )). Combination of these 
asymptotic expansions for E and E leads to 

E{z) = -2 2 (5 C)- 1 (z -()-' + 2- 7 a -22 (5 C)~ 2 + 0((z - ()/(), (3.3) 

where the remainder term is uniform for bounded an d a. It follows that the energy E(z), 
although it has a pole of order one at z = Q, is close to 2~ 7 a if \z — Q\ is large compared to 1/|C|, 
and E(z) = 2~ 7 a + 0(C _1 ) if z — £ and [z — C) _1 are bounded. That is, for large \z\, E(z) is well 
approximated by 2~ 7 a as soon as z leaves the disc centered at z = £ with radius of small order 
1/ICI) where the approximation improves when \z — £| increases to order one. 

4. The solutions near the set where the vector field is infinite 

In this section, we consider the vector field near the infinity set I and show that it is repelling for 
the flow. We also construct the limit set for each solution and show that it is a non-empty, compact 
and connected subset of Sg that remains invariant under the autonomous flow. As a corollary, we 



prove that every solution of Equation (1.1) must have an infinite number of poles in the complex 
plane. 

Let S denote the fiber bundle of the surfaces Sg = Sg(z), z G C\{0}, in which the time-dependent 
Painleve vector field v z , in the Boutroux scaling, defines a regular (= holomorphic) one-dimensional 

vector subbundle V of S. For each z € C \ {0}, let I(z) := Ui=o % \ z ) ^ e the infinity set, the 
set of all points in Sg{z) where v z is infinite, that is, where V is "vertical" (or tangent to the fiber). 
If X denotes the union in S of all I(z), z £ C \ {0}, then S \ I is Okamoto's "space of initial 
conditions", fibered by the surfaces Sg(z) \ I(z), the open subset of S of all points in S where V 
is transversal to the fibers, and therefore defines a regular infinitesimal connection in the bundle of 
the Sg(z) \ I(z), z € C \ {0}. Instead of using the coordinate-invariant description of a bundle of 



surfaces with a connection, we will analyse the asymptotic behavior, for \z\ —> oo, of the solutions 
of the Painleve equation in the Boutroux scaling, by studying the z-dependent vector field in the 
coordinate systems introduced in Section [2] The solution curve in S will we denoted by 7 = 7(2), 
whereas the corresponding solution of the Boutroux-Painleve differential equation is denoted by 
u(z). Note that u(z) is equal to the first coordinate u\(z) of 7(0) in the (u±, U2) coordinate system, 
the 01-coordinate system. 

In this section we begin with an asymptotic description of the solutions near the locus X where 
the vector field is infinite. In the notation we often drop the dependence on z of the surfaces 
Sg(z). All order estimates will be uniform in z for z bounded away from zero. Near the part 
I \ Lg 1 ^ = U^=o U (Lj \ Lg 1 ^) of I we will use the function 1/E, where E is the energy, as 

an indicator for the distance to /, whereas near the remaining part Lg of / we switch to wg2 in 
the 92-coordinate system. See the first statement in Lemma 4.2 The function 1/E is no longer 
useful as an indicator function near because contains the lift (hf 11 )^ to Sg of the base 
point 6| of the anticanonical pencil, and E takes all finite values near 

(6 8 ell ) (1) - One of the points 
of the proof is that wg2 is approximately constant when the solution runs closely along 



Lemma 4.1. Let 

1=0 

For every e > there exists a neighborhood U of I 6 in Sg such that \ {E/E)/{— 6/5 z) — 1| < e in U 
and for all z 6 C\ {0}. For every compact subset K of Ll \ there exists a neighborhood V of 
K in Sg and a constant C > such that \ (E/E) z\ < C inV and for all z 6 C \ {0}. 

Proof. Because J 6 is compact, it suffices to prove that every point of it has a neighborhood in Sg 
in which the estimate holds. The quantity r := (5 z {E/E) + 6)/8 = —2u\/E is equal to 

rm = «021 2 /(4 + 2U021 2 - UQ21 U022 2 ), 

ros = u 3i 2 u 032 /(-u 3i + 2n 3i 2 -u 32 + 4-U032 3 ), 

2 2 11 1 a 2,r, 2 2\ 

ni = win U112 /(-n m + 4mi2 +2u m u\vi ), 
ri2 = ^i2i 2, "i22/(-l + 2ni2i 2 ni22 + 4ni2i 2 tii22 3 ), 

^21 = «211 2 1i212 3 /(-^211 +4n 2 12 + 2n 2 n 2 U212 3 ), 

r 22 = n 22 i 3 ti222 2 /(-l +4n 2 2iU222 2 + 2n 2 2i 3 n 2 22 2 ), 
T3i = n 3 n 2 U3i2 4 /(4 - n 3 n + 2n 3 ii 2 n 3 i2 4 ), 

r32 = W 32 l 4 U 3 22 3 /(-l +4u 322 + 2u 3 21 4 n 322 3 ), 

^41 = «412 3 (4 + U411 M 4 i2) 2 /(-n 4 n + 2u 4 i2 3 (4 + U411 n 4 i 2 ) 2 ), 

^42 = ""421 3 (4 + U 4 2l) 2 'U422 4 /(-l + 2U 4 21 3 (4 + n 4 2l) 2 n 4 22 4 ), 

^51 = ^5i2 2 (4 + U511 u 512 2 ) 2 / (-u 511 + 2n 5 i2 2 (4 + U511U512 2 ) 2 ), 

T52 = W521 2 U522 3 (4 + U521 2 n522) 2 /(-l + 2n521 2 n522 3 (4 + li521 2 n 5 22) 2 ), 

?~6i = U612 (4 + u 6U U612 3 ) 2 / (-uen + 2n 6 i2 (4 + u 6U u 6 i2 3 ) 2 ), 

^62 = ^621^622 2 (4 + n 6 2l 3 1i622 2 ) 2 /(-l + 2n 6 i2li622 2 (4 + U 6 21 3 «622 2 ) 2 ), 

r n = (4 + n 7 n u 7 i2 4 ) 2 /(-u 7 n + 2 (4 + -urn n 7 i2 4 ) 2 ), 

r72 = ^722 (4 + u 7 2i 4 n 722 3 ) 2 /(-l + 2ii 722 ( 4 + n 72 i 4 u 7 22 3 ) 2 ) 

in the coordinate charts which cover I 6 . The part \ of I 6 is equal to the line U021 = on 
which ro2 = 0. The part \ L% of J 6 is equal to the line U121 = on which 772 = 0. The part 



4 7 } \ l£> of J 6 is equal to the line n 22 i = on which r 22 = 0. The part L£ } \ (Lf> U L 2 J ) of I 6 is 
equal to the part u^u / 4 on the line U312 = on which r%i = 0. The part \ (L^ U ^q 9 '') of I 6 
is equal to the part ^3 22 ^ 1/4 on the line U321 = on which r%2 = 0. The part \ of J 6 is 
equal to the line U421 = on which r^i = 0. The part \ of I 6 is equal to the line U521 = 
on which rs 2 = 0. The part \ of I 6 is equal to the line uq21 = on which r^2 = 0. The 

part Lg 3 ^ \ Lg 4 - 1 of J 6 is equal to the line ii7 22 = on which r-j2 = 0. This covers all of I 6 , and the 
proof of the first statement in the lemma is complete. 

For the second statement we observe that iSf 1 \ (L^ U L^) is the line ii7i 2 = 0, U711 7^ 32 on 
which r-j\ = 16/(ii7n — 32) 2 , whereas U721 = 0, U7 22 / 1/32, on which r7 2 = 16 «7 22 /(32 «7 22 — l) 2 , 

(2) (3) (2) 

is an open neighborhood of L\ ' fl L 6 in L\ ' . Note that r becomes infinite when approaching 

(2) (3) (2) . . (2) 

LjDLq on L)j , which is why Lj cannot be included in the first statement of the lemma. □ 



The function \d\ in the following lemma will be used as a measure for the distance to the infinity 
set / of the vector field. 



Lemma 4.2. Suppose z is bounded away from zero. Let q := 2E. There exists a continuous 
complex valued function d on a neighborhood of I in Sg such that d = q^ 1 in a neighborhood in Sg 
of I \ d = W92 in a neighborhood in Sg of the remaining part \ of I, and qd — > 1, 

d/vjQ2 — > 1 when approaching \ 

If the solution at the complex time z is sufficiently close to a point of \ (parametrized 
by coordinate ug 2 i ) ; then there exists a unique ( G C such that \z — £| = 0(|<i(z)| |ug 2 i(z)|) ; where 
d(z) is small and |ttg 2 i(z)| is bounded, and itg 2 i(C) = 0, that is, the solution of the Boutroux- 
-Painleve equation has a pole at z = (. In the sequel we write S := d(() = wg 2 (C) = 2 6 «g 22 (£), 
and consider 5 — > 0. We have d(z)/5 ~ 1. For large finite Rg G M>o, the connected component 
of Q in C of the set of all z G C such that |ug 2 i(-z)| < Rg is an approximate disc Dg with center 
at ( and radius ~ 2 -5 \5\ Rg, and z >->■ ug 2 i(z) is a complex analytic diffeomorphism from Dg onto 
{ueC\\u\< Rg}. 

For i decreasing from 7 to 4 we use the coordinate ^+1)21 in order to parametrize iJf ^ \ 
l \ where U(j +1 ) 21 = corresponds to the intersection point of l\ 9 ^ with l!f + ^ . The point 

on ^ with coordinate U( i+ 2)2i funs to the same intersection point when |u(j +2 ) 21 | — > 00. 

For large finite Ri G M>o, the connected component of Q in C of the set of all z G C such that 
the solution at the complex time z is close to l!f ^ \ 4-°i ^ > w ^ \ u (i+i) 2i( z )l ^ Ri> but not 
close to L^ + ^\ is the complement of Dj+i in an approximate disc Di with center at ( and radius 
~ {2^ i \5\R i ) 1 l^~ i \ where we note that |5|V(9-0/|,5|V(9-(i+i)) = |<5|-V(9-*) (8-i) » 1. More 
precisely, z *— > ti(j+i) 2 i defines a (9 — i)-fold covering from the annular domain Di \ Di + \ onto the 
complement in {u G C | \u\ < Ri} of an approximate disc with center at the origin and small radius 
- (2~ 6 \8\ R i+1 % - l ) l l^\ where u (l+1)21 (z) ~ -2^ 3 5 (z - C) 9 ^. 

For all z G D^, the largest approximate disc, we have \z — (\ « \(\ and d(z)/5 ~ 1. 



Proof. Recall that L 8 \ L 7 is determined by the equation ug 22 = and is parametrized by 
^921 G C. Moreover, Lg minus one point not on corresponds to ug 2 i = and is parametrized 

by u$22- For the study of the solutions near the part \ of /, we use the coordinates 
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(«92i, 1*922)- Asymptotically for 1*922 — > and bounded 1*921, z 1 we have 

tt 92 i ~ -2 -1 U922 -1 , (4.1) 

w 92 ~ 2 6 n 9 22, (4.2) 

w 92 /w 92 = 6(5z)- 1 + 0(u 92 2 2 ) = G(5z)- 1 + 0(w 92 2 ), (4.3) 

qw 92 ~ l-2 8 (5z)- 1 n 92 r 1 . (4.4) 



It follows from ( |4.3[ ) that, as long as the solution is close to a given large compact subset of 
4 1} \4 2) , w 92 (z) = (z/(f/ 5 w 92 (() (l + o(l)), where z/£ ~ 1 if and only if \z-(\ « \(\. In view 



of (4.2), in this situation, 1*922 is approximately equal to a small constant, when (4.1) yields that 



1x921(2) ~ 11921(C) - 2~ l 1x922"* (z—Q, and it follows that 1*921(2), the affine coordinate on L%\Lj, 
fills an approximate disc centered at U921 (C) with radius ~ R if z runs over an approximate disc 
centered at £ with radius ~ 2 | W922 1 R- Therefore, if 1 1x922(0 1 << 1/|CI> the solution at complex times 
z in a D centered at £ with radius ~ 2 | U922 1 R has the properties that along it u 922 {z) / u 922 {Q) ~ 1 
and that z 1— > 1x921(2) is a complex analytic diffeomorphism from D onto an approximate disc 
centered at 1*921 (C) with radius ~ R. If R is sufficiently large, we have £ 1x921(1?), that is, 
the solution of the Boutroux-Painleve equation has a pole at a unique point in D. After having 
established this fact, we can arrange that 1*921 (C) = 0) that is, the center £ of D is equal to the 
pole point. As long as \z — £| << |C|, we have d{z)/d{C) ~ 1, i.e., 2 6 1*922(2)/^ ~ w 92 (z)/5 ~ 1 and 
1*921(2) ~ —2~ 1 u 922 ~ 1 (z — C) ~ —2 5 5~ 1 (z — (), where for a large finite i? 8 G M>o the equation 
1^921(^)1 = Rs corresponds to \z — (\ ~ 2~ 5 \5\Rs, which is still small compared to |£| if \5\ is 
sufficently small. It follows that the connected component D% of £ of the set of all z G C such 
that 1 1*921 (z) I < i?8 is an approximate disc with center at £ and small radius ~ 2~ 5 |<5| R$. More 
precisely, z h-> 1*921(2) is a complex analytic diffeomorphism from Dg onto {ix G C | |ix| < Rs}, and 
d(z)/5 ~ 1 for all 2 G Dg. The function g(^) has a simple pole at 2 = £, but it follows from ( |4.4[ ) 
that 5(2) 1092(2;) ~ 1 as soon as 1 » \z^ ugai^)" 1 ! ~ 1^21X922(0 — C)" 1 ! = 2~ 5 |<J|/|C (2 - 01, 
that is, when \z - (\ » 2~ 5 \S\/\(\. As the approximate radius of D s is 2 -5 |<5| i? 8 >> 2~ 5 |<5|/|C| 
because R$ » 1/\C\, we have q(z)w 92 {z) ~ 1 for z G Z?8 \ -D9, where D9 is a disc centered at ( 
with small radius compared to the radius of D%. 

The set L 7 \ L 6 is visible in the coordinate system (1*821, 1*822), where it corresponds to the 
equation 1*822 = and is parametrized by 1*821 G C. The set Lg 1 '' minus one point corresponds to 
1*821 = and is parametrized by 1*822 G C. It follows from the equations which express (i*gn, 1*812) 
and (1*821, 1*822) in terms of (1*711, 1*712) that 1*822 = l/i*8ii- Also, Ugn = 1*921 — 2 8 (5z)~ 1 which 
implies that 1*921 — > 00 if and only if 1*822 — > 0. That is, the point near Lg 1 ^ approaches the 

intersection point with when ( |4.4| implies that qw 92 — > 1. Therefore the functions q^ 1 and 

1092 can be glued together by means of a continuous interpolation to a continuous function d as 
asserted in the lemma. 

Asymptotically for 1*822 — > and bounded 1*821 and 2 _1 , we have 

U 8 21 ~ -«822~ i (4.5) 
1i 82 2 ~ 2~ 1 1i 82 l~ 1 , (4.6) 
W82 ~ 2 6 1*821 U822 2 , (4.7) 

qw 82 ~ 1, (4.8) 

q/q = E/E ~ -6 (5 z)" 1 - 2 7 (5 z)' 1 i* 82 r 1 . (4.9) 
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It follows from (4.9) and (4.6) that q/q ~ —6 (5z) 1 — 2 8 (5z) 1 u$22, hence 
log{q(zi)/q(zo)) ~ log^iM,)" 6 / 5 ) 

r-zi 

-(2 8 /5) (zi' 1 u 8 22{zi) - zo' 1 u 82 2(zo) + / z~ 2 u 822 (z) dz) . 



2() 



Therefore q(zi)/q(zo) ~ 1, if for all z on the segment from zq to zi we have \z — zo\ « \zq\ and 

1^822(^)1 << \ z o\- We choose zo on the boundary of D$, when ^(zo)^ 1 5 ~ g(z o) <5 ~ i(zq ) ^92 (^0) ~ 

1, and 1^921(^0)! = -^8 implies that |u822(zo)| ~ Rs 1 << 1- Furthermore, (4.7) and (4.8) imply 

that |it82i(-zo)| ~ 2 -6 1 1^82 (-^o) I I ^822(^0) I _2 ~ 2~ 6 \5\ Rs~ 2 , which is small when \S\ is sufficiently 

small. Because Dg is an approximate disc with center at £ and small radius ~ 2 -5 \S\ Rg, and 

i?8 >> |CI~\ we have that |t*92i(z)| > Rs » 1 hence |us22(-z)| << 1 if z = ( + r (zo — C)> r > lj 

and |z — zo|/|zo| = ( r ~ 1 ) | 1 — C/zqI << 1 if r — 1 is small compared to the large number 1/|1 — C/^oJ- 

Th en e quations (4.8), (4.7), and q ~ yield «822 1 ~ (<5 _1 2 6 «82i) 1//2 ) which in combination 

with JZI} leads to 2 d(« 82 T T ^)/dz = -2 3 ^ 1 /2 5 hence w 8 2i(^) 1/2 ~ u 8 2i(^o) 1/2 - 2 2 (z - z ), 

and therefore ug 2 i(z) ~ 2 4 <5 _1 (z — zo) 2 if |z — zo| >> |^82i(^o)| 1 ^ 2 - For large finite i?7 G M>o the 

equation |us2i(^)| = -R7 corresponds to |z — zo| ~ (2~ 4 |<5| Rj) 1 / 2 , which is still small compared to 

|zo| ~ |C|, and therefore \z — £| < \z — zq\ + \zo — £| << |C|- This proves the statements about the 

(2) (3) 

behavior of the solution near L 7 \L & . 

The statements for 4 < i < 6 about the behavior of the solutions near the part Lf ~ i] \ L^~ i] of 
I will be proved by induction over decreasing i. The set L^f ^ \ I^-°i *' is visible in the coordinate 
system (uu + u2ii ^(£+1)22)1 where it corresponds to the equation U(j + i) 22 = and is parametrized 
by it(j_|_i)2l £ C. The set i^+i minus one point corresponds to 21 = and is parametrized by 
n (i+i)22 S C. It follows from the equations which express (i*(i+i)n, ^(i+i) 12) and 21, M (i+i)22) 

in terms of (itai, UU2) that «(i + i) 2 2 = V w (i+i)iij u 7ii = ^821 +32, u 6U = u 72 i, and u 5U = u 62 i- 



This shows that u 



(i+2)21 



00 if and only if u 



the intersection point with L\ 
Asymptotically for uu+x} 22 



(9-0 



(i+l)22 



0, that is, the point near ££-1 approaches 



and bounded ti(i+i) 2 i and z \ we have 

«(i+l)21 - " 
^(i+l)2 



2 6 u 



i)2 1 U( i+1 )22~ 

8-i „. 9-j 
(t+l)21 u (i+l)22 , 



9/? = ^/-E 



(4.10) 
(4.11) 
(4.12) 
(4.13) 



Assume that |«(j+2) 2i( z o)| = Ri+i » 1, where the induction hypothesis yields that \zq 
C| « ICI and l/(q(z )5) ~ d(zo)/<5 ~ 1. It follows from gill ), dfefc, and \u {i+l) 22 (z )| 
Vl u (i+2)2i(*b)l> that 

l«( i+ i)2i(^o)| 8 - i ~ 2- 6 l^zo))" 1 |n (m)22 r- 9 ~ 2- 6 |5| i^+x 9 -, 

which is small if \ 5\ is s ufficiently small. 



It follows from (4.13) that q{z)/5 ~ q(z)/q(zo) ~ 1 along a solution near L^ 9 ^ \ -L-i^ , as long 
as |z - z | << \zq\- Then ( |4.12[ ) and ( |4.11| ) yield 



- 1 fA-!o6„ 8-i\l/(9-i) 
n (i+l)22 ~ (t> 2 «( i+ l)21 J , 



which in combination with (4.10) leads to 



(9 - •) d(n (4+1)21 1 /( 9 -))/dz = -(9 - i) 2 - 1+6 /( 9 -) r 1 /^), 



llGIlCG 

and therefore 21 (2) ~ 2 1-3 J -1 (z — zo) 9_l if |z — 20 1 >> \uu+i)2i( z q)\ • For large finite 

i?i E M>o the equation |«(i+i) 2 i(.z)| = ^ corresponds to |z — zo| ~ (2 3_ * |£| Ri) l /^~ % \ which is still 
small compared to \zg\ ~ |£|, and therefore |z — C| < |z — zq\ + \zg — (\ « \(\. □ 

The following corollary implies that the infinity set I of the vector field is repelling. This in turn 
implies that every solution which starts in Okamoto's space Sg\ I remains there for all complex 
nonzero times. 

Corollary 4.3. For every e\ > 0, < e 2 < 6/5, and < £3 < 1, there exists a 5 £ M>o such that 
for every solution we have that if \zq\ > ei and \d(zo)\ < S, we have the following conclusions. Let 
p denote the supremum of all r > \zq\ such that \d(z)\ < 5 whenever \zo\ < \z\ < r. Then 

i) p is bounded above by the inequality 5 > |d(zo)| (p/\ z o\) 6 ^ 5 ~ e2 (1 — £3)- 

ii) If \zo\ < \z\ < p, then d(z) = d(z ) (z / z )^ 5+£ '^ (1 + s 3 (z)), where \e 2 (z)\ < e 2 and 

\£3(z)\ < £3- 

iii) If \z\ > p then \d(z)\ >8(l-e 3 ). 



Proof. It follows from Lemma |4.2| that for every solution close to / the set of all not too small 
complex times z such that the solution is not near I 6 is a union of approximate discs of radius of 
order \d\ 1 ^ 2 where the distance between the discs is at least of order Id) 1 / 3 , where |d| 1/,3 /|(f| 1//2 >> 1. 
Therefore, if the solution is near I e at the complex times zq and z\, and is near / for all complex 
times z such that \z\ is between \zg\ and \z\\, there is a path 7 from zq to z\ such that for all z on 7 
we have that the solution at time z is near J 6 , and 7 is C 1 close to the path [0, 1] 3 1 1— > e lo+t 



where li = logZj. Then Lemma 4.1 implies that q{z\) = qjz p) (zx/zq) 6 / 5 +°( 1 ) (1 4- o(l)), hence 



d(z) = d(zo) (z/ zo) 6 ^ 5+ °^ (1 + o(l)). Because Lemma 4.2 implies that the ratio between the 



values of d remains close to 1 if the solution stays close to U Lg , it follows that d{z\) = 
d(zo) (zi/zo) 6//5+ °^ 1 ^ (1 + o(l)) if the solution is close to / at all complex times z such that \z\ is 
between \zg\ and \z\\. The corollary follows from these estimates. □ 

Remark 4.4. The substitutions 

Ul (z) = ((5/4)z)- 2 / 5 yi(((5/4)z) 4 / 5 ), 

u 2 (z) = ((5/4)z)- 3 / 5 y 2 (((5/4)z) 4 / 5 ) 

in the beginning of Section [2] lead in combination with (2.3) to 

E(z) = 2- 1 ((5/4)z)- 6 / 5 y 2 (((5/4)z) 4 / 5 ) 2 

-2 ((5/4) z)-^ yi(((5/4) z)^f - ((5/4) z)^ yi(((5/4) z)^). 

Because the solution (yi(x), y 2 (x)) of the Painleve system is single valued, we have the analytic 
continuation formula 

E(z e 5 ^/ 2 ) = -E(z), (4.14) 



analogous to (2.6). Because also (z e 57ri//2 ) 6//5 = — z 6 / 5 , the asymptotic formula E(z) / E(zq 

-(: 

J 7 



(z/zq) 6 / 5 along solutions close to the part / \ (L^ U Ly ) of the inifinity set / is consistent with 



fl4T4) . 



Because the substitutions of coordinates in Section [2] depend in a polynomial way on z -1 , the 
bundle of the complex projective algebraic surfaces Sg(z), z £ C\{0} extends to a complex analytic 
family Sg = Sg(z), z 6 P 1 \{0}, where the complex projective line P 1 is identified with the Riemann 
sphere <CU{oo}. The surface 59(00) over the point 00 6 P 1 is obtained by blowing up P 2 nine times 
as in the definition of Sg(z), where in the formulas for the base point b$(z) and the coordinate 



systems (ugn, M912) and (1*921, ^922) the coefficient 1/z is replaced by zero. Because 63(00) = 6g U , 
the base point of the anticanonical pencil defined by w and Ew, the limit surface 59(00) is equal 
to the rational elliptic surface obtained by blowing up the base points of the anticanonical pencil. 
The Boutroux-Painleve vector field converges for z — > 00 to the vector field of the autonomous 
Hamiltonian system u\ = 112, U2 = 6 u\ 2 + 1 with Hamiltonian function equal to the energy function 



E in (2.3). The function u\(z) satisfies the Weierstrass equation (ui) — 4u\ — 2u\ = 2E, which is 
why in the sequel we will use the function q := 2 E instead of the energy function E. The function 
q defines the elliptic fibration q : Sg(oo) — > P , where the fiber I(oo) = (^({oo}) = lim^oo I(z) 
over q = 00 is a singular fiber of Kodaira type II* . The vector field of the autonomous Hamiltonian 
system is regular in the limit fiber 69(00) \ 1 (00) of Okamoto's space of initial conditions, and 
infinite on I (00). The function q is constant on its solution curves, and each non-singular fiber is 
an elliptic curve where the time parameter of the solution leads to an identification of the fiber 
with C/P(q), where P(q) denotes the period lattice of the flow at the level q. The —1 curve £9(00) 
which appears at the last, the ninth blowup is a global holomorphic section for the elliptic fibration. 
Starting at the complex time z = 0on the unique intersection point of the level curve with Lg(oo), 
the period lattice P(q) is equal to the set of all z G C such that the solution of the autonomous 
Hamiltonian system hits Lg(oo). In view of the Weierstrass equation (u\) 2 — 4«i 3 — 2u\ = q, and 
the fact that U\(z) has a pole at z = ( if and only if £ G P(q), the wi-coordinate of this solution is 
equal to the Weierstrass p function of the lattice P(q), and we recover the fact that hitting Lg(oo) 
corresponds to u\ having a pole. The equilibrium points of the autonomous Hamiltonian system are 
the points in the affine (ui, U2)-charts determined by the equations 112 = and 6u\ 2 + 1 = 0. The 
corresponding singular values are q = — 4u± (—1/6) — 2u\ = — (4/3) ui = ±i ^8/27. For each 
of these two finite singular values of q we have a singular fiber of Kodaira type Ii. Therefore the 
configuration of the singular fibers of the rational elliptic surface 59(00) is II + Ii +Ii, the second 
item in Persson's list [23] pp. 7-14] of configurations of singular fibers of rational elliptic surfaces. 
It also occurs on p. 121 in the classification of Schmickler-Hirzebruch [23] of all elliptic fibrations 
over P 1 with at most three singular fibers. 

The following definition is a complex version of the concept of limit sets in dynamical systems. 

Definition 4.5. For every solution C \ {0} B z 1— > U(z) G Sg(z) \ I(z), let Qjj denote the set of all 
s G 5*9(00) \ I(oo) such that there exists a sequence Zj G C with the property that Zj —> 00 and 
U(zj) — > s as j — > 00. The subset £lu of S^oo) \ I(oo) is called the limit set of the solution U. 



Corollary 4.6 below is analogous to Coddington and Levinson [H Th. 1.1 and 1.2 in Ch.16]. 



Corollary 4.6. There exists a compact subset K of Sg(oo) \ 1 \oo) such that for every solution 
U the limit set Q,u is contained in K. The limit set is a non-empty, compact and connected 
subset of K , invariant under the flow of the autonomous Hamiltonian system on Sg(oo) \ I(oo). 
For every neighborhood A of in Sg there exists an r > such that U(z) G A for every z G C 
such that \z\ > r. If Zj is any sequence in C \ {0} such that Zj — > 00 as j —> 00, then there is a 
subsequence j = j(k) — > 00 as k — >• 00 and an s G such that U(zj^) — )■ s as k —> 00. Finally, 
for every solution U the limit set Vtu is invariant under the transformation T of Sg (00) which in 
the coordinate system (u\, U2) is given by (u\, U2) > (—u\, 1U2), when q 1— > —q and E 1— > —E. 

Proof. For any 5, r G M>o, let K$ r denote the set of all s G Sg(z) such that \z\ > r and |<i(s)| > 5. 
Because Sg is a complex analytic family over P 1 \{0} of compact surfaces Sg(z), z G P 1 \{0}, K$ iT . is 
a compact subset Sg. Furthermore Ks )T is disjoint from union of the infinity sets I(z), z G P 1 \ {0}, 
and therefore K& r is a compact subset of Okamoto's space Sg \ Sg i00 , where the latter is viewed as 
a complex analytic family of non-compact surfaces over P 1 \ {0}. When r f 00, the sets Ks, r shrink 
to the set Ks t oo of all s G £9(00) such that \d(s)\ > 5, which is a compact subset of 59(00) \ I(oo). 
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It follows from Corollary 4.3 that there exists 5 G M>o such that for every solution U there exists 
ro G M>o with the property that U(z) G Ks,r for every z G C such that |z| > ro. In the sequel, 
let r > ro, when it follows from the definition of K 5 r that U(z) G K$^ r whenever \z\ > r. Let 
Z r := {z G C I \z\ > r} and let flu,r denote the closure of U(Z r ) in Sg. Because Z r is connected 
and U is continuous, U(Z r ) is connected, hence its closure flu,r is connected. Because U(Z r ) is 
contained in the compact subset Ks )T , its closure flu,r is contained in Kg ir , and therefore fi[/,r is 
a non-empty compact and connected subset of Sg \ 5g i00 . Because the intersection of a decreasing 
sequence of non-empty compact and connected sets is non-empty, compact, and connected, and 
the sets flu,r decrease to flu as r f 00, it follows that flu is a non-empty, compact and connected 
subset of 59. Because flu,r C Ks )T for ah r > ?"o> anci the sets L^ shrink to the compact subset 
iT^oo of 59(00) \ /(oo) as r t 00, it follows that flu C Kg >00 . This proves the first statement in 
the corollary with K = Kg oo- Because flu is the intersection of the decreasing family of compact 
sets flu,r, there exists for every neighborhood A of flu in Sg an r > such that flu r C A, hence 
U(z) G A for every z G C such that |z| > r. If Zj is any sequence in C\{0} such that \zj\ — > 00, then 
the compactness of Kg r , in combination with U(Z r ) C Kg r , implies that there is a subsequence 
j = j(k) — > 00 as k — > 00 and an s G i^<5 jr such that U(zj^) — > s as k — > 00, when it follows from 
the definition of Qu that s G 

We next prove that is invariant under the flow 3>* of the autonomous system Hamiltonian 
system. Let s G Of/, when there is a sequence G C \ {0} such that Zj — > 00 and U(zj) — > s as 
j — > 00. Because the z-dependent vector field of the Boutroux-Painleve system converges in C 1 to 
the vector field of the autonomous Hamiltonian system as z — > 00, it follows from the continuous 
dependence on initial data and parameters for first order ordinary differential equations, see for 
instance Coddington and Levinson [H Th. 7.4 in Ch. 1], that the distance between U(zj + 1) and 
<& l {U(zj)) converges to zero as j — > 00. Because $*(C/(zj)) — > <&*(s) and Zj — > 00 as j — > 00, it 
follows that U(zj + t) — > $*(s) and Zj + 1 — > 00 as j — > 00, hence 3>*(s) G Qjj. 



The invariance of under the transformation T follows from ( |2.6[ ) and (2.3). □ 

Corollary 4.7. Every solution of the first Painleve equation has infinitely many poles. 

Proof. Let u(z) be a solution of the Boutroux-Painleve equation with only finitely many poles, 
U{z) the corresponding solution of the system in Sg \ Sg t00 , and Qjj the limit set of U. According 
to Corollary |4.6[ flu is a compact subset of Sg \ Sg j00 . If Qu intersects the pole line Lg in a point 
p, then there exist z with \z\ arbitrarily large such that U(z) is arbitrarily close to p, when the 
transversality of the vector field to the pole line implies that U(() G Lg for a unique £ near z, 
which means that u{z) has a pole at z = Q. As this would imply that u(z) has infinitely many 
poles, it follows that fljj is a compact subset of Sg \ (6*9,00 U Lg). However, 5g i00 U Lg is equal to 
the set of all points in Sg which lie over the line Lq at infinity in the complex projective plane, and 
therefore Sg \ (Sg >oa U Lg) is the affine (u\, U2) coordinate chart, of which is a compact subset, 
which implies that u\{z) = u(z) and U2(z) remain bounded for large \z\. In view of the theorem on 
removable singularities it follows that u\{z) and ui{z) are equal to holomorphic functions of 1/z in a 
neighborhood of 1/z = 0, which in turn implies that there are complex numbers U\ (00), 1^2(00) such 
that u\{z) — > ui(oo) and ^(z) — > 7x2(00) as |z| — > 00. In other words, Qu = {("1(00)3 ^2(00))}- 
Because the limit set Qjj is invariant under the autonomous Hamiltonian system and contains only 
one point, this point is an equilibrium point of the autonomous Hamiltonian system. That is, 
7x2(00) = and ui (00) is equal to one of the two zeros c of 71 1— > 6 v? + 1 . According to the last 
statement in Corollary |4.6[ (— c, 0) G flu if (c, 0) G flu, where (— c, 0) 7^ (c, 0) because c 7^ 0. This 
contradiction with flu = {(c, 0)} completes the proof. □ 



In Lemma 5.18 more information will be given about the asymptotic distribution for large |£| of 
the poles £ of the solutions y(x) of the first Painleve equation. In the remainder of this section we 
discuss, for the solutions U(z) close to the infinity set /, the asymptotic behavior of the set of z 



such that U (z) is close to L^f *' for < i < 3, extending the description for 4 < i < 8 in Lemma 



As in Lemma 4.2 one finds concentric approximate discs D±, D2, and approximate discs Dq of 



small radii such that the connected component of the set of all z G C such that the solution in Sg is 
close to Lj 8 ^ \ Lip, iSp \ (L® U £3 )> and \ is equal to D\, D2 \ D\, and Do, respectively. 

More precisely, L\ \L 2 is visible in (1/121, 7*122) chart, where it is defined by 7x121 = and 
parametrized by 7*122- We have 7x122 ~ ^m -1 = w\2^ 1 ^ 2 ~ Q 1 ^ 2 = d -1 / 2 , where q/q = E/E ~ 
—6 (5 z) _1 , hence d is approximately constant. Therefore each connected component of the set of all 
z G C such that 1 1^122 (^) | < R\ is an approximate disc D\ of radius ~ \d\ 1 ^ 2 R\, and z h-» 7*122(2) is a 
complex analytic diffeomorphism from D\ onto {u G C | \u\ < R\}. Furthermore, L^\L^ is visible 
in (7*221, U222) chart, where it is defined by 7x221 = and parametrized by U222, whereas the part of 
in this chart is defined by 7*222 = and parametrized by 7*221 = txiii = ^031^032 1 = ^122 _1 5 
hence | U122 1 > Ri corresponds to |u22i| < Ri 1 - We have 7x222 ~ 2 7*22i _1 = 2 ti?22 -1 / 4 x*222 1//2 , where 
W22 ~ Q 1 = d and q/q = E/E ~ — 6(5z) _1 yields that (tx222 1//2 )* ~ d~ l / A with approximately 
constant d, hence 1*222(2) ~ (u222(zo) 1 ^ 2 + d -1 / 4 (z — zq)) 2 . Let R2 be a large finite positive real 
number. As |u22i| < Ri^ 1 corresponds to |tx 2 22| ~ |d| 1//2 l^il" 2 > \d\ 1/2 R\ 2 , and |d| 1/4 /\d\ l l 2 » 
1, the mapping z *— > U222( z ) is a twofold covering from the complement of D\ in an approximate 

disc D 2 of radius ~ \d\^ 4 R 2 1/2 onto {n G C \d\ 1 ' 2 R x 2 < \u\ < # 2 }- Finally L ( Q 9) \ L,f ] is visible 
in the (7x021, 7x022) chart, where it is defined by 7x021 = and parametrized by 7*022- We have 
7x022 ~ 6txo2i _1 = — 6tx>02 -1 / 3 , where wq2 ~ 4g _1 = 4d, and q/q = E/E ~ —6 (5z) _1 yields that d 
is approximately constant, with 7x022 ~ — 2 1//3 3 <i _1 / 3 . Let Rq be a large finite positive real constant. 
Then each connected component of the set of all z G C such that 1 7*122 (z) | < Ro is an approximate 
disc Dq of radius ~ 2~ 4 / 3 3 _1 |d| 1//3 Rq, and z 1— > 7x122(2) is a complex analytic diffeomorphism from 
Dq onto {?x G C I |«| < Ro}. 

In order to understand the location of the concentric discs -Dg C Dj C Dq C D5 C D4, the 
concentric discs D\ C D2 and the discs .Do in the complex z-plane, we first describe the situation 
for the solutions of the autonomous Hamiltonian system, which have a similar behavior near the 
infinity set /, but in addition has the function q = 2 E as a constant of motion. Recall that for 
each q G C \ {±i -y/8/27} the level set in Sg(oo) of the function q is an elliptic curve C q such 
that z 1 — y & z (a(q)) defines an isomorphism from C/P(q) onto C q . Here <3? 2 denotes the flow of 
the autonomous Hamiltonian system defined by the function E = q/2, a(q) is the unique point in 
C q n Lg(oo), and P(q) is the period lattice of the flow at the level q. In the (?xi, 7x2) coordinate 
system, the first coordinate tx(z) = ui(Q z (a(q))) of <fr z (a(q)) is the solution of the autonomous 
differential equation d 2 7x/ dz 2 = 6tx 2 + 1 such that tx 2 — 4tx 3 — 2 7* = q and u(z) has a pole at 
z G P(q), and no other poles. Therefore u(z) = p p .Jz), the Weierstrass p function defined by the 
lattice P(q). 

Locally the period lattice has a Z-basis pi(q), P2(q) depending in a complex analytic fashion 
on q. The period functions p(q) = Pi(q) satisfy the homogeneous linear second order differential 
equation 

^1+ Mq + ^ p = (415) 

dq 2 + 8 + 27 q 2 dq 4(8 + 27 q 2 ) P ' 1 ' 

as an application of Bruns [HI p. 237, 238] to g2(q) = —2 and 53(9) = —q. 
Lemma 4.8. For large \q\ the period lattice P(q) has a Z-basis of the form 

Pl(q) = q- 1 ' G a{l/q) + q^/H{l/q), 

p 2 (q) = g-V6 e 2Ti/6 a (l/ g ) +g -5/6 e^i/6 
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where a and 6 are complex analytic functions on an open neighborhood of the origin in the complex 
plane such that a(0) = - i 2" 1 vr" 1 T(l/3) 3 and 6(0) = i 2 4 3" 3 / 2 vr 2 r(l/3)~ 3 . Here T denotes 



Euler's Gamma function. The differential equation (4-15) in combination with the explicit values 
of a(0) and 6(0) leads to a successive determination of the coefficients in the Taylor expansions at 
the origin of the functions a and b. 

For < i < 8 and i / 3 the point Q z (o~(q)) is near L^f *' if and only if z belongs to the 
aforementioned approximate discs Di with 5 = q^ 1 . The centers of the concentric discs D% C D- ( C 
Dq C -D5 C -D4 are at the points of the period lattice P(q) that are the pole points of the solution 
u(z) = p P , (z) of the autonomous differential equation u = Qu 2 + 1. The centers of the concentric 
discs D\ C Z?2 are at the zeros of u(z) = p p .Az). The zeros of u(z) are close of order smaller 

than |g| -1 / 6 to the points (pi(q) +P2(<?))/3 and 2 (pi(q) + P2{q))/3 modulo P(q). The centers of the 
discs Dq are at the zeros of the derivative u(z) = dp p ,Jz)/ dz, and close of order smaller than 

q~ l l & to the points p\{q)/2, p2{q)/2, and (pi(q) + f?2(<?))/2 modulo P(q). 

Proof. Any period along a closed path j q on the curve u 2 = 4 u 3 + 2 u + q is equal to 

p = j{Au 3 + 2u + q)~ l/2 du = 4~ 1/3 q 1/3 q~ 1/2 J (v 3 + 2 • 4~ 1/3 g~ 2/3 v + 1)~ 1/2 dv. 

Asymptotically for q — > oo we have 

(v 3 + 2 ■ 4- 1 / 3 q- 2 ' 3 v + I)- 1 ' 2 = (v 3 + I)- 1 ' 2 (1 + 2 • 4- 1 / 3 v /( v * + 1))^ 

= ( w 3 + 1} -l/2 _ 4 -l/3 g -2/3 {v 3 + jj-3/2 „ + O^" 4 / 3 ). 

Furthermore d(v 3 + l) -1 / 2 / du = — \ (v 3 + l) -3 / 2 3 v 2 , and therefore an integration by parts yields 

I (v 3 + I)" 3 / 2 v dv = -(2/3) j> {{v 3 + I)- 1 / 2 )' v- 1 dv = -(2/3) j> (v 3 + l)" 1 / 2 v~ 2 dv. 

Therefore p(q) = q^ 1 ^ a + q~ 5 / 6 (3 + 0(g~ 3 / 2 ), where 

a = 2- 2 ' 3 j> (v 3 + 1)~ 1/2 dv and /3 = 2" 1 / 3 3" 1 ^ (u 3 + 1)~ 1/2 v~ 2 dv, 

where the integration is over a closed path on the elliptic curve v 2 = v 3 + 1 homotopic to jg. 
Here the elliptic curve is the compact one obtained by adding one point at infinity to the curve 
v 2 = v 3 + 1 in the affine (v, v) plane. In the next computations we use the well-known formulas 

/•oo 

r(p) := / e~* tP- 1 dt, (4.17) 
Jo 

r(p + i) = P r(p), (4.18) 

n 1 -) = vr 1 / 2 , (4.19) 



B(pi,p2) := f pi-Hl-ty*- 1 dt 
Jo 



00 



sP^- 1 (s + ds = r M r Xggl ) (4.20) 

o r(pi+p 2 ) 

r(2p) = 2 2p - 1 vr- 1 / 2 r(p)r(p+ i), and (4.21) 

r(p)r(i- P ) = — ^ (4.22) 
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for Euler's Gamma function, where (4.20), (4.21), and (4.22) are Euler's Beta function, Legendre's 
duplication formula, and the reflection formula for the Gamma function, respectively. The second 
identity in (4.20) follows from the substitution of variables t = s/(s + 1) 

As our first loop we take the closed path which doubly covers the real ^-i nterv al from — oo to 
—1. The substitution of variables v = — (s + l) 1 / 3 then yields in view of (4.20) and the other 
identities for the Gamma function that a = ±i 2 1 / 3 3~ 1 B(l/2, 1/6) = ±i 2" 1 vr" 1 T(l/3) 3 and 
0=Ti 2 2/3 3- 2 B(l/2, 5/6) = 2 4 3- 3 / 2 tt 2 T(l/3)- 3 . 

As our second loop we take the closed path in v 2 = v 3 + 1 which doubly covers the real ^-interval 
from —1 to followed by the straight line in the v -plane from to e~ 27ri / 6 , where the substitution 
of variables w = e 27ri//3 v shows that the integral of (w 3 + l) -1 / 2 dw over the the second interval 
is equal to minus the integral of (v 3 + 1) _1//2 dv over the first interval. The first and the second 
loop in v 2 = v 3 + 1 intersect each other once, at the point (v, v) = (—1, 0), where the intersection 
is transversal, and therefore the intersection number of the first loop with the second loop is equal 
to ±1, where the sign depends on the choices of the orientations of the loops, as the elliptic curve 
is a real two-dimensional torus, its first homology group is isomorphic to Z 2 , when the fact that 
the intersection number of the two loops is ±1 implies that the homology classes of the two loops 
form a Z-basis of the first homology group of the elliptic curve. This implies in turn that the 
corresponding periods, asymptotically equal to the integrals of 4 -1 / 3 g -1 / 6 v^ 1 dv over these loops, 
form a Z-basis of the period lattice P(q)- 

The substitution of variables v = — s 1 ^ 3 yields that the integral of (v 3 + l) -1 / 2 dv over the real 
^-interval from —1 to is equal to 3 _1 J Q (1 — s) -1 / 2 s~ 2 ^ 3 ds = 3 _1 B(l/2, 1/3), and therefore the 
integral over the second loop leads to a = (1 - e 2n i/3 ) 2 1 / 3 3 _1 B(l/2, 1/3). As 

g(l/2, 1/3) = r(l/2)r(l/3)r(2/3) sin(7r/6) 1/2 
B(l/2, 1/6) r(l/2)r(l/6)r(5/6) sin(7r/3) 

and (1 — e 2?r '/ 3 ) 3 -1 / 2 = — i e 2,ri//6 , we arrive at the conclusion that, if in the above ± we choose 
the minus sign, the second period is asymptotically equal to e 27ri / 6 times the first period. 

Because the differential equation (4.15) for the periods has q = oo as a regular singular point, 
and the two solutions A = —1/6 and A = —5/6 of its indicial equation A (A — l) + 2A + 5/36 = do 
not differ by an integer, each solution of (4.15) is of the form g _1//6 a(l/q) + g~ 5 / 6 b(l/q), where a 
and b are complex analytic functions on a neighborhood of the origin, and the solution is uniquely 
determined by a(0) and b(0). See for instance Coddington and Levinson [H Ch. 4, Sec. 4]. It 
follows that our first period fits the description for pi(q). The analytic continuation of pi(q), when 
the small 1/q runs around the origin once, is equal to P2(q)- Because P2(q) asymptotically agrees 
with our second period, it is equal to it. Therefore the periods pi(q) and P2(q) described in the 
lemma form a Z-basis of the period lattice P(q). This completes the proof of the first paragraph 
in the lemma. 

For the second paragraph we observe that we took as the center of the discs Dg the pole points. 

(8) (7) 

The line L\ \ L 2 is visible in the (um, 1*122) chart, where it is defined by U121 = and is 
parametrized by ni22- As U121 = U2~ l and U122 = u\, it follows that the centers of the discs Z?i, 
which correspond to M122 = 0, correspond to the zeros of u(z) = p p ,Jz). The line \ 
is visible in the (^021, ^022) chart, where it is defined by U021 = and parametrized by «022- As 
^021 = ui^ 1 and U022 = ""I -1 U2-, it follows that the centers of the discs Dq correspond to the zeros 
of the derivative U2(z) = ii(z) = dp p{q) {z)/ dz of the solution u(z) = p p(?) (z) of the autonomous 
differential equation u = 6 u 2 + 1 . 

For large \q\ and z not in one of the aforementioned small discs Di, the solution of the autonomous 

Hamiltonian system in Sg \ I is close to the part of /. In the (1*311, 1*312) chart, is 
defined by U312 = and parametrized by U311. The base point (1x311, H312) = (4, 0) corresponds to 



L 3 , the origin (1*311, 11312) = (0, 0) to L 3 nL , whereas (14311, 14312) = (00, 0) corresponds 

to (14321, 11322) = (0, 0) hence to L3 PlLg . We have 14311(2:) = u(z)~ 3 u(z) 2 and u 2 —Au 3 —2u = q = 
constant, which suggests the rescaling z = z$ + q~ l l & t and u(z) = q 1 ^ 3 v(q 1 ^ e (z — zq)). Then 
11311 = i>~ 3 (v') 2 , v" = 6v 2 +q~ 2/s , and {v') 2 = 4 v 3 + l+2 q~ 2/3 v, hence i4 3 n = 4+v~ 3 + 2 q~ 2/3 v~ 2 . 
In the limit q = 00 this leads to u 3 n(i) = v(t)~ 3 v'(t) 2 = 4+i>(t) -3 . The equation (1/) 2 = 4i; 3 +l has 
a regular hexagonal period lattice P, and we arrange that v(t) = p(t) is the solution with its poles 
at the points of P. As the poles have order two, it follows that the mapping c/p -> 4 6) -t+P 
ii3ii(£) is a sixfold branched covering, where near the point t + P = + P the mapping behaves 
as 14311 (t) — 4 ~ i 6 . At this ramification point £ + P all the six branches come together, where the 
image point 11311 = 4 corresponds to the centers of the discs D$. There are three ramification points 
t where v G C \ {0}, v' = 0, at each of which i4 311 = and n 3 ' 11 = 2v~ 3 v" = 12 v" 1 ^ 0, which 
means that at each of these ramification points two branches come together. Both ramification 
points ly over 14311 = 0, corresponding to the centers of the discs Dq. The only other ramification 
points t occur when v = 0. There are two of these and at each one three branches come together. 
Both these ramification points ly over 14311 = 00, corresponding to the centers of the discs D\. 

The group of deck transformations of the aforementioned branched covering is isomorphic to 
Z/6Z and generated by T : z + P 1— > e 27ri//6 z + P, which is an automorphism of C/P because P 
has a Z-basis consisting of p\ £ C \ {0} and The ramification points of order 6 are 

the fixed points in C/P of T, which is the single point + P. The ramification points of order 3 are 
the fixed points of T 2 which are no fixed points of T, which are the two points (pi +P2) /3 + P and 
2 (pi+j>2)/3+P. The ramification points of order 2 are the fixed points of T 3 : t+P 1— > —t+P which 
are no fixed points of T, which are the three points pi/2 + P, P2/2 + P, and (p\ + P2)j2 + P. □ 



Figure 4.1 exhibits the complex times z for which the solutions of the autonomous Hamiltonian 
system for large q are near the various irreducible components of the singular fiber I (00) over q = 00. 
The solid dots represent the concentric discs Di, 4 < i < 8, centered at the points of the period 
lattice P(q) where 111 has a pole, with respective radii ~ |g| -1 ^ 9- ^ r{ for large finite where the 
distance between the points of P(q) the period lattice is of order |f7| -1 ^ 6 3> |q| -1 ^ 5 3> • • • ^> ■ 
The solution is near when z G Dq, and near 

L (9 O^go when z € Di\D i+1 for 4 < i < 7. The 
double circles represent the concentric disks D\ and D2, centered at the zeros of u\, approximately 
equal to the points (pi{q) +P2(<?))/3 and 2 (p\(q) +P2(q))/3 modulo P(q), and of radius ~ |g| _1 ^ 2 ?"i 
and ~ Igl -1 ^ 4 ^ for large finite n, respectively. The solution is near when z G D±, and near 
L 2 \L\ when z G D2 \ D\. The single circles represent the concentric discs Do centered at the 
zeros of 142 = tti, approximately equal to the points pi(q)/2, p2(q)/2, (pi(q) +P2(q))/2 modulo P(q) 
with radius ~ |(7| _1 / 3 r for a large finite r. The solution is near when z G -Do- When z is in 
the complement of all the aforementioned discs, which happens for most of the complex times, the 
solution is near the component of I (go). 

Figure 4.2 is a contour plot of the absolute value of the Weierstrass p function defined by the 
regular hexagonal lattice generated by 1 and 5 + 5 V$ i. The areas close to the lattice points = 
the pole points of p(z) should have been blacker than black, because the second order poles are 
very big. Instead the computer program has left these areas blank. The solid level curve is the one 
of the level at the zeros of p'(z), the saddle points of |p(z)|. For our u(z) = p PM (z) the lattice 

P(z) is asymptotically equal to g _1//6 times a standard regular hexagonal lattice Po, and therefore 
14(2;) ~ q 1 / 3 p Po (z). It follows that the pits in the mountain landscape of \u(z)\ in which the zeros 
of 14(21) ly are separated from each other by mountain ridges in which the heights of the passes, 
situated at the zeros of du(z)/ dz, are of large order |g| 1//3 . 
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The following lemma implies that the solutions near the part L^' \ (L Q a; U L 2 " U L^>) of / 
of the non-autonomous Boutroux-Painleve system closely follow the solutions of the autonomous 
Hamiltonian system. 



(9) 



-(7) 



•(5) A 



Lemma 4.9. Let K be a compact subset of \ (L K J> U L y 2 n U L^') and R G 
exists a neighborhood U of K in Sg and a constant C such that the distance between x(z) and xq{z) 
is < C|zo| 1 \q\~ if \ z o\ > C, \z — zq\ < R, xq is a solution of the autonomous Hamiltonian 
system such that xq{zq) £ U and xq(z) 6 U , and x is the solution of the Boutroux-Painleve system 
such that x(zq) = xq{zq). 

Proof. Lf ] \ (4 9) u4 7) ULf } ) is visible in the («3U, U312) chart, where it is determined by U312 = 
and parametrized by u 3 i%. As U311 = 00, U311 = 0, and u 3 n = 4 correspond to L 2 9 \ 

^(5) 



(9) 



(7) 



•(5), 



Then there 



and 



L\ , respectively, we keep U311 bounded and bounded away from and 4. We multiply the vector 



field by the scalar factor q 



-1/6 



-U.311 1/2 1*312 (l*311 



2u3n 2 U3i2 4 ) 1 ^ 6 , which amounts to a 



time reparametrization along each solution curve. Lemma 4.8 implies that the dista nce between 
the periods of the solutions of the autonomous system is of order |(7| -1 ^ 6 . Lemma 



4.1 



implies that, 

along the solutions near L 3 ®\ we have q(z) = q(zo) (z/ zq)^^^ + °^ , which is asymptotically constant 



for 



Z — Zq\ « 

du 3U / dt 
du 3 i 2 / dt 
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This leads to the time rescaled system 
i 311 V2 (12 - 3 u 3 ii + 2 u 311 2 n 3 i2 4 ) (n 3 n 
-3U311 1 / 2 (^3ii-4) 5 / 6 , 



A-2u 3U 2 u 3 i2 4 



-1/6 



"311 ' «312 (-6 + «311 



-U311 2 ^312 4 ) 



(u 3U 



+(5 z) 1 U311 1 / 2 U312 2 (n 3 n - 4 - 2 ii 3 ii 2 ii 3 i2 



4 

4\-l/6 



2%i 2 ll3 i2 4 )- 1/6 



(4.23) 



Actually, we have six vector fields, one for each of the sixth roots of 1/q, corresponding to the 



fact that the curves q equal to a large constant pass six times near L 3 \ (L^ 1 U L 2 ; U Vf 1 ) . The 



(9) 



(7) 



(5h 



autonomous Hamiltonian system is obtained from (4.23 by deleting the z-dependent term from 

-1/6 



the equation for dti3i2/ dt. As U312 ~ q 
vector fields is ~ (5 z^~ l rr 1 /^ i;.qh _1 / 2 



U 3 \l 



-1/2 



^311 — 4) ' , the difference between the two 
Q "'"u 3 u [u 3 \\ — 4) 1 / 6 = 0(z~ l q^ 1 / 3 ), and therefore \x(z) — xq(z)\ = 
where t - t = 0{q 1/6 (z - z )). □ 

(5) A 



It follows from the first equation in (14231) that, for the solutions near L f \ (L^> UL^U Lf) of 



■(9), , r(7), 



both the non-autonomous and the autonomous system, the f-times needed to go from a position with 



first coordinate M311 to a position near L^nL^\ L^(lL 2 7 \ and L^'nL^' are asymptotically equal 
to the integral of -3" 1 U~V 2 (lf-4)" 5 / 6 over the [/-interval from u 3 n to 0, 00, and 4, respectively. 
Therefore the distances between the discs D{ for x and the discs Di for xq are of smaller order than 



(6) or (5) 



1 1/6 



\q\"'". As the distances between the discs for xo are of order |^| 1//6 , where near Lg we have q = d, 
see Lemma 4.2 the discs for x have asymptotically the same relative position as the discs for xq. 



Therefore Lemma 4.9 leads to a description of the solutions of the Boutroux-Painleve system near I 
es the description of the the solutions of the autonomous Hamiltonian system 
implies for the solutions of the Boutroux-Painleve system near we 



which closely resemb 
near Sg oo . Lemma 



4.1 



-1/6 



- k(2 )r 1/6 (M/N) 1/5+o(1) of the 
> \zq\ » 1, until the solution has 



have q(z) ~ q(zo) (z / zo)~ & ^ + °^ hence the order \q{ 
distances between the discs increases slowly with growing 
left the neighborhood of / in Sg where the estimates hold. 

Remark 4.10. Boutroux |2, bottom of p. 310] claimed that the quantity (Y 1 ) 2 — 4Y 3 + 12 Y, see 



Remark 2.1 remains bounded when X runs to infinity along any path with bounded argument 
not passing through any pole point of Y{X). As the function q has a simple pole at every pole 
poin t of u(z), where 1*1(2:) U2 (z) and u\(z) 2 have poles of order 5 and 4, respectively, it follows from 
( 2.5 ) that (Y'f - 4 Y 3 + 12 Y has a pole of order 5 at every pole point of Y(X), and therefore the 
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claim of Boutroux can only be valid if the final X stays sufficiently far away from the pole points 
of Y(X). Boutroux referred for the proof to [21 §10], which in turn refers back to the estimates on 
[21 p. 296, 297]. No proof is given for the existence of paths along which estimates of the form [21 
(17) on p. 296] hold, and of which the concatenation is the desired path in the X-plane running 
to infinity. In particular no analysis is given of the existence of paths along which estimates of the 
form [2 ( (17) on p. 296] hold when \D\ is large. Boutroux did not give an explicit 6/5 power law 



as in Corollary 4.3 although the last term in the right hand side of [21 (36) on p. 321] contains a 
factor X -6 / 5 which Boutroux used in order to argue that the quantity (Y' ) 2 — 4 Y 3 + 12 Y remains 
bounded. 

The asymptotic formula of Joshi and Kruskal [16, (5.18) with c = 6/5] implies that the solution £ 
of the averaged equation for their quantity E, which corresponds to our energy function E, satisfies 
£(z) ~ £(zq) (z/zq)~ 6 ^ 5+ °^ for large Joshi and Kruskal |16j did not provide estimates for E 



in terms of £. It follows from (3.3) that E has a pole of order one at the pole points z = £, and 
Lemma |4.9| together with the paragraph preceding it imply that for the solutions near the infinity 
set the pole points form an approximate regular hexagonal lattice with distance between the pole 
points of order l-E^ 1 / 6 . Therefore we cannot have E(z) ~ E(zq) (z/ zq)~ g ^ 5+ °^ for all large \E\. 



However, Lemma 4.2 and Corollary 4.3 imply this estimate if z stays away from the pole points at 



a distance of larger order than |e?|, when |d| ~ \2E\ 1 << \E\ 1 ^ 6 . 

5. Near the equilibria of the limit system 

In this section, we consider the Boutroux-Painleve system near the equilibria of its autonomous 
limit and prove several results about the solutions near these equilibria. These are the solutions 
called tronquee by Boutroux that are asymptotically free of poles near infinity in certain sectors. 
We use classical methods to determine these properties and end with a determination of their 
sequences of poles near the boundaries of pole- free sectors. 

The limit system 

111 = u 2 , 

112 = 6«i z + l, 



2 , i (5-1) 



obtained from (2.2) by substituting z = oo, has two equilibrium points, determined by u\ = e i /V6, 
U2 = 0, where e S { — 1, 1}- The linearization of the vector field at these equilibrium points is given 
by (5ui, 5u2) *->• {Su2, 12ui8u\) = (8u2, e2\/6 i Sux), which has the eigenvalues 

A± = ±(24) 1 / 4 e" i ( 1 / 2 - e / 4 \ (5.2) 

with the corresponding eigenspaces determined by 5u2 = X± 8u\. 

Because the vector field in 5g(oo) \ /(oo) has no zeros on the pole line Lg(oo), and 59(00) \ 
(I(oo) U Lg(oo)) is equal to the coordinate neighborhood where u\ and U2 are finite, these are the 
only equilibrium points o f th e autonomous Hamiltonian system in 59(00) \ /(oo). The values of 



the function q := 2E in (2.3) at the equilibrium points are equal to — e (2/3) 3//2 i, and the curves 
q = — e (2/3) 3 / 2 i are the only singular level curves of the function q in Sg(oo) \ /(oo). Both these 
singular fibers of the elliptic fibration q : 5*9(00) — > P 1 are of Kodaira type Ii. The fiber q = 00, 
equal to the infinity set /(oo) of the Hamiltonian vector field defined by the Hamiltonian function 
E = q/2, is of Kodaira type II*. The configuration of the singular fibers of the elliptic fibration 
q : 5 9 (oo) ^P 1 is II* +I1+I1. 

5.1. Perturbation of a system with a hyperbolic equilibrium point. With the substitution 
t = A+ z and an affine change of coordinates (u± , 112) >— > (p + , p~) which maps (d/\/6, 0) to (0, 0) 



and the eigenvectors of the linearization of (5.1) at (e i /\/G, 0) for the eigenvalues A+ and A_ to 



(1, 0) and (0, 1), respectively, the system (2.2) is transformed to a system 

dp/dt = u(t -1 , p) (5.3) 
20 



such that the right hand side satisfies the conditions in Lemma 5.1 below. 



this section we discuss arbitrary systems (5.3) which satisfy the conditions in Lemma 5.1 



In the remainder of 
In the 



lemmas 5.2 - 5.11, we describe the solutions p{t) of (5.3) which remain bounded for all t in an 



unbounded domain in the complex plane, where the domain is increased step by step. 

follow from 



The conclusions of the lemmas 5.2 - 5.11 



O. and R. Costin [2 Th. 1 and 2], 
n, where v is complex analytic on an 
0, and less special assumptions are made on 



which deals with systems (5.3) in arbitrary dimensions 
open neighborhood of (0, 0) in C x C n , v(0, 0) 
Lq := dv(0, p)/dp\ p =Q. The proof of [71 Th. 1] uses Borel summation. Our proofs, a sequence of 
variations on the method of Cotton, are more classical, where our step by step approach may be 
helpful in understanding all the aspects of the final picture. As we only need the case that n = 2 
and Lq has the eigenvalues ±1, and a wide generalization would require more elaborate notations 
and proofs, we did not attempt to write down the latter. In the next paragraph we summarize the 
results of the lemmas 15.21 - 15.111 



In Lemma 5.2 we prove that the solutions p(t) of (5.3) which remain small on half lines to + ^>o 

yields 



5.3 



form a one-parameter family parametrized by the complex number p(to)~ ■ Lemma 
unique solutions p&(t) and p±(t) which remain small on horizontal axes in the upper and lower 
complex half plane, respectively. In Lemma 5.5 it is established that for every e > there exists 
an r > such that the p^(t) have a common complex analytic extension to the domain defined 
by the inequalities \t\ > r and —tt/2 + e < argt < 37r/2 — e, and that this common extension 

oo. 



has an asymptotic expansion in strictly negative powers of t as |t| 



Lemma 5.6 states that 



there exists an a G C such that for each solution p{t) in Lemma 5.2 there exist r/, r > such that 
p(t) has a complex analytic extension to a small solution of (5.3) on the domain r defined by 

r(t)\ 



t\ > r, — 7r < argi < n, and |r(i)| < n, where r(t) := e t a . Furthermore, there exist unique 
C G C and dj G C 2 , j G Z>o, such that on every subdomain E where e~* is of smaller order 
than every negative power of t as |t| — > oo the function p(t) — p^{t) has the asymptotic expansion 



p(t) 

p(t) in Lemma 5.2 



Here the dj do not depend on the choice of the solution 
rgument as in the proof of Lemma 5.6 we obtain that in 
a subdomain where C\ \t\~ ei < | e - ' | < C2 \t\~ €2 for strictly positive C\, C2, e%, €2, we have an 
asymptotic expansion of the form p{t) ~ Ylh iez> T (f) h t~ l Ph, % a s |£| — > 00. Here the coefficients 
Ph,i satisfy ph^ = C h Ch,i, where the coefficients c/^j do not depend on the choice of the solution 
p(t) in Lemma 



Lemma 
I, when 
as 



5.8 



states that for each i the series F{(t) := 



h>0 



r h Ph,i 



lemma 



converges for 

5.11 implies that we have an asymptotic expansion of the form 
00 in a domain where Imt > and Fo(r(t)) remains bounded, 



5.3 

sufficiently small |' 
P(t) ~ £i> W<))* 

extending far into the domain where p(t) is bounded away from zero. This expansion will lead 
to the asymptotic determination of [7, Proposition 15] of the sequence of poles of the truncated 
solutions of (1.1) in Lemma 5.18 At some places our statements are more precise and our proofs 



more complete than those in [7j. 

Lemma 5.1. Assume that v = (v + , v~) is a 

neighborhood D of the origin in C 3 such 
dw(0, p)/dp\ p= Q = 0. Here p = (p + ,p~) G 
strictly positive real numbers So, eq, C\, C2, C3. 
and \\dw(u, p)/dp\\ < C3 ||p|| + C4 |n| if \\p\\ < 5q 



2 -valued complex analytic function on an open 
that v^{u, p) = ±p^ + w(u, p), u>(0, 0) = 0, and 
2 and \\p\\ := max{|p + |, Then there exist 

and C4 such that \\w(u, p)\\ < C\ ||p|| 2 + C2 \u\, 
and \u\ < €q. Here the last condition implies the 



preceding one for C\ = C3/2 and a suitable C2. In the sequel we will t ake D equal to the set of all 



(u, p) G C x C 2 such that \u\ < eo and \\p\\ < 5q. For solutions p of (5.3) we will always require 



that \t\ > 1/eo an d 



< 5q for all t in the domain of definition of p. 



Proof. We have < \\w(u, p) — w(0, p)\\ + ||u>(0, p)\\ and, with the notation d2w(u,p) = 

dw(u,p)/dp, \\d2w(u, p)\\ < \\d2w(u, p) — d2w(0, p)\\ + ||c?2u;(0, Because w(0, 0) = and 



d 2 w(0, 0) = 0, an application of a Taylor expansion with estimate for the remainder term to each 
of the terms between the norm signs yields the estimates in Lemma 5.1 □ 



Lemma 5.2. In the situation of Lemma 5.1, let < 5 < do, < e < eo, to, a G C, \t\ > l/e for 
every t G T := to + H^>o> | a_ | + C\ 5 2 + C 2 e < 5, and C3 6 + C^e < 1. Then there exists a unique 



a 



a 



to, a 



G C such that the solution p(t) of the differential equation (5.3) withp(to) = (a 



>s 



defined for allt £ T and satisfies \\p(t)\\ < 6 for every t £ T. For any < 6 < min{(5o, 1/Ci, I/C3}, 
the set of all (to, a~) for which the above conclusions hold is a non-empty open subset ofC 2 on which 
the function (to, a" ) 1— > a. 



t , a 



is complex analytic. Similar conclusions hold with to + M>o replaced 
by a curve r + i ct(t), where r G [Reio> oo[, a : [Reto, 00 [ — > R is continuously differentiate with 
a bounded derivative, and the curve stays sufficiently far away from the origin. 

Proof. We apply the method of Cotton [8]. 



The system (5.3) can be viewed as an inhomogeneous system of linear differential equations with 
w(t , p(t)) as the inhomogeneous term, and therefore it is equivalent to the system of integral 
equations 

rt 

e ± (*~ s ) wis' 1 , p(s)) ± da. 



p(i)± = e ± ('- T ) p(r) ± + 



(5.4) 



Let X denote the set of all continuous functions p : T — > C 2 such th at (t , p(t)) G D for every 
t G T and t h-» w^ 1 , p(t)) + is bounded on T. If the solution p[t) of (5.3) belongs to X, then we 
can let r G T run to infinity in the equation (5.4) for ± = +, and obtain 

wis' 1 , p{s))~ 



J-s 



ds, 



where f. 00 indicates the limit for T 3 r — > 00 of J*. T . Conversely, if 
equations (|5.5|) and 



P (ty 



e*°-* a 



+ 



w(s , p(s)) ds 



(5.5) 

G C and p G X, then the 
(5.6) 



for all t G T imply that pit) is a solution of (5.3) such that p(to)~ = a~ ■ 

Let F denote the integral operator which assigns to each p G X the function F(p) : T — > C 2 such 
that, for each t G T, F(p)(t) + and F(p)(t)~ are equal to the right hand side of (5.5) and (5.6), 
respectively. Let X denote the set of all continuous functions p : T — > C 2 such that ||p(t)|| < 6 
for every t G T. Then X C X, X is a complete space with respect to the metric d{p\, p^j := 
sup tgT — P2{t)\\, and the assumed estimates imply that 



\npi)(ty 



\F(p)(t) + \ 

\pmr\ 

F{P2){t) ± \ 



< 



< 



< 



e~ T (Ci 5 2 + C 2 e) dr < 5, 



+ 



e~ T (Ci <5 2 + C 2 e) dr < 5, 



(5.7) 
(5.8) 



e" T (C 3 5 + C 4 e) || Pl (t + r) - p 2 (i + r)|| dr 



< (C 3 6 + C±e)d(p 1 ,p 2 ) 



(5.9) 

for every p,p\,p 2 G X and t G T. Therefore F{X) C X and i 7 : X — )• X is a contraction, with 
contraction factor < C3 <5 + C4 e < 1, when the contraction mapping theorem implies that F has a 
unique fixed point, a unique p£l such that -F(p) = p. It follows that there is a unique solution 
p(t) = p to a -(t) of (5.3) on T such that p(to)~ = a~ and ||p(i)|| < 5 for every t £ T. As the 
arbitrary solutions p(I) of (5.3 ) ar e uniquely determined by their value p(to) at t = to, it follows 
that for every solution p(t) of (5.3) on T with p(to)~ = a~ and p(to) + 7^ a + (io 5 «~) := Pt ,a-(^o) + 
there exists t G T such that ||p(t)|| > 5. The assumptions remain verified upon small perturbations 



of to and a , and an application of the implicit function theorem yields that the solution p = p t 



depends in a complex analytic way on (to, a ). This completes the proof of Lemma 5.2 

C,5 2 + 



□ 

Lemma 5.3. In the situation of Lemma \5.l\ let < 5 < 5q, < e < eg, C\5 2 + C^t < 6, 
and C3 5 + C4 e < 1. TTien t/iere exists a unique solution p(t) = p^(t) and p = Pi(t) of (5.3) on 
E + i/e and E — i/e suc/i that < 5 for every t G E + i/e and t G E — i/e, respectively. 

Similar conclusions hold with E ± i/e replaced by a curve r + i ct(t), where r G E ; cr : E — >• M is 
continuously differentiate with a bounded derivative, and the curve stays sufficiently far away from 
the origin. 

Proof. We apply the variation on [8| p. 483] of Cotton's method. 

Let I = E + i/e, and 3^ the set of all continuous functions p : I — > C 2 such t hat (t~ , p(i)) G -D 
for every t £ I and t 1— )■ ^(t^ 1 , p(i)) is bounded on /. If the solution of (5.3) belongs to y, 
then we can let I B t — > +00 in the equation (15.41) for ± = + and let I 9 r — )• —00 in (15. 4D for 



and obtain (5.5) and 



p(ty 



w(s , p(s)) ds, 



(5.10) 



respectively. Conversely, if p G y satisfies ( |5.5| ) and (5.10), then it is a solution of (5.3). 

Let G denote the integral operator which assigns to each p G y the function G(p) : I 
such that, for each t G I, G(p)(i) + and G(p)(t) 
(5.10), respectively. Let y be the set of all continuous functions p : / 



C 2 

are equal to the right hand side of (5.5) and 

2 such that \\p(t)\\ < 5 



for every t E I. Then Y C y, Y is a complete space with respect to the metric d(p\, P2) '■= 
sup tg [/ — P2{t)\\, and estimates analogous to (5.7), (5.8), and (5.9) imply that G(Y) C Y 

and G is a contraction in Y with contraction factor < C 3 (5 + C4 e < 1 . This time the contraction 
mapping theorem yields a unique solution p : / — > C 2 of (5.3) such that < 5 for every t E I. 

This completes the proof of Lemma |5.3| □ 



Remark 5.4. According to Anosov [1, §4], "Every five years or so, someone "discovers" the theorem 
of Hadamard and Perron, proving it either by Hadamard's method or by Perron's." The theorem 
alluded to is the stable manifold theorem in dynamical systems, "Hadamard's method" refers to 
Hadamard [12], and "Perron's method" to Perron [22] . However, Perron |22[ p. 130] made perfectly 
clear that he used the method of Cotton [8]. The paper [8] seems to be little known, for no good 
reason. 



Lemma 5.5. 

7/(r) = (r 2 



For each m, n, r G M>o such that m > n, let n 
2 ) 1 / 2 when \t\ < rm(l + m 2 ) -1 / 2 , 77(7") = r (1 + m?) 



Vm, n, r 
2M/2 



R — > E be defined by 
m\r\ when rm(l + 



rrr 



; ) X ^ 2 < |t| < r (1 + m 2 ) 1 / 2 /(m — n), and ry(r) = — n |r| w/ien |r| > r (1 + m 2 ) 1 / 2 /(m — n). Let 



V := {t G C I Imt > n(Ret)}, see Figure 5.1 



If r is sufficiently large, then all solutions p^(t) in Lemma 5.3 have a common extension to a 



unique solution of (5.3), also denoted by p\{t), on the union V of V with all their domains of 



definition. There exists j 1— > Cj 
the asymptotic expansion 



Z 



>o 



' , such that the solution p(t) = Pf (t) of (5.3) on V has 



pit) 



(5.11) 



as t G V, \t\ 
that 



00. The asymptotic expansion (5.11) can be differentiated termwise in the sense 
dp(t) 



dt 



(5.12) 



=1 
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as t 6 V, \t\ —> oo. All the coefficients Cj are uniquely determined from the equations obtained by 



substituting (5.12) and (5.11) in the left and right hand side of (5.3), respectively, and using the 



power series expansion 

v(t-\ p) 

We have 



J l,3+,3-> 



J *,3+,3- 



1,3+,3- 



Cl 



dv(0,p) 



dp 



dv(u, 0) 



p=0, 



(5.13) 



(5.14) 



u=0 



Applying the complex conjugation t i— >■ t, we obtain similar conclusions for the solutions Pi(t) in 
Lemma\5. 



Proof. For each t G V we define half-lines Lf emanating from t, as follows. If Ret > 0, then 
L t + := t + R> . If < - Ret < m Imi, then L t + := {t + r - i(Ret/Imt)r | r G M> }. If Ret < 
and m Imt < — Ret, then := {t + r + i tot | r G M>o}. Finally, L(( := i?(L^^) for every 
t G V, where i? is the reflection 1 1— )• — Re t + i Im t. 

It follows that, for every t G V and s G I/f, \t\ > r and |s| > (to — ra) (1+to 2 ) -1 / 2 |t|. See Figure 

ED 

On the space of functions p : V — > C 2 such that p(t) = 0(t~ 1 / 2 ) as t G V, |t| — > oo, we 
use the norm := sup tg y \t\ 1 ^ 2 \\p(t)\\. Let r > 1/eo and < 6 < r 1 / 2 5q, and Z the space 
of continuous functions p : V — > C 2 , complex analytic on the interior of V, such that \\p\\ < 6. 



For each p G Z we define as above, where the integrals in the right hand sides of (5.5) 

and (5.10) are over Lf and , respectively. Using the Cauchy integral theorem, we obtain that 



G(p) is a complex analytic function on the interior of V, with derivative equal to dG(p)(t) =t / dt = 
±G(p)(t) ± + ur^i -1 , p(t)). Because the velocities of the parametrizations of Lf in the previous 
paragraph are at most (1 + to 2 ) 1 / 2 , we have the estimates 



\G(px)(t) 



± 



|G(p)(t) ± | 
G( P 2)(t) ± \ 



< (i + m 2 ) 1 / 2 sup (Ci <5 2 + c 2 ) i^r 1 , 

< (1 + to 2 ) 1 / 2 sup (C 3 |sr 1/2 + ^4|s| 



ILPi -P2 \s\ 



-1/2 



- 1 (1 + to 2 ) 1 / 2 Itl" 1 and Isl -1 < r" 1 for 



for every p, pi, p2 £ Z and t £ V. Because |s| 1 < (to — n 

every s G L± it follows that \\G(p)\\ < (1 + to 2 ) 3 / 4 (m - n)" 1 / 2 (Ci 5 2 + C 2 )r- 1 / 2 and \\G( Pl 
G(p 2 )\\ < (l+m 2 ) 3 / 4 (m-n)- 1 / 2 (C 3 r-- 1 / 2 +C 4 r- 1 ) ||pi-p 2 ||- Therefore, if r > l/e , < 5 < r l l 2 5$ 
(1 + to 2 ) 3 / 4 (m - n)- 1 / 2 (Ci <5 2 + C 2 ) r^ 1 / 2 < 5 and (1 + to 2 ) 3 / 4 (m - n)" 1 / 2 (C 3 r" 1 / 2 + C 4 r" 1 ) < 1, 
which for any m, n, 5 G M>o can be arranged by taking r sufficiently large, we have G(Z) C Z, 
and G : Z — > Z is a contracti on. This time the contraction theorem implies that there is a unique 
solution p = pv : V -> C 2 of (5.3) such that ||j>(t)|| < 5 |t|~ 1/2 for every t G V. 



For any a G C, successive integrations by parts yield the asymptotic expansions 

A) k 7T k t~ a 



- a ds 



oo 

E< 

k=0 



and 



~ a ds 



oo 



k=0 



for |t| — > oo, where 7Tfc := Yij=o ( a + j)- This is one of the oldest examples of asymptotic expansions 

2 + rv-' 1 



which do not converge. Substituting the inequality ||u>(,s _1 , p(s))\ 
{C\5 2 + C2) in the integrals in the right hand sides of p(t) =t 
p(t) = 0(t _1 ) as t G V, |t| — )• 00. Substituting 



< Cib(s)r + OfeM -1 < 
= G(p)(t) ± , we obtain that 



Pit) 



N-l 



t- j Cj + 0(t 



—N\ 



(5.15) 



and u = t 



u i ( P + y+ {p-y- 

G(p)(t) ± , we obtain (5.15) with 



for t € V, the power series expansion w(u, p) = E» J+ 3 -_ 
p = p{t) in the integrals in the right hand sides of p{t)^ 
N replaced by N + 1 and a unique Cjy G C 2 . This inductive procedure yields that pit) has 
an asymptotic expansion of the form (5.11) for t 6 V, \t\ 
and p 



oo. The substitution of it = £ 



in the power series expansion of v(u, p) implies that there exist bi G C such that 
p> it) = v iir\ p{t)) ~ E»=i * _< bi as t £ V, \t\ -> oo, when j£ p'(s) ds+p(t ) = p(t) ~ E£i c 3 - 
yields that hi = — (i — 1) q_i for every 2 £ Z>i. This proves (5.12) for p = py. 

We will prove next that all solutions in Lemma 5.3 and have a common extension to a 
solution of (5.3) on V. The conditions in Lemma 5.3 hold if and only if 



e < min{e , 1/C 4 , 1/4 Ci C 2 } and <5_(e) < (5 < min{£ , <5+(e)}, 



where <5±(e) := (1 ± (1 — 4Ci C2 e) 1 ^ 2 ) /2C\. If in the proof of Lemma 5.3 we replace Y by the 
space Y e> § of all continuous functions p on the upper half plane U e := {t £ C | Imt > 1/e}, complex 
analytic in the interior of U e and with ||p(t)|| < 5 for all t £ U e , then G is a contraction on Y £i( 5 
and therefore has a unique fixed point p et s £ Ye, 5- The Cauchy integral theorem implies that the 



restriction of p 6i s to the interior of U e is a complex analytic solution of (5.3). As the restriction of 
p 6j s to R + i/e is a fixed point of G : Y — >• Y, and the latter is unique, and equal to the solution 
p-j- in Lemma 5.3, it follows that the solution p^(t) in Lemma 5.3 extends to a solution p = p^ tt of 
(5.3) on U e , complex analytic on the interior of U e , and satisfying ||p(t)|| < 8 for every t £ U e . Note 



that if < e' < e, then the conditions in Lemma 5.3 hold with e replaced by e', and the restriction 
to U e i of p-f- i € is equal to t > . 

If pit) is any solution of (5.3) on some upper half plane U, and p[t) — > as t £ U, \t\ — > 00, 
then there exist 5, e as above such that U e C U and ||p(i)|| < 5 for every t £ U e , and therefore 
p\u t = p^, e - As the sector-like domain V contains an upper half plane U e , and pv(t) — > as t £ V, 
\t\ — > 00, it follows that (pv)\u e = Pt, e- This completes the proof that all solutions p+ in Lemma 



5.3 extend to a common solution pit) of (5.3) on V, of which the restriction to V is equal to py, 
and therefore p has all the properties mentioned in Lemma 5.5 □ 



Lemma 5.6. In the situation of Lemma\5.1\ let a denote the right lower corner in the 2 x 2-matrix 



Li 



d 2 viu, p) 



dudp 



+ 



<9 2 u(0, p) 



n=0,p=0 



dp 2 



ex, 



p=0 



with c\ defined by (5.14)- For rj, r £ M>o, let 



R 



rj,r ■- 



{t£C \ \t\>r and | e~* t a \ < rj), 



(5.16) 



(5.17) 



where t a = e alogi ; logi = log|i| +iargi, and —it < argi < it. For every solution pit) = Pt a -(t) 
andpit) = p&it) of (5.3) in Lemma 5.2 and Lemma 5.5, there exist r/, r £ M>o such that pit) extends 
to a solution of (5.3) on R^.r, which extension is again denoted by p to a - andp^, respectively. There 
exist j h-» dj : Z>o —> C 2 with do = (0, 1), and C = C to a - £ C, such that on any subdomain E of 
Rfj^r on which Ret/ log \t\ — > +00 as t £ E, \t\ — > 00, we have the asymptotic expansion 



(5.18) 



5=0 



2.i 



as t £ E, |i| 
that 



oo. The asymptotic expansion (5.18 ) can be differentiated termwise in the sense 



p'(t) - p\(t) ~ C e~* i Q V (-d,- + (a - j + 1) dj-i) 



(5.19) 



as i G S, \t\ —7- oo. j4ZZ the coefficients dj are uniquely determined from the equations do = (0, 1) 

and —dj + (a — j + 1) dj-i = X^i=o ^» / or a ^ — ^> where the 2 x 2-matrices Li are determined 
from the asymptotic expansion 



as t £ V , \t\ —> oo. Here 



L{t) 



Ln 



dv(t 1 , p) 



c?p 



<9u(0, p) 



p=Pt(*) 



(5.20) 

\ (i l /' ( 5 - 21 ) 
d P P =o V - 1 / 

and L\ is given by (5.16). The asymptotic expansion (5.11) for pit) = p\{t) holds for t £ E, 
\t\ —> oo, and can be differentiated termwise there. In combination with (5.18), (5.19), it follows 
that pit) 



-1 



Pt a- CO ^ as ^ e asymptotic expansions (5.11), (5.12) fort £ E ; |i| 

T/ie complex number C 



a 



to, a - ? 



determined by 



oo. 



lim 

teE, 



(p(0-Pt(0) = (o, c), 



depends in a complex analytic way on {to, a ). Conv ersely, for every C £ C there exist rj, r £ M>o 
and a unique solution p = pc : R V)T C 2 o/ (5.3) such that (5.18) holds, where pc{t) depends 



in a complex analytic way on C . 
P = Pt .a- with a~ = p(to)~ , and a" 



domain of definition onto an open subset of 1 
C = if and only if there exists a sequence tj in 

,ti j. .-a 



If we choose to such that ||p(io)|| is sufficiently small, then 
i — y C to a - is a complex analytic diffeomorphism from its open 
In particular we have p(t) = p\{t) if and only if 



such that | arg(ij-)| < ir, e J tj 



and 



oo. 



Finally, for every C £ C and e > i/iere exist n, r £ M>o swc/i i/tai 

||e't- a (pc(0-Pt(0)-(0, C)||<e (5.22) 

/or every t £ -R^r- Li particular, i/C / and we choose < e < |C|, £/ien pc(0 * s bounded 
away from zero on the set of all t £ C swc/i i/iai | e - ' i Q | = n, |t| > r, and Imt > 0, i/ie part of 
the boundary ofV n ^ r in the upper half plane and sufficiently far away from the origin. A similar 
statement holds in the lower half plane with pj- (i) replaced by p±(t). 

Proof. Let W be a suitable subdomain of V and p : W — > C 2 a solution of (5.3). Then := 
p(i) — pf(t) is a solution of the differential equation 

dy/dt = fit-\ y) :=v(t- 1 ,p t {t)+y)-v(f- 1 ,p r (t)), (5.23) 

where (i, y) i— >■ f(t , y) is a C 2 -valued complex analytic function on x B, if 5 := {y £ C 2 | 
||y|| < <5 }, where £ := 6q — sup ig jy ||pf(i)|| > 6q — br~ x l 2 . Furthermore, 



f(t-\y)~ *~ j fM 

j=0 



(5.24) 



-> oo, where each of the functio ns f j is complex analytic on B. We have /(i , 0) = 0, 
for every j. It follows from (5.20) that dfjiy)/dy\ y= o = Lj for every j £ Z>q. All 



for t £ V, \t\ 
hence fj(0) - 

aforementioned asymptotic expansions are uniform in y £ B, and can be differentiated termwise, 
arbitrarily often, with respect to t, y + , and y~ . 
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E-n 1 t~> A 



j, where the Aj are 2 x 2-matrices with Aq 



If A(t) - ^ . =0 
y = A(t) z transforms the linear differential equation dy/dt 
dz/dt = M(t)z, where L(t) A(t) 
where, for each I £ Z>q 



1, then the substitution 
L(t) y into the linear equation 
A'(t) + A{t)M{t). It follows that M(t) ~ J2T=o t~ k M k, 



j+k=i 



Lj A k 



l)A l _ 1 + Y j AjM k . 



j+k=i 



Here Aj = when j < or j > N. It follows that Mq = Lq as in (5.21 ), when the linear mapping 



A\ i — y Lq A\ — A\ Mq is surjective from the space of all 2 x 2-matrices A\ onto the space of all 
antidiagonal 2 x 2-matrices, with kernel equal to the space of all diagonal 2 x 2-matrices. It follows 
by induction on / < N — 1 that, given the Aj and Mj for j < I — 1, there is a unique antidiagonal 
matrix Ai such that the matrix Mi is diagonal. In other words, with the substitution y = A(t) z 



the differential equation (5.23) is equivalent to the differential equation 

dz/dt = g(t-\ z) := A^)- 1 f{t~\ 

where g{t~ l , z) has an asymptotic expansion for |t| — > 

N-l 

g{t-\z) = 



A{t) z) - A^t)- 1 A'{t)z, 
oo of the same nature as f(t , y), 



h(t 
dh(t 



-l 

1 



t- j Mj z + h^- 1 , 

j=0 

0(\t\ 



(5.25) 



(5.26) 



-N 



Z 



+ 0(11*1 



and 



Ofltl 



-N\ 



+ 0( 



z 



Let [i] 



dz 

oo, where Mj is a diagonal 2 x 2-matrix for each < j < N 



(5.27) 
(5.28) 
1. We have 



as t G V and \t\ 

Mq = Lq is as in (5.21), and M\ is equal to the diagonal part of the 2 x 2-matrix L\ in (5.16). 

( and fij denote the left upper and right lower corner of the diagonal matrix Mj, where 
±1 and /uj~ = a. Write 



/4 



e ± (t) 

0±it) 



t^ e ± it), 

N-l 



where 



exp \J2 1$ 



(5.29) 



Note that (5.29) implies that 8 ± (t) 



~ e 



±t 



neous linear differential equation 



Eto 1 t~ j M 3 C are given by C(i) ± = @ ± {t) 9 ± (r)" 1 C(r) 



t^i for large \t\. The solutions C(^) of the homoge- 



>j=0 



Substituting (5.26) in (5.25), we obtain an inhomogeneous linear differential equation for z with 
h(t , <z(i)) as the inhomogeneous term, when Lagrange's method of variation of constants yields 
the integral equations 



z(t)= 



8±(r) 



z(r)± + 



e ± (t) 
e±( s i 



/i(s" i ,z(s)) ± ds. 



The reasoning leading to (5.5) and (5.6) this time leads to the conclusion that z(i) is a uniformly 
small solution of (5.25) on T := to + M>o such that z(io)~ = b~, if and only if z(t) is a uniformly 



small continuous function which satisfies the integral equations 

50 Q+(t) 



zit) + = H(z)(tY 



e- 



^(s^ 1 , «(s)) + ds 



and 



zit)- = H(z)(ty 



®-(t) 

e-(to) 



6~ + 
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to 



e-(t) 



/i(s 1 , -z(s)) ds. 



(5.30) 



(5.31) 



Assume that N > 1, 7 > |6~/0~(to)|. In the next paragraphs we will prove that, if sup igT |t _1 | 
and sup 4gr I e_ * are sufficiently small, then H is a contraction in the set Z of all continuous 
functions z : T — >• C 2 such that ||jz|| := sup tgT | e* t~ a \ \\z(t\\ < 7, with respect to the metric 
(zi, z 2 ) >->• \\zi - z 2 \\. 

If in the integrand in (5.30) we substitute s = t + r, r G M>o, then (5.27) implies that | e* t~ a \ 
times the absolute value of the integrand is estimated from above by a uniform constant times 



-2t \T + T \a-f4 



l(- 



t 



t 



-JV 



+ c" 



-(*+r) 



(* + r) a | 



In view of (|5_28|), 

I e* t~ a \ \H{ Zl )(ty 



a similar estimate with the last factor ||z|| replaced by \\z\ — z 2 \\ holds for 
- H(z 2 )(t) + \. The integral of r H e~ 2r |(^±^) a ^i | over M> is uniformly 

_1 | and I e _s s a \ 
< 7 and 



bounded if the distance from T to the origin is bounded away from zero. As \s 

"«| \H(z)(t) 



are uniformly small for all s G T, there exists a < (3 < 1 such that | e* t 



I e* t~ a \ \H( Zl )(t)' 
In (I5.31D, I e* t~ 



~ — H(z2)(t) + \ < f3 \\zi — % || for all z, zi, z 2 £ Z and all t G T. 

a I times the absolute value of the integrand is estimated from above by a uniform 



constant times ( s 



-A' 



+ I e _s ||z||. A similar estimate, with the last factor ||z|| replaced by 
\\ z i ~ ^2 1|) holds for the integrand in the integral formula for | e* t~ a \ \H(z\){t)~ — H(z 2 )(t)~\. Let 
N > 1 and 7 > 1 1> / CB> C^o) I - As |s _1 | and | e~ s s a \ are uniformly small for all s G T, there exists 
a < (3 < 1 such that | e* t~ a \ \H(z)(t)-\ < 7 and | e* t~ a \ \H(zi)(t)~ - H(z 2 ){t)-\ < j3 \\z x - z 2 \\ 
for all z, z\, z 2 £ Z and all t G T. This completes the proof that H{Z) C Z and H : Z — > Z is 
a contraction with respect to the metric (z\, z 2 ) h-» sup igT | e* t _a | ||zi(t) — ^2 (*) 1 1 - Because Z is 
complete with respect to this metric, it follows that there is a unique solution z = z tQtb - : T — > C 2 
of (5.25) such that z(to)~ = b~ and ||z(t)|| < 7 | e~* t a \ for every t G T. Because the latter estimate 
implies that ||z(t)|| is uniformly small, the solution p(t) = p^(t)+A(t) z(t) of (5.3) is uniformly small 
on T, and therefore equal to the solution p tQ a -(t) in ii), with a~ = (pf(io) - + Mto) z io 5- (io))~ 
and b~ = A(to) -1 {Pt , a -(to) — Pf(*o)) _ - The assumptions remain valid upon small perturbations 
of to and b~ , and an application of the implicit function theorem yields that the solution z = 
z t ,b~ depends in a complex analytic way on (to, b~), when b~ 1— > a~ is a local complex analytic 
diffeomorphism depending in a complex analytic way on to. 

The estimates in the previous paragraph imply that the integral J t °° _ (s) _1 h(s~ 1 , z(s)) ds is 
absolutely convergent, when (5.31) implies that 0~(t) _1 z(t)~ converges to 

C :=Q-{t )b- + G^(s)- 1 /i(s- 1 , z(s)) ds (5.32) 
J t 

as t = to + T > T - > +00. Because 0~(t) ~ e~* t Q , the function e* t~ a z{t)~ has the same limit. The 
previous estimates for z{t) + = H(z)(t) + yield that e* t~ a z(t) + converges to zero, and it follows 



that e t 



h ,b- 



(t) converges to the vector (0, C) G C 2 which depends in a complex analytic way 



on (to, b ) and on (to, a ). Because A(t) converges to the identity matrix, also e t a (p t a -(t) 



^>o 



the 



Pt(t)) = A(t) (e< t" 
As b~ = z(t )~ 
equation (5.31) implies the integral equation 

z(t)~ =H(z)(ty := 0"(t) 

If T~L(z)(t) + := H(z)(t) + , then then for every C G C the integral operator "H = He is a con- 
traction in if 7 > |C| and the numbers sup tgT |t _1 | and sup tgT | e~* t a | are sufficiently small. 
The unique fixed point z = zc of He is equal to the unique solution z(t) of (5.25) such that 



" %(),&-(*)) converges to (0, C tQiCL ~) as t = t + r, t ->■ +00. 
and the previous remains true if we replace to by any t G T = to + 



e^(s)- 1 /i(s- 1 , «(s)) - ds 



(5.33) 



sup tgT I e * 



< 7, and e* t Q z(t) 



(0, i) as t G T, |t| 
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OO. 



Let R = Rr],r be as in (5.17). If 7 > |C| and 77, r are sufficiently small, then the integral 



operator T~L in the previous paragraph defines a contraction in the complete space of all continuous 
functions z : R fl V — > C 2 , complex analytic in the interior of R, such that | e* t~ a \ \\z(t)\\ < 7 
for every t € R fl V, with respect to the metric (z\, Z2) h-> sup igra -,y | e* t~ a \ \\z\{t) — 22 (Oil- For 
every to £ R(lV the restriction of z to T := to + M>o is equal to the unique fixed point zq of the 
integral operator % : Z — > Z in the previous paragraph, and therefore zq ■ T — > C 2 extends to 



a solution of (5.25) on the much larger domain RnV, which is denoted by the same letter and 

3 i 



satisfies sup 4gra -,y | e* t a \ \\z(t)\\ < 7. 

Substituting an asymptotic expansion z(t) 



e -t f a f -j e . for t G s, |t| ^ 00 in the right 



hand side of z(t) = H(z)(t) yields a similar expansion with k replaced by k + 1, which procedure 
stops at k = N—l. In the induction step it is used that in £ we have ||.z(t)|| = 0(| e* t a \) = 0(\t\~ M 
for every M > 0. This leads to an asymptotic expansion y(t) ~ C e~* t a Ylj=o 



r 



3 dj, where the 



dj satisfy the equations do = (0, 1) and —dj + (a—j + l) dj-i = Yll=o L% dj-i for < j < N — l. For 
j = 1 we have — d\ + a do = Lq d\ + L\ do, which implies that (Li — a) (0, 1) is in the image C x {0} 
of Lq + 1, and we recover that a has to be equal to the right lower corner of L\. If j > 1, then the 
equation — dj + (a — j + 1) cfj-i = Lq dj + X^ =1 <ij-i determines dj only modulo the kernel {0} x C 



2 dj- 



=1 

of Lo + 1, but the equation with j replaced by j + 1 implies that (Li — a + j) dj + 
belongs to the image C x {0} of Lq + 1. As the right lower corner of L\ — a + j is equal to j ^ 0, the 
space (Li — a + j)({0} x C) is transversal to C x {0}, and it follows that dj is uniquely determined 
by the j-th equation and the (j + l)-st equation. Because we can take iV arbitrarily large, and the 
dj, < j < N — 1 do not depend on N, this implies (5.18), where (5.19 ) follows from a substitution 
f(t~\ y(t)), y(t) 



of (5.19) in y'(t) 



p(t)- Pt (t). 



mis completes the proof of the lemma in the upper half plane, where (5.22) follows from the 
aforementioned estimates for the integrals =t (s)~ 1 ^(s^ 1 , z(s)) =t ds. In turn (5.22), in com- 
bination with p^(t) = 0(t _1 ), implies that pc{t) is bounded away from zero on the part of the 
boundary of V Vt r in the upper half plane and sufficiently far away from the origin. The statements 
regarding the behavior in the lower half plane follow by replacing p^(t) by Pi(t). Note that in the 
lower half plane p^(t) is equal to one of the solutions Pt Q , a --> an d therefore p^it) extends from V to 



a domain of the form V U R 



r), r ■ 



□ 



Remark 5.7. Because p^it) is a solution of (5.3) as in Lemma 5.2, it follows from Lemma 5.6 that 
there is a unique constant S G C such that 



■Pt(t)~e-* t a (0, S) 



(5.34) 



as Ret/ log |i| — > 00, \t\ —> 00, where 5 = implies that Pi(t) = p-[(t) for all t £ R-q,r with r\ and 
1/r sufficiently small. Because this is a nonlinear analogue of the phenomena described by Stokes 
[26], we follow Costin [51 Th. 1] in calling S the Stokes constant pertaining to the comparison in the 
right half plane of the solutions Pi(t) and p-f of (5.3). Generically the Stokes constant is nonzero. 



This happens in particular for the Boutroux-Painleve system, see Lemma 5.14 below. 

Furthermore, for every solution p{t) = Pt 0: a- (0 °f (5.3) in Lemma 5.2 there are unique C-j-, Cj, G 
C such that 

p(t) - p t (t) ~ e - 1 t a (0, C t ) and p{t) - Pi (t) ~ e" 4 t a (0, Cl) (5.35) 

as Ret/ log |t| — > 00, |t| — > 00, where conversely the solution p(t) of ( |5.3[ ) is uniquely by any of the 
two asymptotic identities in (5.35). Because (p(t) — p^(t)) — (p(t) —pi(t)) = Pi(t) —p^{t), it follows 
that 

C t - C i = S. (5.36) 



Conversely, every complex number occurs as Of in (5.35) for a unique solution p{t) of (5.3) in 
Lemma |5.2| The following statements are equivalent: 



2!) 



i) Cj. = C+ for some solution p(t) of (5.3) in Lemma 5.2 

ii) 5 = 0. 

iii) Ci = for every solution p(i) of (5.3) in Lemma 
i y ) = Pf(0 f° r an i in a domain of the form i?^, 



5.2 



Recall the definition of \i[ as the upper left corner of the matrix L\ in (5.16). In analogy with 
(5.34), the substitution t i->- — t leads to the existence of a unique constant S- G C, the Stokes 



constant pertaining to the comparison of the solutions Pi(t) and Pf(t) in the left half plane, such 
that 

P^(t)-Pt(t)~e*^(5_,0) (5.37) 
as Ret/log |i| — > — oo, |t| — > oo. Here 5_ = implies that P\.[t) = p\{t) for all t G —R ri ,r with 
?7 and 1/r sufficiently small. We have S- = S = if and only if p+(t) and pi(i) have a common 
extension to a single- valued solution p(i) of (5.3) on a neighborhood of t = oo such that p(t) — > 
as t 



oo. 



Write r(t) = e - * t a , and let 



P formal (t) 



ft, «£Z>o 



r(t) h t-*p h)l 



(5.38) 



be a formal power series in r(t) and t , with coefficients p^ j € C 2 . The termwise derivative of the 



right hand side of (5.38), where d{r{t) h t~ 1 )/ dt = -hr(t) h t- 1 + {ha- i) r{t) h t' 1 ' 1 , is a formal 



power series in r(t) and t , which by de finiti on is the derivative of Pformai(^) with respect to t. 
Substitution of p formal (i) in the expansion (5.13) yields a formal power series in r(t) and t -1 , which 
by d efinition is u(t _1 , p formal (*))■ The formal solutions p(t) = p formal (t) of the differential equation 
(5.3) are obtained by equating the coefficients of r(t) h t~ l . 

Let T be an unbounded subdomain in the complex upper half plane of the domain R Vt r in Lemma 
5.6[ and assume that there exist strictly positive constants C%, C2, ei, £2 such that 



Ci\t\~ ei < \r(t)\ < C 2 \t\ 



for every t E T. In such a domain any formal series (5.38) is an asymptotic series. The proof of 



Lemma 5.6 yields that the solutions pc(t) in Lemma 5.6 have a formal series Pformai,c(i) as their 
asymptotic expansion in T, in the sense that for every N there exists an m such that pcif) — 

E 



h<m, i<m 



T(t) h t- Z p hti = 0{t- N ) as t G T, \t\ 



00. The formal series p formal, c{t) is a formal 
series solution of the differential equation (5.3) such that po, = and pni = (0, C). As will 
be verified in the proof of Lemma 5.8 below, a formal series solution (5.38) of ( |5.3| ) is uniquely 
determined by the initial consitions po, = 0, po, 1 = (0, C), and p^j = C h c^^i for uniquely 
determined universal coefficients j £ C 2 . This is the formal solution in O. and R. Costin (4), 
(5)] for n = 2 and C = (0, C), where the latter implies that only the terms with k = (0, h) appear 
in loc. cit. The proof of Lemma 5.8 also yields an independent verification of the existence of the 
formal series solution of (5.3) such that po 5 o = and po,i = (0, C). 
In the subdomains S in Lemma 



5.6 



the function r(t) is of smaller order than t~ l for every 
i G Z>o, when (5.38) only is an asymptotic series in the above sense if all terms with h > are 
deleted, and p formal, c(0 reduces to the right hand side of (5.11), which is independent of C. For 
Pc(t) ~ Pt(*) we have the asymptotic expansion (5.18) in S, where the right hand side is equal to 



( |5.38 ) with only the terms with h = 1 retained. 



If U is an unbounded subdomain of R Vt r in the upper half plane on which r(t) — > and t~ % = 
o(r(t)) as t G U, \t\ — > 00, then Pformai,c(i) is only an asymptotic series for t G U, \t\ — > 00 if 
all terms with % > are deleted, when pc(t) is asymptotically equal to this asymptotic series for 
t G U, \t\ —> 00. For t on the boundary of R<q, r and \t\ > r, the absolute value of r(t) is equal 
to the constant 77, p formal, c(t) is not an asymptotic series in the above sense, and we only have 
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the estimate (5.22). As T is a transitional region between the subdomains S on the one hand and 
the subdomains U and the boundary of R v<r on the other, the formal series p formal, c{t) and the 
asymptotic expansion pc(t) ~ Pform»i,c(i) as t G T, \t\ — > oo could be called the transitional series 
and the transitional expansion, respectively. 



Lemma 5.8 below, which follows from O. and R. Costin Th. 2(i) and Sec. 6.9], yields that 
for each i 6 Z>o the series (5.39) converges for small |r|. This allows to formulate the asymptotic 



expansion (5.49) for the solution pc(t) of (5.3). The domain of i's where the asymptotic expansion 
(5.49) holds extends well beyond the part R„ r of the boundary of R^^ r where Imt > 0, r(i) = rj, 
and \t\ > r. Along R' r , Lemma 



5.6 



yielded that pc(t) remains bounded away from zero, but did 
not provide an asymptotic expansion for \t\ —> oo. 



Lemma 5.8. Let C 6 C. Then there is a unique formal solution (5.38) of (5.3) such thatpo : o = 
and pi t o = (0, C). For each i G Z>o the series 



F l (r)=J2r h Ph, 



h=0 



converges for r in a neighborhood of in C, where the complex analytic functions Fi satisfy 

dF (r) 



— T 



dr 



«0(F (r)) := u(0, F (r)) and 



dFi(r) 



for i G 



~ T dr 
7 , >0 - Here Pi(F , 



d F- i (t) 

Pi(F (r), Ei(r)) - a r ^ J + (i 



l)^-i(r) 



(5.39) 

(5.40) 
(5.41) 



Pi(Fo, ■ ■ ■ , Fi) 



Fi + Qi^Oi • • • 



(5.42) 



Fi) denotes the coefficient of t 1 in the expansion of 

in nonnegative integral powers of t~ l ; a finite sum of weighted homogeneous polynomials of degree 
< i in the Fj with < j < i, where each Fa has the weight j . In particular 

dv(0, F p) 
dF 

where Q%{Fq, . . . , Fi— i) is a similar polynomial in the Fj with < j < % — 1. We have 

Fi{0) = P0 ,i = a, (5.43) 

where the c% are the coefficients in \ 5.11]) . 

Conversely, the system (5.4-0), (5.4-1) has a unique solution Fi{r) which is complex analytic on 
a neighborhood of r = in C such that Fq(0) = and Fq(0) = (0, C). If j G Z>o then, given the 
functions Fk(r) for < k < j — 1, the function Fj(r) is uniquely determined by the conditions that 
it is a complex analytic solution on a neighborhood of t = of the equation {5.41 ) for i = j, and 
that the equa tion (5.41 ) for i = j + 1 admits a complex analytic solution Fj + i(r) on a neighborhood 



of t = 0. If (5.39) denotes the power series expansion of Fi(r), then (5.38) is the formal solution 
of (5.3\ ) in the previous paragraph. 

If Fi,c{ T ) denotes the solution of (5.40), (5.41) such that -Fo,c(0) = and -Fo,c'(0) = (0, C), 
then Fi c(t) = Fi \{C r). Let Tc be the maximal domain of definition of the, possibly multi-valued, 
solution Fq = Fq^c of (5.40). Then 7o = C, Tc = C~ l T\ if C ^ 0, and for each i G Z>o the 



function Fi^c extends to a, possibly multi-valued, solution of (5.41 ) on Tc- 

r s , F (t) 



Proof. With the substitutions r 
autonomous limit system 



dpp(s) 
ds 



= po(s), the equation (5.40) is equivalent to the 
uoCpo(s)) (5.44) 



of the system (5.3). Because t>o(0) = v(0, 0) = and Uq'(0) = Lq, an application of Lemma 5.6 
with v replaced by vq, when ci = 0, L% = 0, and a = 0, leads to Po,t( s ) = Po,i( s ) = an d a 
unique solution po(s) of (5.44) on Res > a » such that po(s) ~ e~ s (0, C) as Res > a » 0, 
|s| —7- oo. The corresponding function Fq(t) = po(s) is complex analytic on the punctured disc in 
the complex plane determined by the inequalities < |r| < e := e~ a , and satisfies Fq(t) ~ r (0, C) 
as t — > 0. It therefore follows from the theorem on removable singularities that Fo(r) extends to 
a unique co mple x analytic function on the disc |r| < e, denoted again by Fq(t), where Fq(t) is a 
solution of (5.40) on the disc |t| < e s uch that Fo(0) = and Fo'(0) = (0, C). As the formal power 
series solution Ylh>o rh Ph,0 of (5.40) with po, o = 0> Pi,o = (0, C) is unique, the above complex 



analytic solution Fq(t) is unique as well. 

Let i G Z>q. For given complex analytic functions Fj, < j < i 



1, in a neighborhood of the 



origin, (5.41) is a linear inhomogeneous differential equations of the form 

dK; 



dr 



A(r) Fj + Gi(r), 



(5.45) 



where A(r) := dv(0, p) / 'dp\ p= p ^ and Gj(r) are complex analytic functions of r in a neighborhood 



of r = 0. With the substitutions A(r) = Ylk>o T A.k> where Ao = Lq as in (5.21), and Gj(r) 
Sfe>o rh Gi,k, the formal power series Fj = ^fc>o 
and only if 



Fj & is a solution of the differential equation if 



i, k 



1=1 



for every k G Z>o- These equations determine the coefficients Fj fe in terms of the Fij with 
j < k and Gi,fc, with the exception of the equation for k = 1, where the + part yields F^ = 
— (1/2) (AiF^o + Gj i i) + , but the resonance in the — part yields no equation for Ff x , but instead 
the solvability condition (Ai Fj^ + Gi t i) _ = 0. As F^o = — Lq~ G^o, this solvability condition 



is equivalent to the equation G il = (Ai Lq 1 Gj i o) for the inhomogeneous term Gj(r) in the 



Qi{Fo,o, 



Fj-i,o) + (i — l)Fj_i i o is determined, the 



differential equation. Because G^o 
solvability condition determines G~ 1 . 

For i = 1 we have Fi(F , Fi) = div{0, F ) + a 2 v(0, F )Fi, hence G?i(t) = 9iu(0, F (r)) - 
3iu(0, Fq(t)) + ai)(0, Fq(t)). Therefore G^o = 9iv(0, 0) and G\,\ is equal to C times 



d 2 v{u, p) 



dudp 



a 



dv(0, p) 



dp~ 



p=0 



u=0,p=0 

because Fq(0) is equal to C times the second basis vector. As the minus part of the second term 
is equal to —a, it follows from the definition of a in Lemma 5.6 that G^ 1 is equal to C times the 
lower left corner of c?2 2 u(0, 0) Lq^ 1 d\v(Q, 0). On the other hand Ai is equal to C times &2 2 v (0, 0) 
applied to the second basis vector, when the symmetry of the second order partial derivatives yields 
the solvability condition G^ l = (A\ Lq^ 1 Gi t o)~ . 

For i > 2 the above computations yield that G~ 1 depends in an inhomogeneous linear way on 
Fj~ i 1 with coefficient equal to i — 1 ^ 0. Therefore the solvability condition is satisfied by a unique 
choice of F^_ 1 1 . It follows that the system (5.40), (5.41) has a unique formal solution such that 
Fo(0) = and Fq(0) = (0, C), where for the unique determination of Fj(t), j > 0, one needs (5.41 ) 
for % = j and i = j + 1. 

If in a neighborhood of r = the functions Fj(t) are complex analytic solutions of (5.40), (5.41 ) 
for < j < i — 1, then (5.45) is an inhomogeneous linear differential equation with complex analytic 
coefficients A(r) and Gj(r), with r = as a regular singular point. Therefore every formal powers 
series solution of (5.45) is convergent. For inhomogeneous linear differential equations this theorem, 
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which sounds classical, seems to be due to Gerard and Levelt [1 H Lemme 4.2]. For nonlinear higher 
order scalar ordinary differential equations, it has been obtained by Malgrange \20\ Remarque 4.1], 
where the proof also works for nonlinear systems near a regular singular point. 

If g {j) = /(Cr), then r 5 '(r) = (Cr)f'(Cr). Therefore the function $ : t t-> F ,i(Cr) 
satisfies ( |5.40[ ) with $ (0) = 0, $ '(0) = C, hence F 0> c(t) = $ ( r ) = F i(Cr ). Furthermore, 
if Fj(r) = F^i(Ct) for < j < t - 1, then : t A ^^Ct) satisfies ( |5.41[ ), hence F;(r) = 
$i(r) = F^i(Ct). This proves F^c^r) = ^^(Ct) by induction on i If the coefficients p/^j for 
C = 1 are denoted by c/^j, then the equations Fi t c(r) = F i(C r) are equivalent to the equations 
mentioned in the text preceding Lemma 



5.8 



The last statement in Lemma 5.8 follows, by induction on i, from the description of F{ as a 



solution of an inhomogeneous linear differential equation (5.45), of which the coefficients A(r) and 
Gi(r) are complex analytic functions of r G Tc- 



Remark 5.9. It would have been more precise to define the domain Tc as the Riemann surface 
of the maximal solution of (5.40) such that Fq(0) = 0, -Fo'(O) = (0, C), as follows. Let Dq := 
{p G C 2 | (0, p) G -D}, where -D is the domain of definition of the vector field v in (5.3). In 
M := C x Dq \ {(r, p) G C x Dq I r = and u(0, p) = 0} we have the regular complex one- 
-dimensional distribution D defined by the equation t dp + v(0, p) dr = 0. For the local solution 
Fq of ( |5.40| ) with F (0) = and F '(0) = (0, C), the set J c>e := {(r, F (r)) | < |r| < e} is an 
integral manifold of T>; let Iq denote the maximal integral manifold of T> which contains Ic e- Then 
Tc is canonically identified with the Riemann surface Ic U {(0, 0)}. The inverse of the projection 
Ic 3 (t, p) — > t followed by the projection (r, p) i->- p is the, possibly multi-valued, maximal 
solution i 7 ^ c °f (5.40) mentioned in Lemma 



5.8 



Lemma 5.11 below follows from O. and R. Costin Th. 2(ii)]. Our proof is different. Lemma 



5.11 implies all the previous asymptotic expansions in Section |5,1[ but not the explicit estimates 
such as in Lemmal5.2landl5.31 



The following Lemma 5.10 will be used in the last part of the proof of Lemma 5.11 and in 
Lemma 5.18| about the poles of the truncated solutions. Lemma |5.10| is a detailed version of the 
estimate in (24)]. 



A 



= t. Then the 
argt| < ir, and 



Lemma 5.10. Let r G C \ {0} and logr a given solution A of the equation e 
solutions t of the equation e _i t a = r, where t a = e a logt ; logt = log \ t\ + i argt, 
\t\ is large, form two sequences t n , where n G Z, n >> and n « respectively, such that 

t n = 2tt i n + a log(27r i n) — log r + s(u, w) (5.46) 

as \n\ —> oo. Here log(27rin) := log(27r|n|) + i sgn(n)7r/2, 

1 log(27rin) 



w 



-u logr lav, 



(5.47) 



27rin 27rin 

and s(u, w) is a convergent power series in (u, w). More precisely, s(u, w) is a complex analytic 
function of (u, w) in a neighborhood of (u, w) = (0, 0), equal to the unique small solution s of the 
equation 

s = f(u, w, s) := a log(l + w + a u log(l + us + w)) (5.48) 
for u and w both small. It follows that, with the substitutions (5.-^7) and A := logr, 



s(u, w) 



aw + a 2 uw — a w 2 /2 + a 3 u 2 w — 3 a 2 u w 2 /2 + a u> 3 /3 + 0(n 4 ) 
—a A n + a 2 v — a A (a + A/2) u 2 + a 2 (a + A) ui> — a 3 w 2 /2 
-a A (a 2 + 3 a A/2 + A 2 /3) u 3 + a 2 (a 2 + 3 a A + A 2 ) u 2 v 
-a 3 (3 a/2 + A) u v 2 + a 4 w 3 /3 + 0(n" 4 ). 
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The t n depend in a complex analytic way on logr, r G C \ {0}. If r runs around the origin 
once in the positive direction, then t n (r) moves continuously to i n _i(r). In this sense the r n (r) for 
n » and n « can be viewed as two multi-valued complex analytic functions of r. 

Proof. The equation e - * t a = r is equivalent to the equation t — a log t = 2ir i n — log r for some 
n G Z. For large log |t| = o(|i|), |n| is large, 2irin/t = 1 + o(l), and i = 2-7rin (1 + o(l)). 

With the notation a := 2nin — logr, the equation for t is equivalent to = 0(t) := (t — 
a) (1 — (t — a)~ l a logi). If 7 : [0, 1] — > C \ {a} is a Jordan curve which winds once around 
a in the positive direction and \a\ | log s| < \s — a\ for all s G 7([0, 1]), then the the winding 
number = (2tt)~ 1 times the increase of the argument of (fio 7 a is equal to one, hence the function 
<j) has a unique zero in the interior of 7, which is simple, and therefore in view of the implicit 
function theorem depends in a complex analytic way on a. If \a\ > 1 and |s — a\ = \a\ — 1, then 
1 < I s I > 1, \a\ I logs| < \a\ (7r + log |s| +ir) < \a\ (7r + log(2|a| — 1)) and therefore \a\ | logs| < \t — a\ 
for all s on the circle 7 around a with radius \a\ — 1 if |q| (it + log(2|a| — 1)) < |a| — 1, which 
happens if \a\ is sufficiently large. Because |n| >> implies that \a\ = |2-7rin — logr| >> 0, and 
t — a = alogt = a log(2-7rin) + o(l) = a log a + o(l) = o(|a|), the conclusion is that every solution 
t of the equation t — a logt = 27rin — logr such that t = 2?rin (1 + o(l)) is equal to the unique 
zero of 4> in the disc around a with radius \a\ — 1. Therefore the solution t = t n (r) of the equation 
t — a logt = 2-7rin — logr with | argt| < n and \t\ >> is unique, depends in a complex analytic 
way on r, and satisfies t = 27rin(l + o(l)). 

If we substitute the latter estimate for t in the right hand side of t = 2ir i n — log r + a log t = 
2irin — logr + a (log(27rin) + log(l + (27rin) _1 (— logr + a logt)), we obtain that t = 2irin + 
a log(27rin) — logr + s with s = 0((log |n|)/|n|) — > as \n\ — > 00. With the definitions of u, v, 
and w in (5.47), s satifies the equation (5.48). Because /(0, 0, s) = 0, it follows from the implicit 
function theorem in the complex analytic setting that there exist open neighborhoods U, W, and 
S of the origin in C such that for each (u, w) € U x W the equation (5.48) has a unique solution 
s = s(u, w) £ S, and that (u, w) 1— > s(u, w) is a complex analytic function on U x W, with 
s(0, 0) = 0. Substitution of (5.47) in s = s(u, w) leads to (5.46). Subsequent differentiations of 
s(u, w) = f(u, w, s(u, w)) with respect to (u, w) at (n, w) = (0, 0) up to the order three and 
subsequent evaluations of the results at (n, w) = (0, 0) lead to the asymptotic formulas for t n {r) 
modulo terms of order n -4 . With the help of a formula manipulation computer program one could 
go on a little bit further, but the complexity of the asymptotic formulas increases rapidly with 
growing order. 



log r moves continuously 

□ 



If r runs around the origin once in the positive direction, then 2tt i n - 
to 2-7ri(n — 1) — logr, and therefore i n (r) moves continuously to i n _i(r) 

Lemma 5.11. Let K be a compact subset of C 2 such that {0} x K is contained in the domain 
of definition D of the vector field v in (5.3). Let To be a bounded subdomain of the domain T in 
Lemma 5.8, such that Fq(t) £ K for every r G 7o- Let V be a t-domain in the complex upper 
half plane where \t\ >> and r(t) £ 7o- Then the solution pc(t) in Lemma 5.6 extends to a, 
possibly multi-valued, solution p{i) of (5.3) on V . This solution has the asymptotic expansion 



p(t) ~ ^j>o t 1 Fi( T {t)) i n the sense that, for every m G Z>q, 



P (t) = J2 r i Fi(r(i)) + o(r 



(5.49) 



ast£V, \t\ - 
Proof. Write 



00. 



m—l 



] (*):= E t-'FMt)), Po (*):=E*~ ijF H°)' and S(t) ■= P(t) ~ P H (t)- (5-50) 



i=0 



34 



According to Lemma 5.6, p(t)—p^(t) ~ C 'r(t) asr(t) — > 0. The expansion (5.18) and equation (5.43) 



imply that n (t)-pP(t) = 0(t~ m ). Finally p^(t)-pP(t) = F o (r(t))-F Q (0)+^ti t~ l (F(r(t)) 
Fi(0)) ~ Cr(t) + 0(r 1 r(t)) ~ Cr(t) as r(t) -> 0, because F (r) ~ Ct as r -> and r(t) -> 

implies that t -> 0. Therefore yN(t ) := - pt(*)) - (p [m] (*) -P(j mI (*)) = °( r 0))' as T ( t ) ~> 
and <J(t) = i/H(f) + 0(t- m ), hence jjUgjt 4» yN(t) = 0(r m ). 
The differential equations (5.40) and (5.41) imply that 



dpM(i) 
dt 



m—l 



;(0, F (r(t))) + £ r< [L F(r(t)) + Pi{F {r(t)), . . . , iVi(r(t)))] 
dF„ 



+t~ 



a r- 



i=l 

-i(r) 



dr 



m 



IF, 



m—l 



(r) 



r=r(t) 



hence 



dyH(t) 



dt 



-m— 1 \ 



d P (t) dpH(t) d^^-^kt)) 

dt dt dt 

= ^t" 1 , p(t)) - (v(t-\ P H(*)) + 0(t" m )) + 0(t- 
= ^t" 1 , pH( t ) + y H( t) ) _ ^-i, p N (t)) + 0(t -m } 

as long as r(t) G 7o, where in the last identity we have used that p(t) = p' m '(t) + y' m '(^) + (Pt(*) 
p H(t)) = p H(t) + y H(t) + 0(t~ m ). 



We start with the proof for t in a domain as in (5.17) with 77 > sufficiently small, when 
( |5.49| ) is [3 (20)] in our situation. It follows that dy^(t)/dt = L(t, p [m] (t)) y H ft)+0(||y [m] (t)|| 2 ) + 
0(t" m ), where L(t, p) := fluff 1 , p)/dp. Because pM(t) = ^(t) - (p t (t) -p[, m] (t)) +pM(t) - 

PcT 1 (*)) = Pt(*) + 0(t~ m ) + 0(r(t)), we have L{t, pM(t)) = L(t) + 0(f m ) + 0(r(t)) with L(t) as 
in ( |5.20 ), hence 



dyH(t) 
dt 



L(t)y [ra] (*) + 0(|r(t)|||yN( t )||) + o(||y[« 



! ) + 0(t" 



Applying a linear substitution of variables y = A(t) z as in the proof of Lemma 5.6 we arrive at an 
integral e quatio n z = K(z) for z : t i-> A(t)" 1 yH (t), where ~H(z)(t) + = H{z)(t) + and U{z)(t)- 

in view of z(t) = A(t)~V m] (*) = o(r(t)), and 
zj = U(jr(tj| ||z||j + 0(||z|| 2 ) + Off m ). If 7/ and 1/r are sufficiently small, then Ti is a 
contraction on the space i? of all continuous functions z : R^r — > C 2 that are complex analytic on 



are as in (5.30) and (5.33), respectively, C 

h(t~\z) = 0(|r(t)|||z||) + 0(||z|| 2 ) + Or 
contraction on the space Z of all continuou 
the interior of Rn >r and satisfy a uniform bound ||z(t)|| < C\t\~ m , where % is a contraction with 
respect to the metric (z\, Z2) ^ sn PteR v r 1*1™ ll z i(^) ~~ z 2(t)||. It follows that H has a unique fixed 
point z in Z, when y [m] (t) = A(t) z(t) satisfies y [m] (t) = Off m ) when t £ R Vt r and |t| -> 00. This 
completes the proof for t in a domain where |r(t)| remains sufficiently small. 

For t in a subdomain of V where 7~(t) := e~* t a £ To remains bounded away from zero, we use 
r = r(t) instead of t as the independent variable. The domain 7o can be arranged such that its 
points t can be joined with points tq close to the origin by smooth paths 7 in To, parametrized 
by arclength, with a uniformly bounded length, and staying at a uniform distance away from the 
origin. The solution y[ m '(r) will be estimated along such 7. With r along 7, we view t = t n (r), 
n >> 0, as the multi-valued inverse in Lemma 5.10 of the function t h-> r(t) with large |t| and 
< argt < 7T. Then t = t n (r) = 27rin (1 + o(l)), uniformly for r along 7. Because 



dyt 



dr 



-1 + a/ty 



-1 dyt" 



dt 



Of 



\y L 



+ 0(\Ut) 



3.1 



we have 



d||yM( 7 ( s ) 
ds" 



< 



dyH( 7 ( s )) 



ds 



< 



(y lm] )'(l(s)) < C x ||y H (7(^))ll +C 2 M 7 (s))|- 



for some positive constants C\ and C<i- Or, 



d||y H (7(*))ll 
ds 



i \w 



for a non- negative continuous function r(s) 



^( 7 ( s ))\\ + C 2 \tMs))\- m -r(s) 
Lagrange's variation of constants formula yields 



|yH( 7(s) 



ClS \\y [m] (7(0))\\+ / e c ^ s ^ (C 2 \tMa))\~ m -r(a)) da 



< 



Cls ||y W (7(0))||+C 2 / e 01 ^) |t„(7(<r))|- m da. 



Write t = t n ( 7 (s)), where |t| >> 0. Then, uniformly for < a < s, i n ( 7 (cr)) = 27rin(l + o(l)) = 
t(l+o(l)), hence \t n {l{<?))\~ m = 0(\t\~ m ). Furthermore, because 7 (0) = t is small, ||y H (7(0))|| = 
O(|tn(7(0))r m ) = 0{\t\~ m ). Combination of the estimates yields y [m] (j(s)) = 0(t" m ). That is, 
returning to t as the independent variable, y^ m \t) = 0(t~ m ). □ 

Lemma |5.11| can be used in order to obtain asymptotic information about the large t where the 
solution pc(t) becomes singular, where r(t) is close to points r where the solution Fq(t) of (5.40) 



becomes singular. This will be done in more detail for the Boutroux-Painleve system in Subsection 

r with 



5.2 leading to the asymptotic description in Lemma 5.18 of the poles in the domain i?^, 
large rj and r correspondingly large. These poles correspond to the first sequence of poles of the 
truncated and triply truncated solutions of the first Painleve equation which appear beyond the 
boundary of the truncated domains. 

5.2. Truncated and triply truncated solutions. In this subsection we collect the conclusions 
from Subsection |5.1| for the Painleve equation. 



Boutroux [2, §13] found a family of solutions u(z) of the Boutroux-Painleve system (2.2), depend- 



ing in a complex an alyt ic way on one complex variable, such that u(z) converges to the equilibrium 
point (ei\/6, 0) of (5.1) as t = A+ z runs to infinity in the complex plane in the direction of the 
positive real axis. On [21 p. 346] , he called these solutions truncated in the direction of the positive 
real axis. These solutions correspond to the solutions p to a -(t) and pi(t) described in the lemmas 
ET21 and 1531 

In [21 §14], Boutroux found a solution u(z) which converges to the equilibrium point of (5.1) 
as t = X + z runs to infinity in a sector in the complex plane which does not contain the positive 
imaginary axis. Because such a sector contains the three orthogonal axes equal to the negative 
real, negative imaginary, and positive real one, he called this solution the triply truncated solution. 
This solution corresponds to the solution p±(t) described in the lemmas 5.3, 5.5, and 5.6 In the 
sequel we will follow Boutroux's terminology. 

Remark 5.12. The proofs in Boutroux [21 §13, 14] contain gaps and errors. In the formula for Y\ 
following [2, (51)] the factor X" 1 in the first and second integral should be replaced by —X~ 2 /y/V2 
and X~ 2 /\/l2) respectively. A similar correction is needed in [21 (54 bis)]. In the right hand side 
of [21 (54)] the term 6 £Q Y k Y^ k is missing. More seriously, the inductive assumption that 
Y k (X) — > as X — > 00 for every k < j implies that the second limit relation at the bottom of [21 
p. 343] is automatically satisfied, when the first one is equivalent to the condition that Yj(X) — > 

Therefore 



12X 



as X — > 00. This determines Yj only up to the addition of a constant times e 
the estimate lj < Klj-\X a on [21 p. 345] cannot be proved for the solutions Yj described by 
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Boutroux. It is a bit surprising that Boutroux did not use the method of Cotton |8] , which appeared 
two years earlier in the same journal as the paper of Boutroux. 



Remark 5.13. On [2, p. 261], Boutroux gave a definition of truncated solutions which looks quite 
different from the definition in terms of the convergence to the equilibrium point of the limit 
system. In our notation, we understand his definition on [21 p. 261] in the following way. For every 
zq E C\{0} and c E C there is a unique solution U(z) = U ZQ)C (z) of (2.4) which has a pole at z = zq 

l4 in its Laurent expansion at z = zq. The parameter c corresponds 



and c as the coefficient of (z — zq) 
in an bijective affine manner to the position on the pole line Lg \ in Okamoto's space, see (3.2 ). 
Assume that u ZQtC (z) has another pole at z = Z\ ^ zq. Because the vector field is regular, nonzero 
at and transversal to the pole line, an application of the implicit function theorem yields that for 
every (zq', d) near (zq, c) there is a unique Z\ = Z\{zq ', c') near z\ such that the solution u Zo ' )C '(z) 
has a pole at z = z±, and the function (zq', d) h-> Z\{zq , c') is complex analytic. It seems that 
Boutroux considered this to be evident, as he did not provide further proof. Then, according to 
our interpretation of (2j p. 261], u Zo " jC "(z) is a truncated solution if the function Z\{zq , c') has a 
complex analytic continuation for (zo',c') on a path approaching (zq", c"), but that the superior 
limit of \Zi(zo', c')\ is infinite if (zq\ c') — > (zo", c") along the path. We find this an interesting 
definition of truncated solutions. However, we did not find any statement or proof in Boutroux's 
paper which relates the truncated solutions as defined on [21 p. 261] with the truncated solutions 
defined in term of their convergence to one of the equilibrium points of the limit system. We will 
stay with the second definition, and view it as a challenge to find relations with the first one. 



According to Remark 2.1 the substitutions x = —2 3 / 5 3 x / 5 £, y(x) = 2 4 / 5 3 3 / 5 turn the 



first Painleve equation (1.1) into 

d 2 v /de 

((5/4) t) 4 / 5 , r,(0 -- 



- (v 2 - 6/2, 

e /2 Mt), and rf{£) 



(5.51) 

£ 3 / 4 7r 2 (t) turn the corre- 



when the substitutions £ 
sponding first order system for (77, 77') into 

dvri/dt = 7T 2 -2(5£)~ 

d7T 2 /dt = (7Ti 2 - l)/2 - 6 [i)t) "7T2. 

At the equilibrium point (tti, 1T2) = (1, 0) of the limit system for t = 00, the linearization has the 
eigenvalues ±1, and the substitution of variables 7ri = 1 +p + +p~ , 7r 2 = p + — p~ lead to the system 
(5.3) with 



tti, 

-3(5*)- 



(5.52) 



v(i/t, P y 
v(i/t, P y 



p 



i/5t-p + /2t + p-/m + ( P + +p-y/4, 

- l/5t-p~/2t + p + /10t - (p + +p~) 2 /4. 



It follows that the vector c\ in (5.14) and the matrix L\ in (5.16) are equal to 

1/10 



1/5 
-1/5 



and L\ 



-1/2 
1/10 



-1/2 



respectively. Because c{ + c{ =0, the asymptotic expansion (5.11) implies that for the truncated 
solutions p(t) we have vri(t) = l + 0(i~ 2 ) as t E S, \t\ — > 00. Furthermore, both the left upper corner 
and the right lower corner a of L\ are equal to —1/2, and therefore p{t) — p^(t) is asymptotically 
equal to r(t) = e~* t" 1 / 2 times a series in nonnegative powers of t -1 as t E S, \t\ — > 00, see (5.18). 
The fact that Re a < implies that the large positive and negative parts of the imaginary axis are 
contained in the interior of the domains R v , r in (5.11). More precisely, the part in \t\ > r of the 
boundary of R V) r is given by the equation 

cos(argt) = (— (1/2) log \t\ — log 77) /|i| , \t\ > r, — 7r<argt<7r, (5.53) 

which is to the left of the imaginary axis if and only if |t| > r/~ 2 . See Figure 5.2 
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Because the Painleve property implies that all solutions y(x) of (1.1) are single-valued, the 
analytic continuation of (7i"i(t), ^(i)) along the path t e 10 , where fl£R runs from to 5/4 times 
2n, applied to the substitutions 



x = -2- 3 / 5 3- 1 /5(5t/4) 4 /5, 2/ (x) = 2- 4 / 5 3- 3 / 5 (-2 3 / 5 3 1 / 5 x) 1 / 2 ^ 1 (t) 



(5.54) 



leaves the solution y(x) of (1.1) invariant. We conclude as in (2.6) that the analytic continuation 
of (7Ti(i), 7T2(i)) along the aforementioned path is equal to (— 7Ti(t), 1 7T2 (t) ) - This agrees with the 
symmetry (t, tti, 7^) 1— > (it, —tt\, itv-z) of order four of the system ( |5.52 ). This analytic continuation 
converges for t running to infinity in the direction of the positive imaginary axis to the other 
equilibrium point (jri, ^2) = (—1, 0) of the limit system of (5.52) for t = 00. Applying this analytic 
continuation a second and third time, one obtains a solution which converges to (tti, H2) = (1, 0) 
and (-7T1, 7T2) = ( — 1, 0) for t running to infinity in the direction of the negative real and imaginary 
axis, respectively. In this way each of the four types of truncated solutions, for t running to infinity 
in the direction of the positive and negative real axis with (iri(t), ^(i)) converging to (1, 0), and for 
t running to infinity in the direction of the positive and negative imaginary axis with (iri(t), ^(t)) 
converging to (—1, 0), are obtained by analytic continuation from the solutions p t -(t) = Pi(t) in 
Lemma 15.21 and 15.61 

Similarly, the analytic continuation of p^(t) along the path t e l9 is equal to Pi(t) if 6 £ R runs 
from to 5ir. As this agrees with the symmetry (t, wi, -K2) h-» (— t, 7Tl, — ^2) of the system (5.52), it 
follows that Pi(t) + = p^(—t)~ and Pi(t)~ = p^(—t) + . If 6 runs from to 57r/2 and 157r/2, then the 
analytic continuation of p^(t) is equal to the two triply truncated solutions near (7Q, 7^) = (—1, 0). 
In this way each of the four triply truncated solutions, for t running to infinity in the upper and 
lower half plane with (7i"i(i), ^(i)) converging to (1, 0) and for t running to infinity in the left and 
right half plane with (iri(t), 7^(2) ) converging to (—1, 0), is obtained by analytic continuation from 
the solution p+(t) in Lemma 5.3, 5.5, and 5.6 

Because the system (5.52) is real, it has the symmetry (t, m, 1T2) h-> (t, tTT, P2 - )- Therefore, if 
(ni(t), 7T2(t)) is a truncated solution of (5.52) in the sense that it converges to (1, 0) a s t ru ns to 



infinity in the direction of the positive real axis, then t h-> (7Ti(t), ^(t)) is a solution of (5.52) with 
the same limit behavior, and therefore is truncated in the same way. The triply truncated solutions 
satisfy pi(t) = P[(t), which in combination wit h p±( t) = p^(—t) T implies that P[(t) = p^(—t)T. 



The corresponding triply truncated solution of (5.52) satisfies ni(t) = ni(—t) and %2(t) = —iT2(—t). 



Lemma 5.14. The Stokes constant S in (5.34) ^ nonzero and purely imaginary 



Proof. Because p\.(t) = Pt(t), the asymptotic identity for t G M, t — > 00 implies that S is purely 
imaginary. If S = 0, then Pi(t) = p^(t) in the right half plane, and therefore Pf (i) and Pi(t) would 
have a common extension to a small solution p(t) for — (3/2)7T < arg(t) < (3/2)7r, |t| > r, when the 
corresponding solution y(x) would be bounded by a constant times \x\ 1 ^ 2 for — (6/5)7r < arg(x) < 
(6/5)tt and |x| > R for some R. As this implies that the single valued function y(x) has no poles 
for large it would follow that y(x) has only finitely many poles, in contradiction with Corollary 

113 □ 

Remark 5.15. Because r(i) = e _i i -1 / 2 , the positive and negative imaginary axis are transitional 



domains as in the paragraphs preceding Lemma |5.8| the Stokes constant S is equal to the one in 
Costin [6l (2.8)]. According to P, Note (3) on p. 7], S = i \J 6/57T. The formula Pi(t) — Pf(t) ~ 
i ^/6/5^T e _t t -1 / 2 as — 7r/2 < argi < 7r/2, |t| — > 00 agrees, at least up to the sign, with Kapaev 
PH Cor. 2.4]. 



Figure 5.3 shows a region in the x-plane, the plane of definition of the solution of (1.1), where 
the corresponding truncated solution of (2.2) is close to one of the two equilibrium points of the 
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limit system (5.1). This region is the image under the mapping t \- > —2 3 / 5 3 1//5 ((5/4) t) 4 / 5 of a 



r), r ■ 



domain of the form R 

The triply truncated solution p^(t) of (2.2) is defined on a domain determined by the inequalities 
\t\ > r and 

— arccos((— (1/2) log |t| — \ogn)/\t\) < argt < 7r + arccos((— (1/2) log \t\ — \ogrj)/\t\). 

Here r and r/ are sufficiently large and small strictly positive real numbers, and arccos denotes 
the inverse of the bijective function cos : [0, vr] i — >■ [ — 1, 1]. As the inequalities for \t\ and aigt 
allow points t, t' such that \t'\ = \t\ and argf = argt + 2ir, the function p^(t) is interpreted as 

multi-valued. The properties Pt(*) ± = Pt(~ (5.34), and Pi(t) = 0(t _1 ) when t runs to infinity 
in the direction of the negative imaginary axis imply that for argi = —tt/2 and aigt = (3/2) tt we 
have p-f-(i) ~ e~* t~ l l 2 (0, —S) and p^{t) ~ e* t~ l l 2 (— i S, 0) as |£| — > oo, respectively. As Lemma 



5.14 



implies that S ^ 0, i/ie too branches of p^(t) do not coincide on the overlap. 
The image under the mapping t i— )• x = — 2~ 3 / 5 3~ 1//5 ((5/4) t) 4 / 5 of the aforementioned domain 
where p^(t) is small is a domain in the x-plane where \x\ is large and arg(x) runs from slightly 
smaller than (3/5) n to slightly larger than (11/5) tt. The other truncated and triply truncated 
regions are obtained from these by applying a rotation in the x-plane over &27r/5, 1 < k < 4. 



Because the truncated and triply truncated solutions of (2.2) are bounded in their domains of 



definition, they have no pole there, and therefore the corresponding truncated and triply truncated 



solutions of (1.1) have no poles in the aforementioned truncated and and triply truncated regions 



in the x -plane. 

Figure |5.4| shows the unique triply truncated region in the x-plane which is invariant under 
complex conjugation. If y(x) denote s th e corresponding triply truncated solution of (1.1 
i— > yix) is a solution of (1.1) which is triply truncated in the same domain. 



function 

uniqueness of triply truncated solutions of (2.2) leads to the following observation of Joshi and 
Kitaev [H Cor. 3]. 



Lemma 5.16. Let D denote the triply truncated domain in Fig. 5.4 which is invariant under 
complex conjugation. Then the solution y(x) of (1.1) which is triply truncated on D is real in the 
sense that y(x) = yix) for every x G C. In particular y(x) G R for every i£R not equal to a pole 
point of y. 



The next lemma discusses what happens with Lemma 5.11 in the case of the Boutroux-Painleve 
system. Our results correspond to [7, (113)] with the formulas for Hq(^) and Hi(£) on p. 
38], as vri(x) = Y(x) = 1 - 4/25x 2 + h(x), see [71 p. 36], hence 7Q i0 (£) = iZb(0 + 1, ^1,2(6 = 
H 2 (0 - 4/25 and 7Ti,j(£) = H$ for all I <£ {0, 2}. O. and R. Costin p. 39] wrote: "We omit 
the straightforward but quite lengthy inductive proof that all are rational functions of And 



on p. 40: "For large £ induction shows that H n 
of the proof. 



Const n £ r 



but did not give further details 



Lemma 5.17. With the notation of Lemma 5.11, the solution (iri(t), ^(i)) of (5.52) corresponding 
to pc{t) has the asymptotic expansion 



m—l 



*k(t)= t~ l 7T kJ (Cr(t)) + 0(t- 



(5.55) 



1=0 



astGV, \t\ 00. Here 7^(0 = (£ - 12)"^ 2 Pi, and 7r 2 ,j(f) = (£ - 12)"^ 3 P 2 ,f(0> where 
Pi,i(0 an d P2,i{Q is a polynomial in £ := C t of degree < 21 + 2 and < 21 + 3, respectively. We 



:!9 



have 



Pi, 0(6 = (£-12) 2 + 144£, (5.56) 

P 2)0 (0 = 144^ + 12), (5.57) 

Pi,i(0 = C(216 + 210^ + 3e 2 -e 3 /60), (5.58) 

P 2) i(£) = (497664 - 134784 £ + 266112 £ 2 + 25704 £ 3 - 24 £ 4 + £ 5 )/60. (5.59) 



Proof. It follows from the last paragraph in Lemma 5. 
C 



that it suffices to prove all the formulas for 



1, when t = £. 

The shortest proof of the formulas for 7Tl o(£) and 7T2,o(£) is to verify that these functions satisfy 
the differential equations — £ dP^o/d^ = P2,o> — £ dP2,o/d£ = (Pi,o 2 — l)/2 corresponding to 
(5.40), with the initial conditions Pi i o(0) = 1, P 2 , o(0) = 0, and the derivatives with respect to £ at 
£ = of p + = (Pi,o + P2,o — l)/2 and p~ = (Pi,o — Pa,o — l)/2 equal to and C, respectively. The 
longer proof below explains how the formulas for ^1^(0 and 7T2,o(£) could have been found. 

The system ( |5.52 ) is equivalent to the second order differential equation 

' (5.60) 



d 2 vr 



1 dlT l,o 

t dt 2 v 



1 + 



7T 



for 7r(t) 



dt 2 tdt 2 l ' 25 t 2 
7Ti(i), when ^(t) is given in terms of ir(t) by means of the formula 

dvr(t) 2 



7T 2 (i) 



dt 



+ 



The autonomous limit equation of (5.60) for t 
potential energy 1^(11) : = 



5t 

00 is n" 



7T(t). 



(5.61) 



(n 2 



l)/2, a Newton equation with 
n 3 /6 + II/2. It follows that the total energy £ = (lT) 2 /2 + V(n) is a 
constant of motion, and the solution which converges to the equilibrium point (II, II') = (1, 0) has 
energy E = V(l) = 1/3. This leads to the first order differential equation (IT) 2 = 2 (1/3 — V(II)) = 
(n-1) 2 (H" + 2)/3, hence, if the independent variable is denoted by s, ds/dU = 3 1/2 (II -l)" 1 (11 + 
2) 1 / 2 , when the substitution II + 2 = * 2 leads to 



2\/3y (^ 2 

*(s) + V3 



3) _i d^ + c = log 



noo : 



V3" 
*( S ) 2 



c e 



1 + 



c£ 44> *(s) 
12c£ 



^(s) + \/3 
c£ — 1 



+ c 



where c denotes a constant which at every place might be a different one. The function s 1— > II(s) 
is the "de gener ate elliptic function" of p. 38]. With the substitution £ = e _s , the differential 
equation (5.40) is equivalent to dPo/ds = v(0, Fq). As the derivative of Pq~(£) and Fq(£) at 
£ = have to be equal to and 1, respectively, the derivative of II = tt\ = 1 + p + + pT with 
respect to £ at £ = has to be equal to 1. Therefore c = 1 /12 a nd II = 1 + £/( £/12 — l) 2 , 
which proves ^^(C) = Pl,o(£)/(£ — 12) 2 with Pi,o(£) as in (5.56). The formula (5.61) with 
d/ dt = (-1 - 1/2*) £d/d£ yields 



7T2,0 



-C d7ri,o/d£ = P 2)0 (6/K-12) 



with P2,o(£) as in (5.57). 

Because 7Ti,o and 7^,0 are only singular at £ 



5.11 



12, it follows from Lemma 5.8 and Lemma 
that the functions ttu and 1^2,1 have a complex analytic continuation along any path in C \ {12} 
and that the asymptotic expansion (5.55) holds along these paths. If £ runs around 12 along a 



40 



small circle, then the i n (g) with large modulus, see Lemma 5.10, return to the same value. As the 
Painleve property implies that the function pc{t) is single valued, the asymptotic expansion ( 5.55| ) 
implies by induction on I that ii\ i and 1^2,1 are single valued complex analytic functions on C\{12}. 
For our system (5.52) the differential equations (5.41) for i £ Z>o take the form 

dg 



g 



d7Ti 



"g 



dg 

d7T 2 ,i 

dg 



7T 2 , 



7T1 



71% TTl.i 



. 1 , d7r 2; i_i 
+ 2^^^ 



+ i 



1 



7%j7% 



(5.62) 



(5.63) 



Given wi t j and 7r 2 ,j for j < 
differential equations for (7%i, 
equation 



i, (5.62), (5.63) is an inhomogeneous linear system of first order 
7r 2, j), equivalent to the inhomogenous linear second order differential 



d 2 71%, 



+ g- 



d7T 



l.t 



dg 

d7T 2 , 



dg 



d7Tl 



dg 



7Tl,0 7ri,i + ~ g ' + (l - -)7T 2 ,i-l + " 2^ 7%j 



-g 



dg 



dg 

1 d7i%j_i 
2^ dg 



TTl, 



+ i 



r) TTl, 



for 7Ti j, when 7r 2j i can be solved from (5.62) in terms of 7ri j and tti j_i 



Let y?i(g) and 922(g) be a basis of solutions of the homogeneous linear second order differential 
equation ir" = a(g) 7r'+6(g) ir. Lagrange's method of variations of constants yields that the solutions 
of the inhomogeneous equation ir" = a(g) ir' + 6(g) ir + /(g) are of the form 



tt(0 



ci ¥>i(£) + c 2 (^ 2 (g) + / 



"¥>l(0 Pzfa) + ^(g) VI 



/fa) 



dr/. 



w(rj) 

-aw, hence 



(5.64) 



Here w = ip\ (p' 2 — y>\ v? 2 is the Wronskian determinant, which satisfies w' 

The differential equation ir" = a(g) ir' + 6(g) 7r has a regular singular point at g = S if a(g) and 6(g) 
have a pole of order < 1 and < 2 at g = H. If a(g) = A (g - H) _1 + 0(1) and 6(g) = B (g - H)~ 2 + 
0((g — S) _1 ) as g — > 3, and the indicial equation A (A — 1) = ^4 A + -B has to distinct solutions Ai 
and A 2 , then there is a basis of solutions </>i(g) and </? 2 (g) such that </?fc(g) = (g — H) Aft (l + o(l)), and 
w (0 = (^2 — Ai) (g — H)* 1+ * 2-1 . See for instance Coddington and Levinson [H Chap. 4, Sec. 8]. 
Therefore, if 6(g) = O(g-H) 6 , then the integral in ( 5.64 ) is of order (g-~) b+2 as g -> E. If A x = A 2 , 
then one has to replace 992(g) = (g - H) A2 (1 + o(l)) by tp 2 (£) = ((g - H) Al log(g - H)) (1 + o(l)), 
with a corresponding change in the estimate for the integral in (5.64). 

" 2 (12) ; 



At g = 12 we have A = 
are Ai = —3 and A2 = 4. 



and B = (12) 
If, for every < j < i 



= 12, when the solutions of the indicial equation 
1, ttij and TV2,j are meromorphic at g = 12 with 
a pole of order < j + 2 and < j + 3, respectively, then an inspection of the inhomogenous terms in 
the second order diffrential equation for 7%j yields 6 = — i — 4, and because 6 + 2 = — i — 2 < —3 
it follows that m i has a pole of order < i + 2 at g = 12. For % = 1 we have 6 = —4, but then 
Ai = —3 yields that 71% 1 has a pole of order < 3. Subsequently (5.62) implies that 7^,1 has a pole 
of order < % + 3 at g = 12. It follows by induction on I that, at g = 12, mj and 7T2,; have a pole of 
order < I + 2 and < I + 3, respectively. 

At g = 00 we have ^4 = — 1 and B = 1, when the solutions of the indicial equation are Ai = 1 
and A2 = —1. If, for every < j < i — 1, 7ri i3 - and 7T2, j are meromorphic at g = 00 with exponents 



< j, then an inspection of the inhomogenous terms in the second order diffrential equation for 7Ti , 
yields b = i — 2, and because 6 + 2 = z>lit follows that at £ = oo the function 7Ti j has an 
exponent < i. Subsequently (5.62) implies that also 1^2,% has an exponent < i at £ = 00. It follows 
by induction on Z that the functions 7Ti ; and %il have exponents < I at £ = 00, and therefore 
are rational functions of the form tt m (£) = (f - 12)~'- 2 Pi,j(£) and 7r 2 ,i(0 = (f - 12)-'" 3 P 2 , K0> 
where Pi / and P2 z are polynomials of degree < 2 Z + 2 and < 2 Z + 3, respectively. 

The functions 7ri,i(£) = (£ - 12)~ 3 and 7r 2 ,i(Q = (£ - 12)" 4 P 2 , i(£), with th e respective 



polynomials Pi,i(£) and P 2 , i(£) as in (5.58) and (5.59), have been found with the help of a formula 
manipulation computer program, in the following way. The solutions of the system (5.62), (5.63) 
for i = 1 which are complex analytic in a neighborhood of £ = are of the form 

""1, 1 (0 = £ ((720 c- 84672) + (60 c + 4464) £ + 180 £ 2 -£ 3 )/(60(£-12) 3 ), 

7r 2i 1 (£) = (497664 + (8640 c - 1306368) £ + (2880 c - 124416) £ 2 + (60 c + 17568) £ 3 

-24£ 4 + £ 5 )/(60(£-12) 4 ), 

where c is free constant. With these functions tt\ 1 and 7T2 1, an investigation of the explicit solutions 
7Ti 2, ^2 °f the system (5.62), (5.63) for i = 2 yields that there exist solutions which are complex 
analytic in a neighborhood of £ = if and only if c = 678/5. Therefore 7Ti ) i(£) = (£ — 12) -3 Pi,i(£) 
and 7r 2 ,i(0 = (£- 12)- 4 P 2j i(0 with Pl,i(£) and P 2 ,i(0 as in (|5.58[) and (|5.59|), respectively □ 



As the system (5.52) is just a rescaled version of (|2.2l), passing to the complex projective plane 



and successively blowing up the base points of the vector fields, as in Section [2j leads to surface 
1S9, with a locus / where the vector field is infinite equal to the union of nine complex projective 
lines L^ 9 , < i < 8, and a pole line L9 \ /, where Lg is the complex projective line appearing 
at the last blowup. In the common domain of definition of the coordinate systems (vrjji, ir^) and 
(tti, W2), an application of the birational transformation (wi, 112) ^ [ffiji, ^132) to the asymptotic 
expansion (5.55) leads to an asymptotic expansion 



m— 1 



(5.65) 



1=0 



as t £ V, \t\ —> 00, where the functions itiik,l{ T ) are rational expressions in the functions ttij'{t), 
^2 i'( T ) for < V < Z, and therefore are rational functions of r. Because the differential equations 
for iTijk,o(T), analogous to (5.40), correspond to the autonomous limit system, these differential 



equations are regular and its solutions have a complex analytic extension as long as they remain in 
in the complement of the inifinity set I in the coordinate system (^1, 71^2 )• Because also the non- 
autonomous vector field is regular in S\I, the functions Ttijk,i(T) with Z € Z>o satisfy inhomogeneous 
linear differential equations as (5.41 ), variational equations of the differential equations for 7Tijk,o( T ), 
where the inhomogeneous term is a regular expression in the 7Tijk,m( T ) with m < I. It follows by 
induction on Z that all the rational functions 7Tijk,i{T) are regular, when the perturbation argument 
in the last paragraph of the proof of Lemma 5.11 yields that the asymptotic expansion (5.65) 
extends to the whole complement of / in the coordinate system (Vyi, ^2). 

The pole line is visible in the coordinate system (7rgn, ^912) as the line 7Tgi2 = 0, and therefore 
the poles are the solutions T of the equation 7rgi2(T) = 0, where 7Tgi2 = 7Ti/7r 2 . This leads to the 
following asymptotic resul ts for the poles, where in view of ( |5.54[ ) the poles of the corresponding 
truncated solution of (O) are given by X n = -2~ 3 / 5 3~ 1/5 (5T„/4) 4 / 5 . 



Lemma 5.18. There exist universal sequences of coefficients Cj, dk,i j, k, I G Z>o, with cq = 12, 
ci = 10 9/10, and do,o = 0, such that following holds. Let C / 0, and let pc{t) be the solution in 
Lemma 5.11, of the system (5.3) obtained from (5.52) by means of the substitutions %\ = l+p + +p~ , 



7T2 = p + — p . Let (7Ti(t), TT2(t)) be the corresponding solution of (5.52). Then there is a sequence 
of poles T n , n G Z, n >> ofiri(t), such that T n = 2ir\ n(l + o(l)) as n — > oo, and 

t(X„) := e~ T " T n ~ l/2 ~ - ^ CjT n -> as n -> oo. (5.66) 



Furthermore, with the notations 



1 log(2vrin) C u 

u := , v := , and W := u log , 

2vrin 2vrin ' & 12 2' 

where log(27rin) = log(27rn) + iri /2, we have the more explicit but more complicated asymptotic 
expansion 

1 C 

T n ^ 2nin - -log(2-K in) + log— + ^2 d k ju k W l as n oo, (5.67) 

fc,/>0 

of which the leading terms yield 

1 , / . \ , C v , 1 C 109. 
T n = 2vrm- - log(2vrin) +log— + - - (- log— + — )n 

1 1 C 139 

+ 16 V " ~ ( 4 log 12 + 240 } U W + ° (n ~ 2) (5 ' 68) 

as n — > oo. If C £ C \ {0} runs once around the origin in the positive direction, then T n moves 
continuously to T n+ \. 

Conversely, for every n > there exists an r > 0, such that these T n are the only poles T of 
7Ti(t) such that \T\ > r, ImT > 0, and \ e~ T T~ l l 2 \ < r\. 



Proof. It follows from Lemma 5.17 that the poles T of vri(t) with bounded H(T) := C e~ T T -1 / 2 
satisfy H(T) —> 12 as \T\ —> oo. 

There exists a sequence of rational functions 7rgi2 i, I G Z>o, such that, with the notation 
£ = Ct = C e~ l t- 1 ' 2 , 

/x m— 1 

7T 9 i2(t) = ^ = E ^^i2,Ke)+o(n, 

for every m G Z>o- Lemma |5. 17 implies that 

7r 9 i2,o(0 = (£- 12)(144 + 120£ + £ 2 )/(144£(£ + 12)) and 
7r 912 ,i(0 = (-71663616 - 40310784 £- 248832 £ 2 - 11860992 £ 3 - 1221696 £ 4 
+1224£ 5 - 240£ 6 - £ 7 )/((1244160£ 2 (£ + 12) 2 ). 

Because 71912,0(12) = and 7Tg 12 o(12) = 1/24 7^ 0, an application of the implicit function theorem 
yields that there exist open neighbourhoods A and B of 12 and 0, respectively, such that for every 
f 1 £ A the equation 

m— 1 

E t~ I T912,l(0 = 
1=0 

has a unique solution ^ m — 

(f 1 ) G B, which moreover depends in a complex analytic fashion on 
t , and satisfies S m (0) = 12. Furthermore, E(T) = H m (T _1 ) + 0(T _m ), and as the left hand side 
does not depend on m, it follows that the coefficients Cj for < j < m — 1 in the Taylor expansion 
of the function H m at the origin do not depend on m. Because this holds for every m G Z>o, it 

43 



follows that there is a sequence of complex numbers c,-, j E Z>o, such that the poles T of 7Ti(t) 
with bounded S(T) := C e~ T 7 1-1 / 2 satisfy 

Cr(T) = C e~ T T" 1 / 2 ~ 12 + c,- T _i as \T\ -»• oo. 

i>o 

Because 7Tgi2,i(12) = —109/240, we have ci = — 7rgi2, 1 (12)/7Tg 12 q(12) = 109/10. This completes 
the proof of ( 5.66|). 

The proof of (5.67) is analogous to the proof of ( |5.46[ ). For any m G Z>o, (5.66) yields that 

to— 1 



19 1 "" " 

r = r(T n ) = -(l + -^ c^T^+r), 



0(n" 



hence 



12 



TO— 1 



log ~ - log — + log(l + ^ ^| r n ^ + r). 



1 C 
2tt i n - - log(27r i n) + log — + S 

2vrin(l - -v + ulog(C/12) + uS) = 2vrin (1 + W + uS), 



log T n = log(2^ i n) + log(l + W + uS) 



where r = 0(T n n 
Upon the substitution 

T n = 

which implies 
and 

T n ~ l = u{l + W + uS)-\ 
the equation T n = 2nin— (1/2) logT n — logr is equivalent to the equation 

1 m— 1 

S = F(u, W, r, S) := — log(l + W + uS) - log(l + ^ ^ u j (1 + W + uS)~ j + r). 

j=o 

Because ^(0, 0, 0, S) = 0, it follows from the implicit function theorem in the complex analytic 
setting that there exist open neighborhoods U, W, R, and S of the origin in C such that for every 
(u, W, r) E U x W x R the equation has a unique solution S = S m (u, W, r) 6 S, and that S m is a 
complex analytic function on U x W x R. In our setting 

m— 1 m—k 

S m (u, W, r) = S m (u, W, 0) + 0(r) = J] £ 4,/ u fc ^ + 0(n" m ). 

fc=0 «=o 

Here the coefficients do not depend on m because S does not depend on m. This completes 
the proof of (5.67). 



The equation for m = 2 yields 

1 (W + u S) + - (W 2 + 2 W u S) - ^ u (1 - W) + 0(n~ 2 ), 



S 



hence 



2 v ' 4 v 7 12 



S= -\ W -^2 U+ \ w2 + { \ + T 2 )Wu + ^ 2 ^ 



which implies (5.68). 

If C E C \ {0} runs once around the origin in the positive direction, then (5.66) implies that 
r(T n ) runs once around the origin in the negative direction, when Lemma 5.10 implies that T n 
moves continuously to T n+ \. □ 
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Figure 5.5 illustrates the asymptotic approximations of the poles in (5.68) i n the complex t-plane, 

. t/6/5 71-. 



of Boutroux's triply truncated solution Pi(t), when, according to Remark 



Shown are the points in the right hand side of (5.68) without the remainder term 0(n ), for 



5.15 



C 



1 < n < 20. For clarity of the picture, the imaginary part has been multiplied by 1/24 in comparison 
to the real part. It would be interesting to compare the approximate poles in Figure [575] with the 
numerical approximations of the actual poles of Pi(t). 

Remark 5.19. The intriguing "General comments 2" of p. 40] say: 

"The pole structure can be explored beyond the first array, in much of the same way: For large £ 
induction shows that H n ~ Const £ n , suggesting a reexpansion for large £ in the form 



h 



oo 

E 

k=0 



-i 



(118) 



By the same technique it can be shown that (118) holds and, by matching with 



X>- fc #,(£(*)) 



fc=0 



at £2 ~ x 
i.e. 



-2/3 



we get = Hq with = -1/60. 
£(x s )=6 = 12 + ^ + 0(x- 2 ), 



10s 



the second line of poles is given by the condition 

-3/2 



-XI 



(113) 

Hence, if x s belongs to the first line of poles, 

(116) 

-60 • 12 c 



i.e., it is situated at a logarithmic distance from the first one: 

x\ — x s = — lnx s + (2 n + 1) 7r i — ln(60) + o(l). 

Similarly one finds x Sj 3 and in general x s>n . The second scale for the n-array is x - ™ -1 / 2 e~ x . 

The expansion (113) can however matched directly to an adiabatic invariant-like expansion valid 
throughout the sector where h has poles, similar to the one in Joshi and Kruskal [16J. In this 
language, the successive expansions of the form (118) pertain to the separatrix crossing region. We 
will not pursue this issue here." 

The word "suggesting" preceding (118) indicates that (118) is a conjecture, but in the sequel all 
the statements, including (118), are treated as facts, with only some hints of proofs. The phrase 
"matching with (113) at £2 ~ x- 2 / 3 " sug gests that H n ~ Const n £ n implies that the expansion 
(113) extends to domains where £(x) is of order x 1 / 3 , thus allowing a matching with (118) for 
£2(3;) = CW £(x) x" 1 of order x~ 2 / 3 . 

Write T^(i) := r(t)t~ N = e~* t~ N ~ 1 / 2 , the secon d sca le for the (N + l)-st array of poles. 
Because nk,l(Q ~ Constfc^ for £ — > 00, see Lemma 5.17, the t~ l 7Tk,i(C r(t), I G Z>o form an 
asymptotic sequence for \t\ — > 00, \r(t)\ = o(\t 



and r(t) bounded away from 12/ 'C, in the sense 
that t~ l 7Tkj(C r(t) = 0(r[ 1 l(t) / ), where r^(t) — > 0. A stronger conjecture would be that (5.55) 
extends as an asymptotic expansion for \t\ — > 00 in the aforementioned domain. Because (5.56) 
and (5.57) imply that (ir^ofjt;), tt2,o(£)) - > (1) 0) as £ — > 00, it would follow that (vri(t), TT2\t)) 
converges to the equilibrium point (1, 0) of the autonomous limit system if \t\ — > 00, \r(t)\ — > 00, 



and |r(t)| = o(|t|). In view of ( 5.56| ), (5.57), (5.58), and (5.59), the first two terms of the extended 
asymptotic expansion yield 

7ri(i) = 1 + 144 (C t (t))- 1 - C rW(i)/ '60 + o(r(t)- 1 )+o(rW(t)) ) 

7r 3 (i) = 144(C , r(t))- 1 + CrW(t)/60 + o(r(i)- 1 ) + o(TW(t)). 



-Li 



Therefore, if we restrict to r(t)" 1 = o(rW(i)) = o(r(t)/t), that is \T(t)\/\t\ 1 / 2 -> oo, then 

(7r 1 (t) ! 7r 2 (i)) = (l,0) + ^(l, -l) + o((eW)) if £ [1] (*) = -CrW(t)/60, 
the leading term of an asymptotic expansion 



As in Lemma 



5.8 



(7n(t), 7r 2 (t) ~ £ r' (vrW (eW(t)), vrW (fM(t))). 
rW Mih J 1 ! 



the function £M ^ J (f [1] ), 4 Q (C [1] )) satisfies the same differential equation 
(5.40) as the function £ i— )• (7Ti 5 o(£), 7T2,o(£))> where both functions have the same value a nd fi rst 
order derivative at the origin. Therefore 7r! = 7^0, with 7^0 and tt2,o as in Lemma 5.8 It 
follows, if the aforementioned statements about the asymptotic expansions hold, that the second 
sequence of poles occurs at points 2iri n — (3/2) log(2-7ri n) — log(— 720/C) + o(l), equal to the 
first sequence 2iri n — (1/2) log(27ri n) — log(12/C) + o(l) plus — log(27rira) — log(— 60) + o(l) as 
n — > 00. 

The text in "General comments 2" of [7J, p. 40] continues with the statement that for each N 
there is an asymptotic expansion of the form 

00 

(Mt), **(*)) ~ £ t- 1 (^}( C m T wi m v w ic w T m m 
1=0 

valid, if interpreted in the strong sense, for \t\ — > 00, |r [N] (t)\ = o(t), and \r [N - 1] (t)\ -> 00, where 
the constant depends linearly on C. This would lead to an asymptotic description of the 

(N + l)-st sequence of poles, equal to the iV-th sequence plus — log n + 77V + o(l) as n — > 00, where 
the constant 77V neither depends on n nor on C. For \t\ — > 00 and t between the A^th and the 
(N + l)-st sequence of poles, the solution (vri(t), ^{t)) converges to the equilibrium point (1, 0) of 
the autonomous system, if and only the distance from t to both sequences of poles tends to infinity. 
Furthermore, for every M > we have that the energy E = 7T2 2 /2 — 7ri 3 /6 + tt\/2 converges to 
the critical level 1/3, meaning that the solution converges to the solution of the autonomous limit 
system at the critical energy level, if Imi —> 00, Ret > — M log(Imt), and \t\ times the distance 
from t to the poles tends to infinity. The latter condition is related to the description of the energy 



near the poles in (3.3). 

We would like to prove statements like those in the second paragraph in the "General comments 
2" of [TJ p. 40] by means of the averaging method. This is not a trivial matter, as all the asymptotic 
expansions up till now are near one of the critical values of the energy function, where solutions of 
the averaged differential equation for the energy function are not uniquely determined. One might 
expect that the energy function acquires different limit values from the critical value 1/3 of the 
energy function, if t runs to infinity in the direction of e 1 s with it/ 2 < 9 < 3tt/2. On the other 
hand the truncated solution converges to the equilibrium point of the autonomous limit system 
(with energy equal to the critical value 1/3) if —it/ 2 < 9 < tt/2. 



Appendix A. Okamoto's Space 

We construct Okamoto's space of initial conditions |21j in the Boutroux rescaling. (See also |10j 
for the original Painleve equation ( jl. 1[ ) .) Recall that the notation (w^i, Uij2 

) will be used to denote 

the coordinates in the j-th chart of the z-th blowup and that in each coordinate chart, the Jacobian 
of the coordinate change from (ui,U2) to (1^1,1^2) will be denoted by 

_ duiji duij 2 _ duiji duij 2 
lJ du\ du2 du2 du\ 
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Up to and including the seventh blowup, the function z E is rational with ittj,- in the denominator, 



and we have added the formula for Ewij in each coordinate chart. 



A.l. Embedding into P . Recall the embedding of (u\, U2) into projective space given in Section 



1.2. We have the second affine chart in P 2 : 



^021 = 1*1~\ 

UQ22 = ui^ 1 u 2 , 

Ul = 1*Q21~\ 

U2 = UQ21~ 1*022, 

1*021 = 1*021 (-W022 + 2 (5 Z 



-1\ 



«022 

U>02 = -1*021 



-"021 1 (6 + U021 2 - ""021 M022 2 - (5 Z) 1 U 21 U022), 

3 



^02 = 3n 02 i 3 (n 22 - 2 (5z) _1 ), 
Ew 02 = 2 + u 2i 2 - 2~ x u 2i U022 2 , 
Ew§2 = ~(5z)~ 1 (12 + 2u 2i 2 - 3u 211i022 2 ) 



The line at infinity Lq corresponds to U021 = 0. In this chart there are no base points for the 
Painleve vector field or the anticanonical pencil. 
Third affine chart in P 2 : 



1*031 = U 2 ~ 1 , 

«032 = 1*1 U2~ , 

Ul = U031 _1 ^032, 

U2 = «031~ ) 



U031 = --U031 2 - 6n 03 2 2 + 3 (5z) 1 t^ 3 i, 

^032 = «031 _1 (^031 - «031 2 ^032 - 6«032 3 + (5 2)" 1 ^031 1*032) 



^03 = «031 3 , 



[u>03 1*032 3 ]* = 3ti 31 3 (-1 + 2 (5z) 1 n 32)«032 4 , 
E WQ3 = 2 _1 U 3l - M031 2 -"032 ~ 2 1X032 3 , 

Ew 03 = (Bzy 1 (-3iio3i + 2u 3i 2 1*032 + 12w 32 3 ) 



The line at infinity Lq corresponds to uo 3 i = 0. Both the The Painleve vector field and the 
anticanonical pencil both have a base point 60 given by 11031 = 0, U032 = 0. 



A. 2. Resolution of the flow at bo. Blowing up P 2 at 60 leads to Si. 
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First coordinate chart: 



^031 = WmUo32, 

= UU2, 

uiu = ui' 1 , 

u 112 = itiu 2 _1 , 

u\ = Mm -1 , 

u 2 = Uiii -1 «ii2 -1 , 

win = Uii2 _1 «m (-1 + 2(5z)~ 1 u 112 ), 

U112 = Mm -1 (u m - 6un2 2 - nm 2 uii2 2 + (5z) _1 niiiuii 2 ), 

w u = wiii 3 u.ii2 2 , 

[w n u 112 ~ 2 ]' = 3-uin 3 (-1 + 2 {hz)' 1 U112) un 2 -1 , 

[wn nin -1 ]* = 2um U112 2 (-6W112 - tim 2 u,ii2 + 3 (5z) _1 u m ), 

E w u = 2 _1 u m - 2 U112 2 - urn 2 u n2 2 , 

Ew u = (5Z)' 1 (-3um + 12un2 2 + 2u m 2 uii2 2 ). 



Then u\i 2 = defines L\ and u±u = defines Lq . The Painleve vector field and the anticanonical 
pencil both have a base point 61 given by u±u = 0, u\\ 2 = 0. 
The second coordinate chart after the first blowup is defined by 



■"031 




U121, 


UQ32 




U\22 «031, 


«121 




u 2 _1 = Uin U112, 


ni22 




«1 = Will -1 , 


Mi 




■"122, 


U2 




■"121~\ 


«121 




""121 (-1*121 - 6u m U122 2 + 3 (5 zy 1 ), 


«122 




U121' 1 (1 - 2 (5 z)' 1 U121 U122), 


1012 




■"121 2 , 


1012 




2ni2i 2 (-U121 - 6ni2i U122 2 + 3 (5z) _1 ), 


£'ttJi2 




2 _1 - «121 2 U U2 - 2 m 2 l 2 U122 3 , 


Ewi 2 




(5 z)" 1 (-3 + 2 u 121 2 ni22 + 12 u m 2 wi 22 3 ) 



The equation u\ 2 \ = defines L\. The line is not visible, and there are no base points in this 



chart. 
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A. 3. Resolution of the flow at b\. Blowing up S\ at b\ leads to 52- First coordinate chart: 



Will = 1*2111*112, 

1*112 = 1*212, 

1*211 = Ui~ 2 U 2 , 

n 2 i2 = u 2 ~ 1 u 1 , 

1*1 = U2U~ 1 U212~ 1 , 

1*2 = 1*211 _1 U 2 12~ 2 , 

1*211 = U212' 1 (-2^211 + 6u 2 i2 + n 2 n 2 u 2 i2 3 + (5^) _1 1*211 1*212), 

1*212 = 1*211 _1 (U211 - 6li212 - M211 2 1i212 3 + (5 z)~ X 1i 2 ll 1*212), 

H>21 = 1i211 3 U212 4 , 

[u>21 1i212~ 1 ]* = 3ll211 3 n212 2 (-l + 2(5z)~ 1 1l212), 

[1W21 ii2ii -1 ]* = 2u2iiU2i2 4 (-6-U2ii 2 n 2 i2 2 + 3(5z) _1 u 2 ii), 

E W21 = 2 _1 11211 - 211212 - U 2 11 2 11212 3 , 

Ew2\ = (5z) _1 (-3u 2 n + 12ii 2 i2 + 2u 2 ii 2 n 2 i2 3 )- 



Then 11212 = defines L2 and 11211 = defines the proper transform L Q ' of Lq . The proper 
transform of L\ is not visible in this chart. The Painleve vector field and the anticanonical 
pencil both have a base point 62 given by U211 = 0, 11212 = 0. 

The second coordinate chart after the second blowup is defined by 



lilll 




1*221, 


1*112 




1*222 1*111, 


1*221 




U-C 1 = 11211 11212, 


1*222 




1*1 2 U2 _1 = 1*211 _ \ 


U\ 




1*221~\ 


1*2 




-2 -1 
1*221 1*222 , 


U221 




1*222 _1 (-1 + 2 (5 zY 1 U 2 21 11222), 


11222 




1*221 _1 (2 — 611221 1*222 2 — 1*221 3 U222 2 — (5^) _1 11221 1*222), 


1*>22 




4 2 
1*221 H222 , 


li>22 




2 1*221 4 1*222 2 (-6 1*222 ~ 1*221 2 1*222 + 3 (5 z)' 1 ), 


-Elf 22 




2' 1 - 2 1i 2 21 1*222 2 - 1*221 3 **222 2 , 


EW22 




(5 z)- 1 (-3 + 12 U 2 21 1*222 2 + 2 1*221 3 «222 2 )- 



The equations 1*221 = and U222 = define L2 and respectively. The line is not visible, 
and there are no base points in this chart. 
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A. 4. Resolution of the flow at hi- Blowing up S2 at 62 leads to 53. First coordinate chart: 



1*211 = U311M212, 

^212 = «312, 

M311 = ^l _3- "2 2 , 

W.312 = flM2 _1 , 

"1 = ^311 _1 «312 _2 , 

^2 = ^311 _1 ^312~ 3 , 

«311 = ^312 _1 (12 - 3l/ 3 n + 2u 3 ii 2 n 3 i2 4 ), 

^312 = ^311 _1 (-6 + ti311 - n 3 n 2 ti3i2 4 + (5z) _1 U311 U312), 

W31 = «311 3, U312 6 , 

[w 3 i («sii - 4)" 1 ] * = 2 u 311 3 U312 6 (1x311 - 4)- 2 (-ix 3 n 2 x/ 3 i2 3 + 3(5 z)- 1 (U311 - 4)), 

Etu 3 i = -2 + 2 -1 U311 -n 3 n 2 n 3 i2 4 , 

Ewai = (5z)- 1 (12-3u 3 ii + 2u3ii 2 «3i2 4 ). 



Then 1/312 = defines L3 and 1/311 = defines the proper transform Lq of Lq . The proper 

transforms of L2 and of are not visible in this chart. The Painleve vector field and 
the anticanonical pencil both have a base point 63 given by 1/311 = 4, 1/312 = 0. 
The second coordinate chart after the third blowup is defined by 



«211 


= U321, 




■"212 


= U322U211, 




■"321 


= Ul~ 2 U 2 = U311 U312, 




■"322 


= Ui 3 U 2 ~ 2 = -U311~\ 




Ul 


-2 -1 
= U32I U322 , 




U2 


—3 —2 
= I/321 1/322 , 




U321 


= Xi322 _1 (-2 + 6 U322 + U321 4 M322 3 + (5 z) 


1 ^321 ^322), 


U322 


= «321 _1 (3 - 121/322 - 2l/321 4 U 3 22 3 ), 






6 4 
— «321 «322 , 






= 2 U321 6 1/322 4 (1 - 4 1i 3 22)~ 2 (--U321 3 «322 2 


+ 3 (5 ZY 1 (1 - 41/322)), 


Ew 3 2 


= 2 _1 - 2l/ 3 22 - n 32 i 4 ii322 3 , 




EW32 


= (5z)- 1 (-3 + 12ii322 + 2i/32i 4 n 3 22 3 ). 





The equations 1/321 = and 1/322 = define L3 and L\ , respectively. The lines Lq and L\ are 
not visible. The Painleve vector field and the anticanonical pencil both have a base point 63 given 
by 1/321 = 0, 1/322 = 1/4 in this chart. 
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A. 5. Resolution of the flow at 63. Blowing up S3 at 63 leads to S4. First coordinate chart: 
1*311 -4 = U411U312, 

1*312 = M412, 

W411 = ui~ A u 2 (-4i*i 3 + i* 2 2 ), 

U412 = 1*1 U2" 1 , 

Hi = -U412" 2 (4 + U 4 iiU4i 2 ) _1 , 

1*2 = U412~ 3 (4 + U411^412) _1 , 

1*411 = -U412 -1 (4 + U411^412) _1 

x (—10 1*411 - 4u 4 n 2 U412 + 128n 4 i2 3 + 112 ii 4 n u 4 i2 4 + 32n 4 n 2 ti 4 i 2 5 
+3n 4 ii 3 u 4 i2 6 - (5 U411U412 (4 + U411U412)), 
«4i2 = (4 + U411 U412) -1 

X (-2 + U411 U 412 - 16 U 4 12 4 - 8 U411 U412 5 - 1*411 2 1*412 6 
+ (5 1*412 (4 + 1*411 U412)), 

U>41 = U412 5 (4 + 1*411 W412) 3 , 

[i«4i i*4ii _1 ] = 2 U412 5 (4 + U411 U412) 3 i*4ii~ 2 (--U412 2 (4 + 1*411 W412) 2 + 3(5 z)~ l l*4n), 

Ewti = 2' 1 1*411 - 1*412 3 (4 + 1*411 1*412) 2 , 

EW41 = (5z) _1 (-3 1*411 + 2U412 3 (4 + 1*411 li412) 2 )- 

Then 1*412 = defines L4 and 4 + 1*4111*412 = defines the proper transform of The 
proper transforms of the other lines on which the Painleve vector field is infinite are not visible in 
this chart. The Painleve vector field and the anticanonical pencil both have a base point 64 given 
by 1*411 = 0, U412 = 0. 

The second coordinate chart after the fourth blowup is defined by 

11311-4 = 1*421, 

11312 = U422 (ll311-4), 

1*421 = Ul~ 3 (-4 Ui 3 +U 2 2 ) = 1*41111412, 

1*422 = 1*1 4 1*2 _1 (-4l*l 3 + 1*2 2 )~ 1 = 1*411 ~\ 

1*1 = l*42l" 2 (4 + U421)" 1 H422 -2 , 

1*2 = U421~ 3 (4 + U42l) _1 H422~ 3 , 

11421 = 11422" 1 (-3 + 32li42l 3 n 42 2 4 + 16l*421 4 l*422 4 + 2l*421 5 l*422 4 ), 

1*422 = 1*421 _1 (4 + U421) -1 (10 + 4l*42l - 128n421 3 U422 4 - 112u 4 2l 4 ll422 4 

-32 l*4 2 i 5 ll4 22 4 - 3li 42 i 6 l*42 2 4 + (5z) -1 li 42 l (4 + li 42 l ) 1*422) , 

WA2 = U421 5 (4 + 1*421 ) 3 H422 6 , 

li>42 = 2l*421 5 (4 + lt42l) 3 l*422 6 (-l*421 2 (4 + li42l) 2 li422 3 + 3 (5 2;)^ 1 ), 

Ew A 2 = 2 _1 - li 42 i 3 (4 + U 4 2l) 2 U422 4 , 

-Eli'42 = (5z) _1 (-3 + 2u 4 21 3 (4 + li 4 2l) 2 l*422 4 ). 

The equations 1*421 = 0, 4 + 1*421 = 0, and u 42 2 = define L4, and and respectively. The 
proper transforms of the other lines on which the Painleve vector field is infinite are not visible in 
this chart. Both the Painleve vector field and the anticanonical pencil have no base point in this 
chart. 
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A. 6. Resolution of the flow at 64. Blowing up £4 at 64 leads to £5. First coordinate chart: 



1*411 = 1*5111*412, 
U412 = 7*512, 

■"511 = U2 2 (~4ui 3 + t/ 2 2 ) Ui~ 5 , 

1*512 = 1*l1*2~\ 

Ml = M512" 2 (4 + U511 U 5 i2 2 )~\ 

1*2 = M512 _3 (4 + U511 U 5 12 2 ) -1 , 

1*511 = M512 _1 (4 + U511 U 5 12 2 ) -1 

x (-8 U511 + 128 U512 2 - 5 n 5 n 2 u 512 2 + 128 7/511 u 5 i 2 4 + 40 u 51 i 2 u 512 6 

+4 U511 3 U 5 i2 8 - 2 (5 U511 U512 (4 + U511 U.512 2 )), 

1*512 = (4 + W511 U 5 12 2 ) -1 

x (-2 + 7/511 u 5 i2 2 - 16 U512 4 - 8 7/511 U512 6 - U511 2 U512 8 
+(5 z)" 1 7/512 (4 + U511 u.512 2 )), 

H>51 = M512 4 (4 + n 5 iiU 5 i 2 2 ) 3 , 

[w 5 i i*5ii _1 ] * = 2 n 5 i2 4 (4 + U511 n 5 i2 2 ) 3 u 5n ~ 2 (-1*512 (4 + U511 i*5i2 2 ) 2 + 3 (5 z)' 1 u 511 ), 
Ew 51 = 2 _1 U511 - U512 2 (4 + 7/511 U512 2 ) 2 , 

£^7^51 = (5z) _1 (-3 7/511 + 27/512 2 (4 + 7/511 7/512 2 ) 2 )- 

Then 7/512 = defines L5 and 4 + 7/511 U512 2 = defines the proper transform Lq 5 ^ of The 
proper transforms of the other lines on which the Painleve vector field is infinite are not visible in 
this chart. The Painleve vector field and the anticanonical pencil both have a base point 65 given 
by 1x511 = 0, 7i 5 i2 = 0. 

The second coordinate chart after the fifth blowup is defined by 

7/411 = 1*521, 

1*412 = 1*522 1*411, 

1*521 = 1*1~ 4 (-4 Til 3 + 7t 2 2 ) 7/2 = 1*511 7/512, 
1*522 = 1*1 5 (-47/1 3 + U 2 2 y 1 7*2~ 2 = 1*511 _1 , 
1*1 = 1*521~ 2 1*522~ 2 (4 + 7/521 2 7i 5 22) _1 , 
1*2 = 1*521 _3 1*522~ 3 (4 + 7/521 2 7t522) _1 , 
1*521 = 1*522 _1 (4 + 7Z521 2 7i 5 22)~ 1 

X (-10 - 4 7/521 2 7/522 + 128 7i 5 21 2 1*522 3 + 112 7/ 5 21 4 1*522 4 + 32 7i 52 l 6 7/ 5 2 2 5 
+3t/ 5 21 8 1*522 6 - (5^) _1 1*521 7/ 5 22 (4 + 1* 52 1 2 7/ 5 2 2 )) , 
1*522 = 1*521 _1 (4 + 7/521 2 T*522)~ 1 

X (8 + 5 7i 52 l 2 7/522 ~ 128 Tf 5 21 2 7/ 52 2 3 ~ 128 7Z 52 1 4 7/ 52 2 4 ~ 40 7t 5 21 6 1* 52 2 5 
-4//521 8 1*522 6 + 2 (5 zY 1 7/521 7/ 52 2 (4 + U^ 2 2 1*522)), 
1i>52 = 1*521 4 1*522 5 (4 + 7/ 5 21 2 1*522) 3 , 

77)52 = 2 7/521 4 7/522 5 (4 + 7t 52 l 2 T*522) 3 (-1*521 1*522 2 (4 + U^ 2 7/ 52 2) 2 + 3(5z) _1 ), 

E 77)52 = 2" 1 - 7/521 2 7i522 3 (4 + U 5 21 2 1/522) 2 , 

Ew 52 = (5z)- 1 (-3 + 2 7i521 2 7/ 5 22 3 (4 + 7i521 2 1*522) 2 )- 
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The equations U521 = 0, 4 + 11521 2 1*522 = 0, and U522 = define L5, L , and and L4 , respectively. 
The proper transforms of the other lines on which the Painleve vector field is infinite are not visible 
in this chart. Both the Painleve vector field and the anticanonical pencil have no base point in this 
chart. 



A. 7. Resolution of the flow at 65. Blowing up £5 at 65 leads to Sq. First coordinate chart: 



U511 = ^611 "512, 

1*512 = ^612, 

U6H = uf 6 U 2 3 (-4ui 3 +« 2 2 ), 

U 612 = U1U2' 1 , 

Ul = U 612 ~ 2 (4 + U6H W612 3 )" 1 , 

U2 = U 612 ~ 3 (4 + U 6 iiU 6 i2 3 ) _1 , 

■"611 = -U612" 1 (4 + 1X611 M612 3 )" 1 

x (-6 it 6 n + 128 u 612 - 6 u 6 n 2 u 6 i2 3 
+144 u 6 n -U612 4 + 48 i* 6 ii 2 n 6 i2 7 + 5 i* 6 n 3 u 6 i 2 
-3 (5 z)' 1 u 6 n i* 6 i2 (4 + n 6 n t*6i2 3 )), 

1*612 = (4 + 1*611 ^612 3 ) _1 

X (-2 + 1*611 1i 6 12 3 - 16 1*612 4 - 8 U 6 11 1* 6 12 7 
-1*6H 2 n 6 12 10 + (5 z)' 1 1*612 (4 + U 6 11 H612 3 ))) 
W 6 1 = 1*612 3 (4 + 11611 U612 3 ) 3 , 
[wqiUqh^ 1 ]' = 2u 6 12 3 (4 + U 6 ii1l612 3 ) 3 1i611~ 2 

X (-(4 + 1*611 U612 3 ) 2 + 3 (5 ZY 1 1*61l), 
Ewqi = 2 _1 1*611 -1*612 (4 + U611U612 3 ) 2 , 
Ewqi = (5z) _1 (-31*611 + 2u 6 i2 (4 + 1*611 i*6i2 3 ) 2 )- 



o (b)(5) 

Then 1*512 = defines Lq and 4 + uqu uqi 2 = defines the proper transform L of L . The 
proper transforms of the other lines on which the Painleve vector field is infinite are not visible in 
this chart. The Painleve vector field and the anticanonical pencil both have a base point b$ given 
by 1*611 = 0, U612 = 0. 
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The second coordinate chart after the sixth blowup is defined by 



■"511 = ^621, 
^512 = «622«511, 

UQ21 = U2 2 (-4 U! 3 + U 2 2 ) = U 6 H U612, 

■"622 = ni 6 u 2 " 3 (-4ui 3 + u 2 2 y 1 = Uqh^ 1 , 

Ul = UQ21~ 2 U622~ 2 (4 + U 6 21 3 U 6 22 2 ) _1 , 
U2 = n 6 21~ 3 «622~ 3 (4 + UQ21 3 UQ22 2 y 1 , 
UQ21 = «622 _1 (4 + ti 6 21 3 «622 2 ) _1 

X (-8 + 128 n 62 i -U622 2 - 5 U 6 21 3 U622 2 
+ 128u 62 l 4 U 6 22 4 + 40u 62 l 7 n 6 22 6 + 4«621 10 ^622 
-2 (5 zY X U 6 21 ^622 (4 + U 6 21 3 «622 2 )), 
■"622 = ^621 _1 (4 + 1i621 3 ^622 2 )~ 1 

X (6 - 128 U 6 21 ^622 2 + 6 U 6 21 3 «622 2 
-144u 6 21 4 -U622 4 - 48u 62 l 7 n 6 22 6 - 5li 6 21 10 ^622 
+3 (5 zY X U 6 21 ^622 (4 + U 6 21 3 «622 2 )), 
U> 6 2 = ^621 3 W 6 22 4 (4 + ti 6 21 3 «622 2 ) 3 , 
W622 = 2u 6 2l 3 ti 6 22 4 (4 + U 6 21 3 «622 2 ) 3 

x (-u 6 22 (4 + n 62 i 3 n 622 2 ) 2 + 3 (5 z)' 1 ), 

EWQ2 = 2" 1 - U 62 l n 62 2 2 (4 + W 6 21 3 1i622 2 ) 2 , 

Ew 62 = (5 2) _1 (-3 + 2n 6 2iu 6 22 2 (4 + n 6 2i 3 u 62 2 2 ) 2 ). 



The equations uq21 = 0, 4 + U621 3 ^622 2 = 0, and ^622 = define Lq, and and respectively. 
The proper transforms of the other lines on which the Painleve vector field is infinite are not visible 
in this chart. Both the Painleve vector field and the anticanonical pencil have no base point in this 
chart. 
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A. 8. Resolution of the flow at be- Blowing up Sq at b$ leads to S 7 . First coordinate chart: 



■"611 = «711«612, 
^612 = «712, 

""711 = ni _7 u 2 4 (-4ui 3 + u 2 2 ), 

"U712 = UlU2 _1 , 

ui = u 712 ~ 2 (4 + uni ^7i2 4 )~ 1 , 

«2 = ^712~ 3 (4 + U711 U712 4 )" 1 , 
""711 = ^712 _1 (4 + U 7U M712 4 )" 1 

x (128 - 4 um + 160 U711 U712 4 - 7 u 7 n 2 U712 4 
+56 u 7 u 2 u n 2 S + 6 U711 3 u 7 i2 12 
-4 (5 u 7 ii (4 + u n i ^7i2 4 )), 

""712 = (4 + U 7U U 7 12 4 ) _1 

X (-2 - I6M712 4 + U711 U 7 i2 4 - 8U7U U712 8 
-u 711 2 u n2 12 + (5 zy 1 u 712 (4 + U711 U712 4 )), 

W71 = U 7 12 2 (4 + U7HU712 4 ) 3 , 
Kl («711 - 32)" 1 ] * = -2 U712 2 (4 + U711 U712 4 ) 2 {U711 - 32)~ 2 

x (U712 3 (1024 - 64 u 7U + 512 u 7 n u 7 i 2 4 + 12 u 7 n 2 u 7 i 2 4 

+64li 7 ii 2 -U712 8 + U 7 n 3 U712 8 ) 

+(5 z)- 1 (32 - 3 «7ii) (4 + U711 «7i2 4 )), 
Ew 71 = 2~ 1 u 711 - (4 + n 7 nu 7 i2 4 ) 2 , 
Ew 71 = (5z)- 1 (-3u 7 ii + 2(4 + U7im7i2 4 ) 2 ). 



Then U712 = defines L 7 and 4 + u 7 u U712 4 = defines the proper transform Lq of Lq . The 
proper transforms of the other lines on which the Painleve vector field is infinite are not visible in 
this chart. The Painleve vector field and the anticanonical pencil both have a base point b 7 given 
by u 7 n = 32, U712 = 0. Remarkably this base point in the Boutroux coordinates does not depend 
on the independent variable z, whereas the seventh base point in the unsealed system is given by 
7/711 = 32 x, y n2 = 0. 

The second coordinate chart after the seventh blowup is defined by 



«611 = «721, 

^612 = U722U6U, 

U 7 21 = U!~ 6 U 2 3 (-4lii 3 + U 2 2 ) = U 711 U 712 , 

U 7 22 = U 1 7 U 2 ^ 4 (-4 til 3 + -U2 2 )" 1 = -"711~\ 

Ul = U 7 21~ 2 U 7 22~ 2 (4 + U721 4 -U722 3 )~\ 

U2 = U 7 21~ 3 U 72 2~ 3 (4 + U72i 4 U722 3 ) _1 
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""721 = U722 1 (4 + U 72 1 4 U 7 22 3 ) 1 

X (-6 + 128 U 7 22 - 6 U 7 21 A U 7 22 3 + 144 n 72 i 4 U 7 22 4 + 48 u 72 i S u 72 2 7 
+5 U 72 l 12 U 72 2 10 - 3 ( 5 zY 1 U 72 l U 72 2 (4 + U 7 2\ 4 U 72 2 3 )), 

U 7 22 = -U121' 1 (4 + U 7 2l 4 U722 3 )" 1 

x (4 - 128 u 7 22 + 7 M 72 i 4 U722 3 - 160 U721 4 u 72 2 4 - 56 u 72 i 8 «722 7 

-6U 72 1 12 U 7 22 W + 4:(5 zy 1 U 72 lU 7 22 (4 + M 7 21 4 ^722 3 )) , 
■^72 = «721 2 «722 3 (4 + ti721 4 «722 3 ) 3 , 

[w 72 (32 u 722 -I) -1 ]" = 2n 72 i 2 u 7 22 3 (4 + n 7 2i 4 n 72 2 3 ) 2 (32u 7 22 - l)" 2 

X (u 72 i 3 U 722 3 (-64 + 1024 U 722 + 12 M 72 l 4 t* 7 22 3 
+512 U 72 l 4 M 7 22 4 + ""721 8 ^722 6 + 64 U 72 l 8 U 7 22 ? ) 

+(5 z)- 1 (-3 + 32 u 722 ) (4 + u 72 i 4 u 72 2 3 )), 

Ew 72 = 2 _1 - U 72 2 (4 + U 7 21 4 «722 3 ) 2 , 

Ew 72 = (5z)- 1 (-3 + 2n 7 22 (4 + u 7 2i 4 u 722 3 ) 2 ). 

The equations ii 7 2i = 0, A + u 7 2i 4 u 7 22 3 = 0, and -u 7 22 = define L 7 , Lq \ and and respectively. 
The proper transforms of the other lines on which the Painleve vector field is infinite are not visible 
in this chart. The Painleve vector field and the anticanonical pencil have the base point M721 = 0, 
1*722 = 1/32 in this chart, which is equal to the previously found base point b 7 . 



A. 9. Resolution of the flow at b 7 . Blowing up S 7 at b 7 leads to S 8 . First coordinate chart: 
u 7n -32 = U8ii«7i2, 

^712 = U812, 

U8U = -ui~ s u 2 (32 ni 7 + 4ui 3 u 2 4 - U2 6 ), 

U 812 = «1-U2 _1 , 

Ul = ^812 _2 (4 + 32u 8 i2 4 + M811^812 5 ) _1 , 

u 2 = n 8 i2~ 3 (4 + 32u 8 i2 4 + u 8 iiu 8 i2 5 ) -1 , 

"811 = ^812~ 1 (4 + 32u 8 i2 4 + li 8 llU 8 12 5 r 1 

x [-2 (u 811 + 1024 u 8 i2 3 + 152 u 8 n ^8i2 4 + 4 u 8 n 2 u 8 i 2 5 ) 
+usi2 7 (32 + u 8 ii U812) (1792 + 64 n 8 n u 8 i 2 + 6144 u 8 i2 4 + 416 u 8 n u 8 i2 5 
+7u 8 n 2 n 8 i2 6 ) - (5 z)- 1 (128 + 5« 8 n 1*812) (4 + 32 w 8 i 2 4 + u 8 n u 8 i 2 5 )], 
^812 = -(4 + 32ii 8 i2 4 + itsii w 8 i2 5 ) -1 

x [2 - 16 n 8 i2 4 - u 811 u 8 i2 5 + 256 u 8 i 2 8 + 8 u 8 n u 8 i 2 9 + 1024 u 8 i 2 12 
+64u 8 n n 8 i2 13 + n 8 n 2 u 8 i2 14 - (5z) _1 u 8 i2 (4 + u 8 u u 8 i2 5 + 32u 8 i 2 9 )], 
W 8 l = U 8 12 (4 + 32 u 8 i2 4 + u 8U ii 8 i2 5 ) 3 , 

Ew 8 \ = 2 _1 u 8 n - u 8 i2 3 (32 + ii 8 iiu, 8 i 2 ) (8 + 32n 8 i2 4 + n 8 n n 8 i2 5 ), 

Ew 8 \ = (5z)~ l usu' 1 (-64 - 3 itsii «8i2 

+2u 8 i2 4 (32 + u 8U u 812 ) (8 + 32u 8 i2 4 + usu -U812 5 ))- 
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Furthermore 



[w 81 (u 8 u + 256 (5 z)- 1 )- 1 ] * = 2 u 812 (4 + 32 u 812 4 + u 811 u 812 5 ) 2 (u 8U + 256 (5 z)' 1 )- 2 

x [-n 8 i2 2 (-2 10 - 2 6 u 8 n u 8 i2 + 2 12 • 7 u 8 i 2 4 + 2 8 • 5 u 8 n u 8 i 2 5 + 2 2 • 3 u 8 n 2 u 8 i 2 6 

+2 15 • 3 n 8 i2 8 + 2 10 • 7 u 8 n n 8 i2 9 + 2 5 • 5 u 8 n 2 W + u 8U 3 u 812 u ) 

+(5 z)- 1 (2 2 • 3 «8ii - 2 12 • 5 n 812 3 - 2 5 • 3 3 « 8n n 812 4 + 3 u 811 2 u 812 5 + 2 15 • 5 u 812 7 

+2 10 • 5 « 8 ii u 8 i2 8 + 2 17 • 5 n 8 i 2 n + 2 13 • 5 u 8U u 812 12 + 2 7 • 5 u 811 2 u 812 13 ) 

+768 (5 z)~ 2 (4 + 32 u 812 4 + u 8 n u 8 i 2 5 )]. 



The equation u 8 \ 2 = defines L 8 and 4 + 32n 8 i 2 4 + u 8 n u 8 i 2 5 = defines the proper transform L 

(7) 

of L Q . The proper transforms of the other lines on which the Painleve vector field is infinite are 
not visible in this chart. The Painleve vector field has a base point b 8 given by u 8 u = —256 (5 z)^ 1 , 
u 8 i 2 = 0. In the Boutroux coordinates, this is the first base point which depends on the independent 
variable z. The anticanonical pencil has a base point £>| u given by u 8 n =, ii 8 i 2 = 0. We have 
6 8 n b 8 with a distance between both base point vanishing of order 1/z as z — > oo. 
The second coordinate chart after the eighth blowup is defined by 



u 7U - 32 = u 82 i, 

^712 = u 822 (u ni - 32), 

u 82 i = -ui' 7 (32 ui 7 + 4ui 3 u 2 4 - u 2 6 ) = u 8U u 8 i 2 , 
u 822 = — ui 8 u 2 ~ l (32 ui 7 + 4ni 3 u 2 4 — m 2 6 ) -1 = u 8 ii~ L , 
ui = M82i~ 2 ^822~ 2 (4 + 32u 82 i 4 n 822 4 + u 82 i 5 u 822 4 )~ 1 , 
u 2 = M 82 i _3 n 822 _3 (4 + 32u 82 i 4 n 822 4 + u 82 i 5 u 822 4 )" 1 , 
u 82 i = m 8 22 _1 (4 + 32n 82 i 4 u 822 4 + ti 82 i 5 n 822 4 ) -1 

x [-4 - 2048 n 82 i 3 u 822 4 - 288 u 82 i 4 u 822 4 - 7 u 821 5 u 822 4 

+2 u 82 i 7 (32 + u 82 i) 2 u 822 8 (28 + 96 u 82 i 4 u 822 4 + 3 u 82 i 5 u 822 4 ) 

-4 (5 z)~ x (32 + n 82 i) u 822 (4 + 32 u 82 i 4 u 822 4 + ii 82 i 5 u 822 4 )], 

^822 = -1*821 _1 (4 + 32u821 4 W822 4 + ^821 5 ^822 4 ) _1 

x [-2 - 8 «82i 3 (256 + 38 u 821 + u 82 i 2 ) u 822 4 

+u 8 2i 7 (32 + n 8 2i) u 822 8 (1792 + 64 u 8 2i + 6144 u 8 2i 4 u 822 A + 416u 8 2i 5 ^822 4 

+ 7 «821 6 U 822 4 ) - (5 zY 1 (128 + 5 U 8 2l) U 822 (4 + 32 U 8 21 4 U822 4 + U 82 \ b U822 4 )], 
^82 = M82lW822 2 (4 + 32n 8 21 4 1i822 4 + 1i821 5 «822 4 ) 3 , 
Ew 82 = 2 _1 - U.821 3 (32 + U 8 2l) «822 4 (8 + 32U821 4 U822 4 + «821 5 M822 4 ), 

Ew 82 = (5 z)" 1 u 82 \~ l (-64 - 3«82i 

+2 u 82 i 4 (32 + u 82 i) u 822 4 (8 + 32 u 82i 4 u 822 4 + n 82 i 5 u 822 4 )). 
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Furthermore, 



[w 8 2 ^(s^^ + i)- 1 ]* = 

2 u S22 2 (4 + 32 u 821 4 u 822 4 + u 82 i 5 u 822 4 ) 2 (256 (5 z)~ l u S2 2 + 1)~ 2 

x [-U821 3 u S 22 3 (-2 10 - 2 6 u 8 2i + 2 17 ■ 7u 82 i 4 u 822 4 + 2 8 • 5 u 82 i 5 u 822 4 + 2 2 • 3 u 82 i 6 u 822 4 

+2 15 • 3 n 82 i 8 u 822 8 + 2 10 • 7u 821 9 u 822 8 + 2 5 • 5 u 82 i 10 u 822 8 + u 82 i n u 822 8 ) 

-(5 z)- 1 (-2 s - 2 2 • 3 u 821 - 2 18 m 821 3 u 822 4 - 2 11 • 3 2 u 821 4 u 822 4 

-2 5 • 5 u 82 i 5 u 822 4 - 3 u 82 i 6 u 822 4 + 2 20 • 7 n 82 i 7 u 822 8 

+2 16 • 5 n 82 i 8 u 822 8 + 2 10 • 3u 82 i 9 u 822 8 + 2 23 • 3 n 82 i U u 822 12 

+2 18 • 7n 82 i 12 u 822 12 + 2 13 • 5 u 821 13 u 822 12 + 2 8 u 821 14 u 822 12 ) 

+2 7 (5 z)~ 2 (2 7 + 11 u 821 ) (4 + 32 u 82 i 4 u 822 4 + u 82 i 5 u 822 4 ) u 822 ]. 



The equations -u 82 i = 0, 4 + 32 u 82 i 4 u 822 4 + ii 82 i 5 u 822 4 = 0, and « 822 = define L 8 , L Q , and and 

L^p , respectively. The proper transforms of the other lines on which the Painleve vector field is 
infinite are not visible in this chart. The Painleve vector field has the base point 6 8 defined by the 
equations u 82 i = 0, 256 (5 z)~ l -u 822 + 1 = 0. The base point 6 8 n of the anticanonical pencil is not 
visible in this chart. 



A. 10. Resolution of the flow at 6 8 . Blowing up 5 8 at 6 8 = b%{z) leads to Sg = Sg(z). First 
coordinate chart: 



U 8 n + 256 (5 Z) 1 = U911«812, 

""812 = «912i 

UQii = Ul~ 9 U 2 {-32ui 7 U 2 - Au 1 3 U2 5 +U 2 7 + 256 (5 z)' 1 Ui 8 ), 

Ug 12 = U 1 U 2 ~ 1 , 

ui = n 9 i 2 -2 (4 + 32u 9 i 2 4 + n 9 n ii 9 i 2 6 - 256 (5z) -1 «9i2 5 ) -1 , 
u 2 = n 9 i 2 - 3 (4 + 32u 9 i 2 4 + n 9 iin 9 i 2 6 - 256 (5z)- 1 u 912 5 )- 1 , 
mil = (4 + 32n 9 i 2 4 + ugii u 9 i 2 6 - 256 (5 z)~ l u 9 i 2 5 ) -1 

x [u 9 i 2 (-2 11 - 2 6 • 5«9ii u 912 2 + 2 13 • 7u 9 i2 4 - 3 2 u 9 n 2 u 912 4 

+2 12 ugn u 912 6 + 2 16 • 3 u 9 i2 8 + 2 3 • 3 2 u 9 n 2 u 9 i 2 8 

+2 12 • 5 n 9 ii u 912 w + 2 6 • 11 u 9 n 2 n 9 i 2 12 + 2 3 u 9 ii 3 n 9 i2 14 ) 

-2 (5 z)- 1 (2 2 • 3 ugn - 2 12 • 3 2 u 9 i 2 2 - 2 5 • 3 2 • 7 u 911 u 912 4 

+2 15 -3-5 u 912 6 + 3 u 9 n 2 u 9 i2 6 + 2 10 • 17 u 9U u 912 8 

+2 17 • 19 M 9 i2 10 + 2 13 • 3 • 7 ugn u 9 i2 12 + 2 7 • 23 u 9 n 2 m 9 i 2 14 ) 

+2 9 (5 z)- 2 u 912 3 (-2 6 -3-5 + 3 u 9U u 912 2 

+2 13 n 9 i2 4 + 2 14 • 5 u 912 8 + 2 8 • 11 U9n n 9 i 2 10 ) 

-2 24 -7(5z)- 3 n 912 12 ] 
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n 9 i2 = -(4 + 32n 9 i2 4 + U9iiu 9 i2 6 - 256(5z)- 1 u 912 5 )" 1 

x [2 - 2 4 u 912 4 - ugii ^9i2 6 + 2 8 u 912 S + 2 3 U9ii u 912 10 
+2 10 u 912 12 + 2 6 ugii u 912 u + u 911 2 u 912 w 
-(5z)~ l u 912 (2 2 - 2 5 • 7u 912 4 + n 9 n u 9 i 2 6 
+2 11 u 9l2 8 + 2 14 u 9 i2 12 + 2 9 ugii n 9 i 2 14 ) 
+2 8 (5z)- 2 U9 i2 6 (l + 2 8 n 912 8 )], 
w 91 = (4 + 32u 9 i 2 4 + U9iiU9 12 6 - 256 (5z)- 1 n 9 i 2 5 ) 3 , 
WQi = 3u9i2 3 (4 + 32u9i2 4 + u 9 iin 9 i2 6 - 256 (5z)- 1 U9i2 5 ) 2 ) 

x [-2 6 - 3 «9n u 912 2 + 2 9 u 9 i2 4 + 2 4 u 9 n u 912 6 + 2 11 u 912 8 
+2 7 ngn u 9 i 2 10 + 2 ugn 2 n 9i2 12 

-2 8 (5 z)~ l u 9l2 (-3 + 2 4 u 912 4 + 2 7 n 912 8 + 2 2 u 9n u 912 10 ) 

+2 17 (5z)- 2 n 912 10 ], 
Ew 9 i = -2" 1 ti 9 i 2 _1 [ugi2 (-1*911 + 2 9 n 9 i 2 2 + 2 4 u 9 n n 9 i 2 4 

+2 11 u 912 6 + 2 7 u 9U u 912 8 + 2 Ugii 2 u 912 10 ) 

-2 8 (5 z)- 1 (-1 + 2 4 u 912 4 + 2 7 ngi 2 8 + 4 u 9U W°) 

+2 17 (5zr 2 u 912 9 ], 
Ew 9l = (5z)" 1 n9i 2 " 2 [-2 6 -3M 9 iiiigi 2 2 + 2 9 n 9 i 2 4 + 2 4 n 9 iingi 2 6 

+2 11 u 912 8 + 2 7 ugn u 912 10 + 2 u 911 2 u 912 12 

-2 8 (5 z)- 1 u 912 (-3 + 2 4 u 9 i2 4 + 2 7 n 9 i 2 8 + 2 2 u 911 u 912 10 ) 

2 17 (5z)- 2 i i912 10 ]. 



As the change of coordinates from (ui, u 2 ) to all previous coordinate systems (it^i, Uij 2 ) for i < 8 
do not depend on z the limiting system of differential equations °ui = u 2 , °u 2 = 6u\ 2 + 1 in the 
coordinate systems (lijji, iijj 2 ) for i < 8 is obtained by deleting the term in and itjj 2 which 
have a factor l/z in front. This is no longer true in the coordinate system (1x911, ix 9 i 2 ). However, 
the difference between the Painleve-Boutroux system and the limiting system still has a relatively 
simple expression: 



ugn -° lign = -2(5z) 1 1x912 2 (64 - 640 (5 z) 1 u 912 + 3 u 9U u 9 i 2 2 

U912 -° 1*912 = (5z) _1 1i 9 i2. 



The equation ugi 2 = defines L 9 and 4 + 32 ugi 2 4 + 1*911 ugi 2 6 — 256 (5 z) 1 u 9 i 2 5 = defines the 
proper transform of . The proper transforms of the other lines on which the Painleve vector 
field is infinite are not visible in this chart. The Painleve vector field is regular along L 9 , nonzero, 
and transversal to it. Moreover, the Painlevee vector field has no base points in this chart. On the 
other hand the blowing up of 5s in the point 6g, which is not the base point of the anticanonical 
pencil, causes Ew 9 \ to be infinite along L 9 , the line determined by the equation -u 9 i 2 = 0. The 
image Wftg 11 of 6| n in S 9 is not visible in this coordinate chart. 
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The second coordinate chart after the ninth blowup is defined by 



U 92 \ = U8U + 2 (5z)~ , 

^812 = u 922 (u 8 n + 2 8 (5z) _1 ), 

W921 = ur 8 (-2 5 Ul 7 U2-2 2 ui 3 u 2 5 + u 2 7 + 2 8 (5zy 1 u 1 8 ) = u 9 i 1 u 912 , 
u 92 2 = m 9 U2- 1 (-2 5 Ul 7 u 2 - 2 2 Ul 3 u 2 5 + u 2 7 + 2 8 (5 z)- 1 ii! 8 )- 1 = u 9U -\ 

Ul = U 9 21 -2 U 922 T 2 (2 2 + 2 5 U 9 21 4 U922 4 + «921 6 ^922 5 ~ 2 8 (5 z)^ 1 U 92 i 5 7ig22 5 )~\ 
U 2 = U 921 - 3 U 922 ~ 3 (2 2 + 2 5 U 9 2l 4 U 9 22 4 + n 92 l 6 U 9 2 2 5 - 2 8 (5 Z)" 1 M 9 2l 5 Ug22 5 )~\ 
«921 = -n922" 1 (2 2 + 2 5 U921 4 U922 4 + ^921 6 n922 5 - 2 8 (5z)" 1 U921 5 n 9 22 5 ) _1 

x [2 + 2 11 n 92 i 2 n 9 22 3 + 2 4 • 19 u 9 2i 4 u 922 4 + 2 3 n 92 i 6 u 922 5 

-2 13 • 7U 921 6 U 922 7 - 2 8 • 3 • 5 U 9 2l 8 n 922 8 - 2 6 U 9 2l 10 M 922 9 

-2 16 • 3u92i 10 U922 n - 2 10 • 19n 92 i 12 n 922 12 - 2 7 ■ 5« 92 i 14 n 922 13 
-7n 92 i 16 n 92 2 14 + (Sz)' 1 u 921 u 922 

X (2 2 • 5 - 2 13 3 2 n 92 l 2 W " 2 5 • 7 • 17u 9 21 4 U 9 22 4 

+5 n 9 2i 6 n 92 2 5 + 2 16 • 3 • 5 u 92 i 6 u 922 7 + 2 15 U921 8 U922 8 

+2 18 • 19u 921 10 u 922 u + 2 16 • 5u 9 2i 12 u 922 12 + 2 s • 3 • 7u 92 i 14 u 922 13 ) 

-2 8 (5 z)- 2 u 921 4 u 922 5 (-2 7 -3-5 + 5 u 92 i 2 ^922 

+2 14 u 9 2i 4 U922 4 + 2 15 • 5 n 9 2i 8 U922 8 + 2 8 • 3 • 7u 92l 10 n 922 9 ) 

2 24 -7(5z)- 3 n92i 13 n 9 22 14 ], 

t*922 = ~(2 2 + 2 5 U92l 4 U922 4 + n 9 21 6 M922 5 - 2 8 (5z)" 1 n 9 21 5 U922 5 ) _1 U 9 22 
X [7/921 U 922 2 (-2 11 - 2 6 • 5 7/921 2 7t 9 22 - 3 2 ti 9 21 4 U 922 2 
+2 13 • 77/921 4 U922 4 + 2 12 7t 9 21 6 U922 5 + 2 3 • 3 2 7t 9 21 8 n 922 6 

+2 16 • 3 u 921 8 u 922 8 + 2 12 • 5 7/921 10 «922 9 + 2 6 • 11 U921 12 n 92 2 10 
+2 3 t/921 14 u 922 n ) - 2 (5 z)- 1 (2 2 • 3 - 2 12 • 3 2 u 92 i 2 u 922 3 

-2 5 • 3 2 • 7 7/921 4 «922 4 + 3 7/921 6 «922 5 + 2 15 • 3 • 5 7/921 6 U 922 7 

+2 10 • 17ti 9 21 8 n 9 22 8 + 2 17 • 19n 9 21 10 U922 U + 2 13 • 3 • 7t/ 9 21 12 U 922 12 

+2 7 • 23 U921 14 n 922 13 ) + 2 9 (5 z)~ 2 u 92l 3 u 922 A (-2 6 -3-5 + 3 u 9 2i 2 ^922 
+2 13 n 92 i 4 U922 4 + 2 14 • 5 n 9 2i 8 U922 8 + 2 8 • 11 u 92 i 10 u 92 2 9 ) 

-2 24 .7(5z)- 3 n 9 21 12 ^922 13 ] 
W 92 = M922 (2 2 + 2 5 7/921 4 7t 9 22 4 + W 92 l 6 ti922 5 - 2 8 (5z) _1 U 9 21 5 7i 9 22 5 ) 3 , 
W 92 = 2U 922 {2 2 + 2 b U 92 i i U 922 4 + U 921 & U 922 b -2 8 (5z)~ 1 7/ 9 2l 5 7t 9 22 5 ) 2 

X [-7/921 U 922 2 (-2 10 - 2 6 7/921 2 7t 92 2 + 2 12 • 7 7t 92 l 4 U 9 2 2 4 

+2 8 • 5 ti 9 2i 6 W922 5 + 2 2 • 3 u 921 8 u 922 6 + 2 15 • 3 n 92 i 8 ^922 8 

+2 10 • 77/ 9 21 10 7/ 9 22 9 + 2 5 • 5 7X921 12 U 922 10 + U 921 U U 922 U ) 

+(5 z)- 1 (2 2 • 3 - 2 12 • 3 2 n 92 i 2 u 922 3 - 2 5 • 3 3 u 92 i 4 7i 9 22 4 

+3 U 9 21 6 U 9 22 5 + 2 15 • 3 • 5 7i 9 21 6 U 922 7 + 2 10 • 11 U 921 8 U 922 8 

+2 17 • 19 7i 9 21 10 7/922 11 + 2 13 • 3 • 5 U 921 12 ^22^ + 2 7 • 11 7i 9 21 14 U 9 2 2 13 ) 

-2 8 (5 Z)~ 2 M921 3 U922 4 ("2 6 -3-5 + 3 7i 9 21 2 ^922 + 2 13 U 92 ^ U 922 4 

+2 14 • 5 t/921 8 n 9 22 8 + 2 11 n 92 i 10 n 92 2 9 ) + 2 23 • 7 (5 z)- 3 u 921 12 u 922 13 } 



Ew 92 = -2 1 U 92 1 1 [« 9 21 ("I + 2 9 U 9 21 2 n 9 22 3 + 2 4 ti 92 l 4 «922 4 

+2 11 u 9 2i 6 n 922 7 + 2 7 w 921 8 u 922 8 + 2 u 9 2i 10 u 922 9 ) 

-2 8 (5 z)- 1 (-1 + 2 4 u 9 2i 4 u 922 4 + 2 7 u 92 i 8 u 922 8 + 2 2 u 92 i 10 u 922 9 ) 

+2 17 (5z)- 2 U 92 l 9 U 9 22 9 ], 
£u> 92 = (5^) _1 M 9 21 _2 ^922 _1 

X [-2 6 - 3 U 9 21 2 U 9 22 + 2 9 U 9 21 4 U 9 22 4 + 2 4 U 9 21 6 W 9 22 5 
+2 11 U 9 2l 8 n 922 8 + 2 7 U 92 l 10 W 922 9 + 2 U 921 12 U 9 2 2 10 

-2 8 (5 z)- 1 u 921 u 922 (-3 + 2 4 u 9 2i 4 u 922 4 + 2 7 u 9 2i 8 u 922 8 
+2 2 n 92 i 10 n 922 9 ) + 2 17 (5 z)^ u 921 10 u 922 10 ]. 

The difference between the system and the fimiting system is given by 

■"921 -° «921 = -(5z) _1 M 92 l _1 n 92 2 _1 (128 - 1280 (5z) _1 U 9 21 U 9 22 + 5li 9 21 2 U 9 22) 

■u 9 2 2 -°n 9 22 = 2 (5z)" 1 u 9 2i~ 2 (64 - 640(5 2;)~ 1 u 9 2iti 9 22 + 3M 92 i 2 u 922 ). 

The equations -u 92 i = 0, 2 2 + 2 5 u 92 i 4 ii 922 4 + «92i 6 w 9 22 5 — 2 8 (5 z)^ 1 U921 5 U922 5 = 0, and w 9 22 = 
define L 9 , and and respectively. The proper transforms of the other lines on which the 
Painleve vector field is infinite are not visible in this chart. 
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Figure 4.1. Near the components of the infinity set 
Figure 4.2. The absolute value of the Weierstrass p function 
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FIGURE 5.1. The domain V, four half-lines Lf C V, and their distances to the origin 
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Figure 5.2. The region iZ^r, to the right of the curved boundary 



Figure 5.3. Truncated region in the x- plane, to the left of the curved boundary 
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Figure 5.4. Triply truncated region in the x-plane, to the left of the curved boundary 
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Figure 5.5. The first twenty of the asymptotic approximations in ( 5.68 1 of the poles of Pi(t) 
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